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I. INTRODUCTION 

The first demonstration of oscillation of a supercon- 
ducting qubit by Nakamura et al. in 1999i can be 
said to represent the " tip of the iceberg" : it rests on 
a huge volume of advanced research on Josephson junc- 
tions (JJ) and circuits developed during the last 25 
years. Some of this work has concerned fundamen- 
tal research on Josephson junctions and superconduct- 
ing quantum interferometers (SQUIDs) aimed at under- 
standing macroscopic quantum coherence (MQC) 2-3-4,5 , 
providing the foundation of the persistent current flux 
qubrfc§i£. However, there has also been intense research 
aimed at developing superconducting flux-based digital 
electronics and computers2*i2*ii. Moreover, in the 1990's 
the single-Cooper-pair box/transistor (SCB, SCT)A2, was 
developed experimentally and used to demonstrate the 
quantization of Cooper pairs on a small superconducting 
island 13 , which is the foundation of the charge qubilii* 1 ^. 

Since then there has been a steady 
developmenti^iSiiiiiLiS, with observation of microwave- 
induced Rabi oscillation of the two-level populations 
in chargo 21 . 22 ! 23 , and flux 24 i 25 i 26 i 27 qubits and dc- 
pulse driven oscillation of charge qubits with rf-SET 
detection 2 ^. An important step is the development of the 
charge-phase qubit, a hybrid version of the charge qubit 
consisting of an SCT in a superconducting loopSLSS, 
demonstrating Rabi oscillations with very long coherence 
time, of the order of 1 /us, allowing a large set of basic 
and advanced ("NMR-like") one-qubit operations (gates) 
to be performed^ 3 .. In addition, coherent oscillations 
have been demonstrated in the "simplest" J J qubits of 
them all, namely a single Josephson junctioriSS* 3 ^ 3 ^* 3 ^, 
or a two- J J dc-SQUIP 34 , where the qubit is formed by 
the two lowest states in the periodic potential of the JJ 
itself. 

Although a powerful J J-based quantum computer with 
hundreds of qubits remains a distant goal, systems with 
5-10 qubits will be built and tested by, say, 2010. Pair- 
wise coupling of qubits for two-qubit gate operations is 
then an essential task, and a few experiments with cou- 
pled JJ-qubits with fixed capacitive or inductive cou- 
plings have been reported 3 ^ 3 ^ 3 ^ 3 ^* 3 ^ 4 ^, in particular the 
first realization of a controlled-NOT gate with two cou- 
pled SCBaSi, used together with a one-qubit Hadamard 
gate to generate an entangled two-qubit state. 

For scalability, and simple operation, the ability to con- 
trol qubit couplings, e.g. switching them on and off, will 
be essential. So far, experiments on coupled JJ qubits 
have been performed without direct physical control of 
the qubit coupling, but there are many proposed schemes 
for two(multi)-qubit gates based on fixed or controllable 
physical qubit-qubit couplings or tunings of qubits and 
bus resonators. 

All of the JJ-circuit devices introduced above are 
based on nanoscale science and technology and represent 
emerging technologies for quantum engineering and, at 
best, information processing. One may debate the impor- 



tance of quantum computers on any time scale, but there 
is no doubt that the research will be a powerful driver of 
the development of solid-state quantum state engineering 
and quantum technology, e.g. performing measurements 
" at the edge of the impossible" . 

This article aims at describing the inner workings of 
superconducting Josephson junction (JJ) circuits, how 
these can form two-level systems acting as qubits, and 
how they can be coupled together to multi-qubit net- 
works. Since the field of experimental qubit applications 
is only five years old, it is not even clear if the field repre- 
sents an emerging technology for computers. Neverthe- 
less, the JJ-technology is presently the only example of a 
working solid state qubit with long coherence time, with 
demonstrated two-qubit gate operation and readout, and 
with potential for scalability. This makes it worthwile to 
describe this system in some detail. 

It needs to be said, however, that much of the basic 
theory for coupled JJ-qubits was worked out well ahead 
of experimentpi 4 * 41 ^, defining and elaborating basic op- 
eration and coupling schemes. We recommend the reader 
to take a good look at the excellent research and re- 
view paper by Makhlin et ali 4 ^ which describes the basic 
principles of a multi-JJ-qubit information processor, in- 
cluding essential schemes for qubit-qubit coupling. The 
ambition of the present article is to provide a both in- 
troductory and in-depth overview of essential Josephson 
junction quantum circuits, discuss basic issues of readout 
and measurement, and connect to the recent experimen- 
tal progress with JJ-based qubits for quantum informa- 
tion processing (JJ-QIP). 



II. NANOTECHNOLOGY, COMPUTERS AND 
QUBITS 

The scaling down of microelectronics into the nanome- 
ter range will inevitably make quantum effects like tun- 
neling and wave propagation important. This will even- 
tually impede the functioning of classical transistor com- 
ponents, but will also open up new opportunities for 
multi-terminal components and logic circuits built on e.g. 
resonant tunneling, ballistic transport, single electronics, 
etc. 

There are two main branches of fundamentally differ- 
ent computer architectures, namely logically irreversible 
and logically reversible. Ordinary computers are irre- 
versible because they dissipate both energy and infor- 
mation. Even if CMOS circuits in principle only draw 
current when switching, this is nevertheless the source 
of intense local heat generation, threatening to burn up 
future processor chips. The energy dissipation can in 
principle be reduced by going to single-electron devices, 
superconducting electronics and quantum devices, but 
this does not alter the fact that the information process- 
ing is logically irreversible. In the simplest case of an 
AND gate, two incoming bit lines only result in a single 
bit output, which means that one bit is in practice erased 



(initialized to zero) and the heat irreversibly dissipated 
to the environment. The use of quantum-effect devices 
does not change the fact that we are dealing with com- 
puters where each gate is logically irreversible and where 
discarded information constantly is erased and turned 
into heat. A computer with quantum device components 
therefore does not make a quantum computer. 

A quantum information processor has to be built on 
fundamentally reversible schemes with reversible gates 
where no information is discarded, and where all inter- 
nal processes in the components are elastic. This issue 
is connected with the problem of the minimum energy 
needed for performing a calculation^ 3 - (connected with 
the entropy change created by erasing the final result, 
i.e. reading the result and then clearing the register). A 
reversible information processor can in principle be built 
by classical means using adiabatically switched networks 
of different kinds. The principles were investigated in the 
1980's and form a background for much of the work on 
quantum computation^^^^i^ 8 ^ 9 *^. Recently there 
has been some very interesting development of reversible 
computers ( se e£iiSi£i£i and references therein). How- 
ever, a reversible computer still does not make a quan- 
tum computer. What is characteristic for a quantum 
computer is that it is reversible and quantum coherent, 
meaning that one can build entangled non-classical multi- 
qubit states. 

One can broadly distinguish between microscopic, 
mesoscopic and macroscopic qubits. Microscopic effec- 
tive two-level systems are localized systems confined by 
natural or artificial constraining potentials. Natural sys- 
tems then typically are atomic or molecular impurities 
utilizing electronic charge, or electronic or nuclear spins. 
These systems may be implanteo&Si or naturally oc- 
curring due to the material growth process^. A related 
type of impurity qubit involves endohedral fullerenes, i.e. 
atoms implanted into C60 or similar cages^, to be placed 
at specific positions on a surface prepared for control and 
readout. 

Mesoscopic qubit systems typically involve geomet- 
rically defined confining potentials like quantum dots. 
Quantum dots (QD) for qubitsS are usually made in 
semiconductor materials. One type of QD is a small 
natural or artificial semiconductor grain with quantized 
electronic levels. The electronic excitations may be ex- 
citonic (excitonsp*2i§i or biexcitonspS), charge-like^ 3 -, or 
spin-like (e.g. sing let-triplet )Si£&. Another type of QD 
is geometrically defined in semiconductor 2DEG by elec- 
trostatic split-gate arrangements. Although the host ma- 
terials are epitaxially layered semiconducting materials 
(e.g.GaAs/AlGaAs), the ungated 2DEG electronic sys- 
tem is metallic. A split-gate arrangement can then define 
a voltage-controlled system of quantum dots coupled to 
metallic reservoir electrodes, creating a system for elec- 
tron charge and spin transport through quantum point 
contacts and (effective) two-level quantum dot o 64 ' 65 ! 66 ' 67 . 

There is also a potential qubit based on liquid He 
superfluid technology, namely "Electrons on Helium", 
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FIG. 1: 3-junction persistent current flux qubit (PCQ) (inner 
loop) surrounded by a 2-junction SQUID. Courtesy of J.E. 
Moot], TU Delft. 




FIG. 2: Single Cooper pair box (SCB) (right) coupled to a 
single-electron transistor (SET) (left) for readout. Courtesy 
of P. Delsing, Chalmers. 

EoH 68 . This is really an atomic-like microscopic qubit: 
a thin film of liquid He is made to cover a Si surface, 
and electrons are bound by the image force above the He 
surface, forming an electronic two-level system. Qubits 
are laterally defined by electrostatic gate patterns in the 
Si substrate, which also defines circuits for qubit control, 
qubit coupling and readout. 

Macroscopic superconducting qubits - the subject of 
this article - are based on electrical circuits containing 
Josephson junctions (JJ). Looking at two extreme exam- 
ples, the principle is actually very simple. In one limit, 
the qubit is simply represented by the two rotation direc- 
tions of the persistent supercurrent of Cooper pairs in a 
superconducting ring containing Josephson tunnel junc- 
tions (rf-SQUID)(fiux qubit shown below in Fig. Q] 
for the Delft 3-junction flux qubit, 

In another limit, the qubit is represented by the pres- 
ence or absence of a Cooper pair on a small supercon- 
ducting island (Single Cooper pair Box, SCB, or Transis- 
tor, SCT) (charge qubit) 1 ! 13 ! 28 ! 69 , as illustrated in Fig. 
Hybrid circuitsSLiS*^ can in principle be tuned between 
these limits by varing the relations between the electro- 
static charging energy Ec and the Josephson tunneling 
energy Ej. 

All of these solid-state qubits have advantages and dis- 
advantages, and only systematic research and develop- 
ment of multi-qubit systems will show the practical and 
ultimate limitations of various systems. The supercon- 



ducting systems have presently the undisputable advan- 
tage of acutally existing, showing Rabi oscillations and 
responding to one- and two-qubit gate operations. In 
fact, even an elementary SCB two-qubit entangling gate 
creating Bell-type states has been demonstrated very 
recently^. All of the non-superconducting qubits are 
so far, promising but still potential qubits. Several of 
the impurity electron spin qubits show impressive relax- 
ation lifetimes in bulk measurements, but it remains to 
demonstrate how to read out individual qubit spins. 



III. BASICS OF QUANTUM COMPUTATION 
A. Conditions for quantum information processing 

DiVincenzo^ has formulated a set of rules and con- 
ditions that need to be fulfilled in order for quantum 
computing to be possible: 

1. Register of 2-level systems (qubits), n — 2 N states 
|101..01) (N qubits) 

2. Initialization of the qubit register: e.g. setting it to 
|000..00) 

3. Tools for manipulation: 1- and 2-qubit gates, e.g. 
Hadamard (H) gates to flip the spin to the equator, 
U H \0) = (|0) + |l))/2, and Controlled-NOT (CNOT) 
gates to create entangled states, UcnotUh\00 >— 
(| 00) + |H))/2 (Bell state) 

4. Read-out of single qubits \ip) = a\0) + 6e 40 |l) —> a, b 
(spin projection; phase (f> of qubit lost) 

5. Long decoherence times: > 10 4 2-qubit gate opera- 
tions needed for error correction to maintain coherence 
" forever" . 

6. Transport qubits and to transfer entanglement be- 
tween different coherent systems (quantum-quantum in- 
terfaces). 

7. Create classical-quantum interfaces for control, read- 
out and information storage. 

B. Qubits and entanglement 

A qubit is a two-level quantum system caracterized by 
the state vector 

|^> =cos||0>+sm|e^il) (3.1) 

Expressing |0) and |1) in terms of the eigenvectors of the 
Pauli matrix <r 2 , 

|o) = (J), l 1 ) = (?)- ( 3 - 2 ) 

this can be described as a rotation from the north pole 
of the |0) state, 




Li) 

FIG. 3: The Bloch sphere. Points on the sphere correspond 
to the quantum states \tp); in particular, the north and south 
poles correspond to the computational basis states |0) and 
|1); superposition cat-states = |0) + e^|l) are situated on 
the equator. 



can be characterised by a unit vector on the Bloch sphere: 

The state vector can be represented as a unitary vector 
on the Bloch sphere, and general unitary (rotation) oper- 
ations make it possible to reach every point on the Bloch 
sphere. The qubit is therefore an analogue object with 
a continuum of possible states. Only in the case of spin 
1/2 systems do we have a true two-level system. In the 
general case, the qubit is represented by the lowest levels 
of a multi-level system, which means that the length of 
the state vector may not be conserved due to transitions 
to other levels. The first condition will therefore be to 
operate the qubit so that it stays on the Bloch sphere 
(fidelity). Competing with normal operation, noise from 
the environment may cause fluctuation of both qubit am- 
plitude and phase, leading to relaxation and decoherence. 
It is a delicate matter to isolate the qubit from a perturb- 
ing environment, and desirable operation and unwanted 
perturbation (noise) easily go hand in hand. It is a major 
issue to design qubit control and read-out such that the 
necessary communication lines can be blocked when not 
in use. 

The state of N independent qubits can be represented 
as a product state, 

= |Vi>|^>...#w) = \^2-4n) (3.4) 

involving any one of all of the configurations |00...0 >, 
|00...1 >, 1 1 1 ... 1 >. A general state of an N-qubit 
memory register (i.e. a many-body system) can then 
be written as a time-dependent superposition of many- 
particle configurations 

|V(t)) = ci(t)|0...00) + c 2 (i)|0...01> (3.5) 
+ c 3 (t)|0.. .10) + ....+ c(t)|l...ll) 
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where the amplitudes Ciit) are complex, providing phase 
information. This state represents a time-dependent su- 
perposition of 27V N-body configurations which in gen- 
eral cannot be written as a product of one-qubit states 
and then represents an entangled (quantum correlated) 
many-body state. 

In the case of two qubits, the maximally entangled 
states are the so-called Bell states, 



|V) = (|00) + |ll))/2 (3.6) 

IV) = (|00) - |ll))/2 (3.7) 

|V^> = (|01) + |10))/2 (3.8) 

IV.) = (|01) - |10))/2 (3.9) 



where the last one is the singlet state. In the case 
of three qubits, the corresponding maximally entan- 
gled ("cat") states are the Greenberger-Horne-Zeilinger 
(GHZ) states^ 

|V) = (|000) ± |H1»/V8 (3.10) 

Another interesting entangled thrcc-qubit state appears 
in the teleportation process, 

|V) = [|00)(o|Q) + 6|1» + |01)(a|0) - 6|1)) (3.11) 
+|10)(6|0) + O |l)) + |ll)(6|0)- a |l))]/V8 

C. Operations and gates 

Quantum computation basically means allowing the 
iV-body state to develop in a fully coherent fashion 
through unitary transformations acting on all N qubits. 
The time evolution of the many-body system of N two- 
level subsystems can be described by the Schrodinger 
equation for the iV-level state vector \tp(t)), 

ind t \^)=H\ip). (3.12) 

in terms of the time-evolution operator characterizing by 
the time-dependent many-body Hamiltonian H (t) of the 
system determined by the external control operations and 
the perturbing noise from the environment, 

|V(t)> = U(t, t )m )). (3.13) 
The solution of Schrodinger equation for U(t, to) 

ihd t U{t, t ) = HU(t, t ) (3.14) 
may be written as 

[/(Mo) = 1 - t / H{t')U{t',to)dt' , (3.15) 
and finally, in terms of the time-ordering operator T, as 

-4- V H(t')dt' 

U{t,to)=Te o , (3.16) 



describing the time evolution of the entire TV-particle 
state in the interval [t, to]. If the total Hamiltonian com- 
mutes with itself at different times, the time ordering can 
be omitted, 

—r V H(t')dt' 

U(t,to)= e fiJ 'o . (3.17) 

This describes the time-evolution controlled by a ho- 
mogeneous time-dependent potential or electromagnetic 
field, e.g. dc or ac pulses with finite rise times but having 
no space-dependence. If the Hamiltonian is constant in 
the interval [to , t] , then the evolution operator takes the 
simple form 

[/(Mo) = e~r t H{t-t ) ) (318) 

describing the time-evolution controlled by square dc 
pulses. 

The time-development will depend on how many terms 
are switched on in the Hamiltonian during this time in- 
terval. In the ideal case, usually not realizable, all terms 
are switched off except for those selected for the specific 
computational step. A single qubit gate operation then 
involves turning on a particular term in the Hamiltonian 
for a specific qubit, while a two-qubit gate involves turn- 
ing on an interaction term between two specific qubits. 
In principle one can perform an 7V-qubit gate operation 
by turning on interactions for all N qubits. In practical 
cases, many terms in the Hamiltonian are turned on all 
the time, leading to a "background" time development 
that has to be taken into account. 

The basic model for a two-level qubit is the spin- 1/2 in 
a magnetic field. A system of interacting qubits can then 
be modelled by a collection of interacting spins, described 
by the Heisenberg Hamiltonian 

H(t) = Mt) S, + \ Jij{t) s * (3.19) 

controlled by a time-dependent external magnetic field 
hj(t) and by a time-dependent spin-spin coupling Jy(t). 

Expressing the Hamiltonian in Cartesian components, 
K(t) = {h x ,h y ,h z ), Si(t) = (S x ,S y ,S z ) and introducing 
the Pauli a-matrices, (S x , S y , S z ) — ^h((T x ,ay,a z ) we 
obtain a general N-qubit Hamiltonian with general qubit- 
qubit coupling: 

H = —- ^(e 4 cr zi -I- ReAicr xl + ImAicr yi ) (3.20) 

i 

+ ^ ^ t ^u,ij{t) <Jui<Jyj 

+ y](/t(*)o*j + 9xi(t)cr xl + g v i(t)a vi ) 

i 

We have here introduced time-independent components 
of the external field defining qubit energy level splittings 
£j, ReAi and TmAj along the z and x,y axes, as well 
as time-dependent components f%{t) and <?j(t) explicitly 
describing qubit operation and readout signals and noise. 
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Inserted into Ea. (|3.17|l . this Hamiltonian determines 
the time evolution of the many-qubit state. In the ideal 
case one can turn on and off each individual term of the 
Hamiltonian, including the two-body interaction, giving 
complete control of the evolution of the state. 

It has been shown that any unitary transformation can 
be achieved through a quantum network of sequential 
application of one- and two-qubit gates. Moreover, the 
size of the coherent workspace of the multi-qubit memory 
can be varied (in principle) by switching on and off qubit- 
qubit interactions. 

In the common case of NMR applied to 
molecules24iSi2&, one has no control of the fixed, 
direct spin-spin coupling. With an external magnetic 
field one can control the qubit Zeeman level splittings 
(no control of individual qubits). Individual qubits can 
be addressed by external RF-fields since the qubits have 
different resonance frequencies (due to different chemical 
environments in a molecule). Two-qubit coupling can 
be induced by simultaneous resonant excitation of two 
qubits. 

In the case of engineered solid state JJ-circuits, indi- 
vidual qubits can be addressed by local gates, control- 
ling the local electric or magnetic field. Extensive single- 
qubit operation has recently been demonstrated by Collin 
et ali^i Regarding two-qubit coupling, the field is just 
starting up, and different options are only beginning to 
be tested. The most straightforward approaches involve 
fixed capacitivei^Siiii or inductive2L2& coupling; such 
systems could be operated in NMR style, by detuning 
specific qubit pairs into resonance. Moreover, there are 
various solutions for controlling the direct physical qubit- 
qubit coupling strength, as will be described in Section 

El 



D. Readout and state preparation 

External perturbations, described by fi(t) and gi(t) in 
Ea. (|3.2U[> can influence the two- level system in typically 
two ways: (i) shifting the individual energy levels, which 
may change the transition energy and the phase of the 
qubit; and (ii) inducing transitions between the levels, 
changing the level populations. These effects arise both 
from desirable control operations and from unwanted 
noise (see^2i2L22i2S for a discussion of superconducting 
circuits, and Grangier et ali^S for a discussion of state 
preparation and quantum non-demolition (QND) mea- 
surements). 

To control a qubit register, the important thing is to 
control the decoherence during qubit operation and read- 
out, and in the memory state. In the qubit memory state, 
the qubit must be isolated from the environment. Op- 
eration and read-out devices should be decoupled from 
the qubit and, ideally, not cause any dcphasing or relax- 
ation. The resulting intrinsic qubit life times should be 
long compared to the duration of the calculation. In the 
operation and readout state, the qubit must be connected 



to the operation fields. This also opens up the system 
to a noisy environment, which puts great demands on 
signal-to-noise ratios. 

The readout operation is a particularly critical step. 
The ultimate purpose is to perform a " single-shot" quan- 
tum non-demolition (QND) measurement, determining 
the state (|0) or |1)) of the qubit in a single measure- 
ment and then leaving the qubit in that very state. Dur- 
ing the measurement time, the back-action noise from the 
"meter" will cause relaxation and mixing, changing the 
qubit state. The "meter" must then be sensitive enough 
to detect the qubit state on a time scale shorter than the 
induced relaxation time T* in the presence of the dissipa- 
tive back-action of the "meter" itself. Under these con- 
ditions it is possible to detect the qubit projection in a 
single measurement (single shot read-out). Performing a 
single-shot projective measurement in the qubit eigenba- 
sis then provides a QND measurement (note the " trivial" 
fact that the phase is irreversibly lost in the measurement 
process under any circumstances). 

The readout/measurement processes described above 
can be related to the Stern-Gerlach (SG) experiment 8 ', 
A SG spin filter acts as a beam splitter for flying qubits, 
creating separate paths for spin- up and spin-down atoms, 
preparing for the measurement by making it possible in 
principle to distinguish spatially between the two states 
of the qubit. The measurement is then performed by 
particle counters: a click in, say, the spin-up path col- 
lapses the atom to the spin- up state in a single shot. 
There is essentially no decohering back-action from the 
detector until the atom is detected. However, after de- 
tection of the spin-up atom the state is thoroughly de- 
stroyed: this qubit has not only decohered and relaxed, 
but is removed from the system. However, if it was entan- 
gled with other qubits, these are left in a specific eigen- 
state. For example, if the qubit was part of the Bell pair 
\tp) = (jOO) + |ll))/2, detection of spin up (|0» selects 
1 00) and leaves the other qubit in state |0). If instead the 
qubit was part of the Bell pair \ip) = (|01) + |10))/2, de- 
tection of spin up (|0)) selects 1 01} and leaves the other 
qubit in the spin-down state |1). Furthermore, for the 
3-qubit GHZ state \ip) = (|000) ± |111))/V8, detection 
of the first qubit in the spin up (|0)) state leaves the 
remaining 2-qubit system in the 1 00} state. 

Finally, in the case of the three-qubit entangled state 
in the teleportation process, Ea. (|3.11|l . detection of the 
first qubit in the spin up state |0), leaves the remaining 
2-qubit system in the entangled state [|0)(a|0) + 6|1)) + 
|l)(a|0) — 6|l))]/2. Moreover, detection of also the second 
qubit in the spin- up state |0) will leave the third qubit 
in the state (a|0) + b\l))/y/2. In this way, measurement 
can be used for state preparation. 

Applying this discussion to the measurement and read- 
out of JJ-qubit circuits, an obvious difference is that the 
JJ-qubits arc not flying particles. The detection can 
there not be turned on simply by the qubit flying into 
the detector. Instead, the detector is part of the JJ-qubit 
circuitry, and must be turned on to discriminate between 
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the two qubit states. These are not spatially separated 
and therefore sensitive to level mixing by detector noise 
with frequency around the qubit transition energy. The 
sensitivity of the detector determines the time scale T m 
of the measurement, and the detector-on back-action de- 
termines the time scale of qubit relaxation T£. Clearly 
T m <C Tj is needed for single-shot discrimination of |0) 
and |1). This requires a detector signal-to-noise (S/N) 
ratio 3> 1 , in which case one will have the possibility also 
to turn off the detector and leave the qubit in the deter- 
mined eigenstate, now relaxing on the much longer time 
scale T\ of the "isolated" qubit. This would then be the 
ultimate QND measurement. 

Note that in the above discussion, the qubit dephas- 
ing time ("isolated" qubit) does not enter because it is 
assumed to be much longer than the measurement time, 
T m <C T$. This condition is obviously essential for uti- 
lizing the measurement process for state preparation and 
error correction. 



IV. DYNAMICS OF TWO-LEVEL SYSTEMS 

To perform computational tasks one must be able to 
put a qubit in an arbitrary state. This is usually done 
in two steps. The first step, initialization, consists of 
relaxation of the initial qubit state to the equilibrium 
state due to interaction with environment. At low tem- 
perature, this state is close to the ground state. During 
the next step, time dependent dc- or rf-pulses are ap- 
plied to the controlling gates: electrostatic gate in the 
case of charge qubits, bias flux in the case of flux qubits, 
and bias current in the case of the J J qubit. Formally, 
the pulses enter as time-dependent contributions to the 
Hamiltonian and the state evolves under the action of the 
time-evolution operator. To study the dynamics of a sin- 
gle qubit two-level system we therefore first describe the 
two-level state, and then the evolution of this state under 
the influence of the control pulses (" perturbations" ) . 



and project onto the basis states 

ff£|m><m|V} = iW 

m 

obtaining the ususal matrix equation 

^ (k\H\m)a m = Ea k 

m 

where = (k\tp), and 



e A 

Hq P = -^( k Wz\m) - — (k\a x \m) 

giving the Hamiltonian matrix 



The Schrodinger equation is then given by 

1 ( e + 2E A \ / oi 

a-2 



(4.4) 



(4.5) 



(H-EM = 



A -e + 2E 



The eigenvalues are determined by 



det(H -E)=E 2 --{e 2 + A 2 ) = 



with the result 



Si,2 = ±iv / e 2 + A2 



The eigenvectors are given by: 



a 2 = — ai 



e + 2E 



(4.6) 

(4.7) 

= 
(4.8) 

(4.9) 
(4.10) 

(4.11) 



After normalisation 



A. The two-level state 

The general 1-qubit Hamiltonian has the form 

H = ~(e a z + A a x ) (4.1) 

The qubit eigenstates are to be found from the stationary 
Schrodinger equation 

HW)=E\$) (4.2) 

To solve the S-equation we expand the 1-qubit state in a 
complete basis, e.g. the basis states of the a z operator, 

|V> = |fc) =co|0)+ci|l) (4.3) 

k 



"-vf^sf-isfi^m <4 ' 12) 
02 = ±v /7^ =± _iy lT _L_ (4 . 13 ) 

We finally simplify the notation by fixing the signs of 
the amplitudes, 

fll = -7=2\Fw\ ; fl2 = : 7tf-W\ ; (4 ' 14) 

and explicitly writing down all the energy eigenstates, 

|£i) = oi|0)+o 2 |l) (4.15) 
|£ 2 } = a 2 |0}-oi|l) (4.16) 



where 

E 1 = -y$ + A* ; E 2 = +y$ + Al (4.17) 

The sign of \E 2 ) has been chosen to give the familiar 
expression for the superposition at the degeneracy point 
ei = where |ai| = |aa| = I/a/2, 

|^) = -L(|0} + |1)) (4.18) 
\E 2 ) = -L(|0) - |1» (4.19) 

B. The state evolution on the Bloch sphere 

To study the time-evolution of a general state, a conve- 
nient way is to expand in the basis of energy eigenstates, 

m)) = ci|B 1 )e- <Blt + c 2 \E 2 )e- iE * (4.20) 

If we know the coefficients at t = 0, then we know the 
time evolution. On the Bloch sphere this time evolution 
is represented by rotation of the Bloch vector with con- 
stant angular speed (E\ — Ei)/% around the direction 
defined by the energy eigenbasis. Indeed, by introduc- 
ing parameterization, C\ = cos 9', c 2 = sind'e 1 ^', we see 
that according to Eq. (|4.20f) the polar angle remains 
constant, 9' — const, while the azimuthal angle grows, 
4>'(t)4>'(0) + (Ei — E 2 )t/Ti. The primed angles here refer 
to a new coordinate system on the Bloch sphere related to 
the energy eigenbasis, which is obtained by rotation from 
the earlier introduced computational basis, Eq. f3.1[l . 

The dynamics on the Bloch sphere is conveniently de- 
scribed in terms of the density matrix for a pure quantum 
statefii, 

p=\1>)(il>\. (4.21) 

This is a 2 x 2 Hermitian matrix whose diagonal ele- 
ments p\ and p 2 define occupation probabilities of the 
basis states, hence satisfying the normalization condition 
Pi +P2 = 1, while the off-diagonal elements give informa- 
tion about the phase. The density matrix can be mapped 
on a real 3- vector by means of the standard expansion in 
terms of cr-matrices, 

p = + p x a x + p.ydy + p z a z ). (4.22) 

Direct calculation of the density matrix Ea. (|4.21|l using 
Ea. (|3.1() and comparing with Ea. H4.22|) shows that the 
vector p = (p x , p y , p z ) coincides with the Bloch vector, 

p = (sin 9 cos 4>, sin 9 sin </>, cos 9) (4.23) 

introduced in Fig. and also shown in Fig. 0J In the 

same a- matrix basis, the general two- level Hamiltonian 
takes the form 

H = (H x a x + Hy<j y + H z a z ). (4.24) 
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giving a 3- vector representation for the Hamiltonian, 

H = {H x ,H y ,H z ). (4.25) 

shown in Fig. 0] 




FIG. 4: The Bloch sphere: the Bloch vector p represents the 
states of the two-level system (same as in Fig. 3). The poles 
of the Bloch sphere correspond to the energy eigenstates; the 
vector H represents the two-level Hamiltonian. 

The time evolution of the density matrix is given by 
the Liouville equation, 

ihd t p=[H,p\. (4.26) 

The vector form of the Liouville equation is readily de- 
rived by inserting EQS. f4.24fl . I4.25| and using the com- 
mutation relations among the Pauli matrices, 

d t p=^[Hxp}. (4.27) 

This equation coincides with the Bloch equation for a 
magnetic moment evolving in a magnetic field, the role 
of the magnetic moment being played by the Bloch vector 
p which rotates around the effective "magnetic field" H 
associated with the Hamiltonian of the qubit (plus any 
driving fields) (Fig. 0}. 

C. dc-pulses, sudden switching and free precession 

To control the dynamics of the qubit system, one 
method is to apply dc (square) pulses which suddenly 
change the Hamiltonian and, consequently, the time- 
evolution operator. Sudden pulse switching means that 
the time-dependent Hamiltonian is changed so fast on 
the time scale of the evolution of the state vector that 
the state vector can be treated as time-independent - 
frozen - during the switching time interval. This implies 



that the system is excited by a Fourier spectrum with an 
upper cut-off given by the inverse of the switching time. 

In the specific scheme of sudden switching of different 
terms in the Hamiltonian using dc-pulses, the initial state 
is frozen during the switching event, and begins to evolve 
in time under the influence of the new 

|0) = |^(0)) = ci|£?i)+02|£a> (4.28) 

To find the coefficients we project onto the charge basis, 
k=0,l 

(0|0)=d(0|£ 1 )+c 2 <0|^) (4.29) 
(l|0)=d(l|S 1 )+ca<l|B 2 ) (4.30) 

and use the explict results for the energy eigenstates to 
calculate the matrix elements, obtaining 

1 = c\a\ + c 2 a 2 (4-31) 

= c\a 2 — c 2 a\ (4.32) 

As a result, 

|0> =a 1 \E 1 ) + a 2 \E 2 ) (4.33) 

This stationary state then develops in time governed by 
the constant Hamiltonian as 

\m) = aie~ 4filt |£i) + a 2 e~ iE ^\E 2 ) (4.34) 

Inserting the energy eigenstates we finally obtain the time 
evolution in the charge basis, 

\m) = 

|0) [a\e- iE ^ + a 2 2 e iElt ] + |1) [a^e-^ - e lBlt )] 

(4.35) 

The probability amplitudes of finding the system in one 
of the two charge states is then 

<0|V(*)> = [a\e~ lElt + a 2 2 e tElt ] (4.36) 
= [aia 2 {e- lElt - e lE ^)} (4.37) 

If the system is driven to the degeneracy point e\ = 0, 
where \a%\ = \a 2 \ = l/v*! then 

(0|VC*)> =<x>sEit (4.38) 
= sin-Eii (4.39) 

In particular, the probability of finding the system in 
state |1) (level 2) oscillates like 

p 2 (t) = =sin 2 E 1 t 

= - [1 - cos (E 2 - E^t] (4.40) 

with the frequency of the interlevel distance. On the 
Bloch sphere, this describes free precession around the 
X-axis. 
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c) d) 

FIG. 5: Qubit operations with dc-pulses: the vector H repre- 
sents the qubit Hamiltonian, and the vector p represents the 
qubit state, a) The qubit is initialized to the ground state; 
b) the Hamiltonian vector H is suddenly rotated towards x- 
axis, and the qubit state vector p starts to precess around 
H; c) when qubit vector reaches the south pole of the Bloch 
sphere, the Hamiltonian vector H is switched back to the 
initial position; the vector p remains at the south pole, indi- 
cating complete inversion of the level population (7r-pulse) ; d) 
if the Hamiltonian vector H is switched back when the qubit 
vector reaches the equator of the Bloch sphere (7r/2-pulse), 
then the p vector remains precessing at the equator, repre- 
senting equal-weighted superposition of the qubit states (cat 
states) \ip) — |0) + e l ^|l); this operation is the basis for the 
Hadamard gate. 



Let us consider, for example, the diagonal qubit Hamil- 
tonian, H = (e/2)a z , and apply a pulse 5e during a 
time r. This operation will shift phases of the qubit 
eigenstates by, ±Ser/2h. If the applied pulse is such 
that e is switched off, and instead, the a x component, 
A, is switched on, Fig EI then the state vector will ro- 
tate around the x-axis, and after the time At/21% = it 
(n- pulse) the ground state, |+), will flip and become, 
|+) — > | — ). This manipulation corresponds to the quan- 
tum NOT operation. Furthermore, if the pulse duration 
is twice smaller (7r/2-pulsc), then the ground state vector 
will approach the equator of the Bloch sphere and pre- 
cess along it after the end of the operation. Such state is 
an equal-weighted superposition of the basis states (cat 
state). 



D. Adiabatic switching 

Adiabatic switching represents the opposite limit to 
sudden switching, namely that the state develops so fast 
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on the time scale of the Hamiltonian that this can be 
regarded as "frozen" , i.e. the time-dependence of the 
Hamiltonian becomes parametric. This implies that en- 
ergy is conserved and no transitions are induced - the 
system stays in the same energy level (although the state 
changes) . 



E. Harmonic perturbation and Rabi oscillation 

A particularly interesting and practically important 
case concerns harmonic perturbation with small ampli- 
tude A and resonant frequency %uj = E 2 — E\. Let us 
consider the situation when the harmonic perturbation 
is added to the z-component of the Hamiltonian corre- 
sponding to a modulation of the qubit bias with a mi- 
crowave field. 

In the eigenbasis of the non-perturbed qubit, \Ei), 
l-EM, the Hamiltonian will take the form, 



H = E x a z + cosut (X z a z + X x a x ) , 



e A 
X Z = X— , A K = A— . 

-C/2 till 



(4.41) 



(4.42) 



The first pcrturbative term determines small periodic os- 
cillations of the qubit energy splitting, while the second 
term will induce inter level transitions. Despite the am- 
plitude of the perturbation being small, X/E2 <C 1, the 
system will be driven far away from the initial state be- 
cause of the resonance. Indeed, let us consider the wave 
function of the driven qubit on the form 

= a{t)e- lE ^ t ' h \E 1 ) + b{t)e- lE ^/ h \E 2 ) . (4.43) 

Substituting this ansatz into the Schrodinger equation, 

ih\i)) = H(t) (4.44) 

we get the following equations for the coefficients, 

iha = X x cosujte l(El ~ E2)t/h b, 

iU = X x cosujte i( - E2 - El)t/n a. (4.45) 

(Here we have neglected a small diagonal perturbation, 
A z .) 

Let us now focus on the slow evolution of the coeffi- 
cients on the time scale of qubit precession, and average 
Eas. H4.45|) over the period of the precession. This approx- 
imation is known in the theory of two-level systems as 
the "rotating wave approximation (RWA)" . Then, taking 
into account the resonance condition, we get the simple 
equations, 



iha = b, ihb-^- a, 

whose solutions read, 

oW(t) = b^(t) = e- iX ^ 2H , 



,(2 



>(t) = - b {2 \t) = e lX * 1 / 



2ti 



(4.46) 



(4.47) 



Thus, the dynamics of a driven qubit is characterized by 
a linear combination of the two wave functions, 

|-0 C1) > = -±= e - iX *V 2R (e-"* 1 */*!^) + e- iE *l h \E 2 )) , 

|^ (2) ) = -L e iA ^/ 2 ' t (e- iElt /' i |£i)-e- iE2 */ R | J B 2 )). 

(4.48) 

Let us assume that the qubit was initially in the ground 
state, I -Ei), and that the perturbation was switched on 
instantly. Then the wave function of the driven qubit 
will take the form, 

IV) = cos ^ e-W^i) + i sin ^ e - iE ^ H \E 2 ). 

(4.49) 

Correspondingly, the probabilities of the level occupa- 
tions will oscillate in time, 



P\ = cos 



X x t 1 



2 



2% 2 



1 + cos ■ 



X x t 



h 

2h 2 V h 



. 2 X x t 1 

Pi = Sin TT-r- = — I 1 — cos 



(4.50) 



with small frequency f2# = X x /h <C w, Rabi oscillations, 
illustrated in Fig. HJ 




FIG. 6: Rabi oscillation of populations of lower level (full 
line) and upper level (dashed line) at exact resonance (zero 
detuning). Tr = 2n/Qfi is the period of Rabi oscillations. 



F. Decoherence of qubit systems 

Descriptions of the qubit dynamics in terms of the den- 
sity matrix and Liouvillc equation are more general than 
descriptions in terms of the wave function, allowing the 
effects of dissipation to be included. The density ma- 
trix defined in Section IIV Bl for a pure quantum state 
possesses the projector operator property, p 2 = p. This 
assumption can be lifted, and then the density matrix 
describes a statistical mixture of pure states, say energy 
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eigenstates, 



= Y,Pi\Ei){Ei 



(4.51) 



The density matrix acquires off-diagonal elements when 
the basis rotates away from the energy eigenbasis. Such 
a mixed state cannot be represented by the vector on the 
Bloch sphere; however, its evolution is still described by 
the Louville equation Q4.26p . 



ihd t p = [H,p\. 



(4.52) 



The density matrix in Eq. I|4.51|) is the stationary solu- 
tion of the Liouville equation. The evolution of an arbi- 
trary density matrix, off-diagonal in the energy eigenba- 
sis, is given by the equations, 



Pi,p2 = const, p u cx e 



i{E 1 -E 2 )t/h 



(4.53) 



Dissipation is included in the density matrix descrip- 
tion by extending the qubit Hamiltonian and including 
interaction with an environment. The environment for 
macroscopic superconducting qubits basically consists of 
various dissipative elements in external circuits which 
provide bias, control, and measurement of the qubit. The 
" off-chip" parts of these circuits are usually kept at room 
temperature and produce significant noise. Examples are 
the fluctuations in the current source producing magnetic 
field to bias flux qubits and, similarly, fluctuations of the 
voltage source to bias gate of the charge qubits. Elec- 
tromagnetic radiation from the qubit during operation is 
another dissipative mechanisms. There are also intrinsic 
microscopic mechanisms of decoherence, such as fluctuat- 
ing trapped charges in the substrate of the charge qubits, 
and fluctuating trapped magnetic flux in the flux qubits, 
believed to produce dangerous 1/f noise. Another intrin- 
sic mechanism is possibly the losses in the tunnel junction 
dielectric layer. Various kinds of environment are com- 
monly modelled with an infinite set of linear oscillators 
in thermal equilibrium (thermal bath), linearly coupled 
to the qubit (Caldeira-Leggett mode&£). The extended 
qubit-plus-environment Hamiltonian has the form in the 
qubit energy eigenbasis^, 

H = -\^Ea z + ^(A^crj + \i±cr±)Xi 

i 

E + ) > (4 - 54) 

where E = E\ — E%, The physical effects of the two 
coupling terms in Eq. (|4.54() are quite different. The 
"transverse" coupling term proportional to Aj^ induces 
interlevel transitions and eventually leads to the relax- 
ation. The "longitudinal" coupling term proportional to 
A z commutes with the qubit Hamiltonian and thus does 
not induce interlevel transitions. However, it randomly 
changes the level spacing, which eventually leads to the 
loss of phase coherence, dephasing. The effect of both 



processes, relaxation and dephasing, are referred to as 
decoherence. 

Coupling to the environment leads, in the simplest 
case, to the following modification of the Liouville 
equatioi 



.83.84 



(oh 



dtPz = -7fr(Pz -P. 
J-l 



a 1 JT 1 

OtP\2 = Pl2 - 7fT Pl2 ■ 



(4.55) 



(4.56) 



This equation is known as the Bloch-Redfield equa- 
tion. The first equation describes relaxation of the 
level population to the equilibrium form, p^ — 
-(1/2; ) tanh(£y2fcT), T x being the relaxation time. The 
second equation describes disappearance of the off- 
diagonal matrix element during characteristic time T2, 
dephasing. 

The relaxation time is determined by the spectral den- 
sity of the environmental fluctuations at the qubit fre- 
quency, 



1 



A 



= E). 



(4.57) 



The particular form of the spectral density depends on 
the properties of the environment, which are frequently 
expressed via the impedance (response function) of the 
environment. The most common environment consists 
of a pure resistance, in this case, S^lo) cx lo, at low 
frequencies. 

The dephasing time consists of two parts, 



1 



1 

2T\ 



1 

Tl 



(4.58) 



The first part is generated by the relaxation process, 
while the second part results from the pure dephasing due 
to the longitudinal coupling to the environment. This 
pure dephasing part is proportional to the spectral den- 
sity of the fluctuation at zero frequency. 



(4.59) 



There is already a vast recent literature on de- 
coherence and noise in superconducting circuits, 
qubits and detectors, and how to engineer the 
qubits and environment to minimize decoherence and 

r: l x:x| .i c n 20.42. 86. 87.88. 89. 90. 91. 92. 93. 94. 95. 96. 97. 98. 99. 100. 101. 102. 103. 104. 105.1 

Many of these issues will be at the focus of this article. 



V. CLASSICAL SUPERCONDUCTING 
CIRCUITS 

In this section we describe a number of elementary su- 
perconducting circuits with tunnel Josephson junctions, 
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r 









FIG. 7: Current-biased Josephson junction (J J) (left), equiv- 
alent circuit (center), and effective (washboard- like) potential 
(right). The superconducting leads are indicated with dark 
color, and the tunnel junction with light color. 



which are used as building blocks in qubit applications. 
These basic circuits are: single current biased Joseph- 
son junction; single Josephson junction (JJ) included in 
a superconducting loop (rf SQUID); two Josephson junc- 
tions included in a superconducting loop (dc SQUID); 
and an ultra-small superconducting island connected to 
a massive superconducting electrode via tunnel Joseph- 
son junction (Single Cooper pair Box, SCB). 



A. Current biased Josephson junction 

The simplest superconducting circuit, shown in Fig. \7\ 
consists of a tunnel junction with superconducting elec- 
trodes, a tunnel Josephson junction, connected to a cur- 
rent source. An equivalent electrical circuit, which repre- 
sents the junction consists of the three lumped elements 
connected in parallel: the junction capacitance C, the 
junction resistance R, which generally differs from the 
normal junction resistance Rn and strongly depends on 
temperature and applied voltage, and the Josephson ele- 
ment associated with the tunneling through the junction. 

The current-voltage relations for the junction ca- 
pacitance and resistance have standard forms, Ic — 
C (dV/dt), and Ir = V/R. To write down a similar rela- 
tion for the Josephson element, it is necessary to intro- 
duce the superconducting phase difference <f>(t) across the 
junction, often simply referred to as the superconducting 
phase, which is related to the voltage drop across the 
junction, 



<K*) = f / vdt- 



where cj> is the time-independent part of the phase dif- 
ference. The phase difference can be also related to a 
magnetic flux, 



2e * „ $ 

= $ = 27T — , 

n $ 



(5.2) 



where $o = h/2e is the magnetic flux quantum. The 



current through the Josephson element has the for 



85 



1.1 



I r sin < 



(5.3) 



where /, 
the 



is the critical Josephson current, i.e. 
maximum non-dissipative current that may flow 
through the junction. The microscopic theory of 
superconductivityiii*ii2iii£ gives the following equation 
for the Josephson current, 



7tA A 

I r = tanh — , 

2eR N 2T' 



(5.4) 



where A is the superconducting order parameter, and T 
is the temperature. Using these relations and expressing 
voltage through the superconducting phase, we can write 
down Kirchhoff's rule for the circuit, 



2e 



Ccp- 



h 

2eR < ' 



I c sin < 



(5.5) 



where I e is the bias current. This equation describes the 
dynamics of the phase, and it has the form of a damped 
non-linear oscillator. The role of the non-linear induc- 
tance is here played by the Josephson element. 

The dissipation determines the qubit lifetime, and 
therefore circuits suitable for qubit applications must 
have extremely small dissipation. Let us assume zero 
level of the dissipation, dropping the resistive term in 
Eq. H5.5fl . Then the circuit dynamic equations, using the 
mechanical analogy, can be presented in the Lagrangian 
form, and, equivalently, in the Hamiltonian form. The 
circuit Lagrangian consists of the difference between the 
kinetic and potential energies, the electrostatic energy of 
the junction capacitors playing the role of kinetic energy 
while the energy of the Josephson current plays the role 
of potential energy. 

The kinetic energy corresponding to the first term in 
the Kirchhoff equation (|5.5|l reads, 



(5.6) 



This energy is equal to the electrostatic energy of the 
junction capacitor, CV 2 /2. It is convenient to introduce 
the charging energy of the junction capacitor charged 
with one electron pair (Cooper pair), 



Er 



[2ef 



(5.1) m which case Eq. (|5 



2C ' 
takes the form 



K{4>) 



4£ c 



(5.7) 



(5.8) 



The potential energy corresponds to the last two terms 
in Eq. I|5.5|) , and consists of the energy of the Josephson 
current, and the magnetic energy of the bias current, 



U(<j>) = Ej(l 



2e 



(5.9) 
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where Ej = h/2el c is the Josephson energy. This po- 
tential energy has a form of a washboard (see Fig. |7J). In 
the absence of bias current this potential corresponds to 
a pendulum with the frequency of small-amplitude oscil- 
lations given by 



2el c 

hC 



(5.10) 



This frequency is known as the plasma frequency of the 
Josephson junction. When current bias is applied, the 
pendulum potential becomes tilted, its minima becom- 
ing more shallow, and finally disappearing when the bias 
current becomes equal to the critical current, I e = Iq- 
At this point, the plasma oscillations become unstable, 
which physically corresponds to switching to the dissipa- 
tive regime and the voltage state. 

Now we are ready to write down the Lagrangian for 
the circuit, which is the difference between the kinetic 
and potential energies. Combining Eqs. (|5.8|l and l|5.9|l . 
we get, 



U 



AEr 



Ej(l 



2e 



(5.11) 




R 



FIG. 8: Superconducting quantum interference device - 
SQUID (left) consists of a superconducting loop (dark) in- 
terrupted by a tunnel junction (light); magnetic flux <E> is sent 
through the loop. Right: equivalent circuit. 




It is straightforward to check that the Kirchhoff equation 
(|5.5|) coincides with the dynamic equation following from 
the Lagrangian (|5.11|) . using 



FIG. 9: SQUID potential: the full (dark) curve corresponds 
to integer bias flux (in units of flux quanta), while the dashed 
(light) curve corresponds to half-integer bias flux. 



d_dL 
dt~d4> 



dL 



(5.12) 



It is important to emphasize, that the resistance of 
the junction can only be neglected for low temperatures, 
and also only for slow time evolution of the phase; both 
the temperature and the characteristic frequency must be 
small compared to the magnitude of the energy gap in the 
superconductor: T, Tim <C A. The physical reason behind 
this constraint concerns the amount of generated quasi- 
particle excitations in the system: if the constraint is 
fulfilled, the amount of equilibrium and non-equilibrium 
excitations will be exponentially small. Otherwise, the 
gap in the spectrum will not play any significant role, 
dissipation becomes large, and the advantage of the su- 
perconducting state compared to the normal conducting 
state will be lost. 



B. rf-SQUID 

The rf-SQUID is the next important superconducting 
circuit. It consists of a tunnel Josephson junction in- 
serted in a superconducting loop, as illustrated in Fig. |S] 
This circuit realizes magnetic flux bias for the Josephson 
junction^ 1 - 3 . To describe this circuit, we introduce the 
current associated with the inductance L of the leads, 



A(^_^ e ), </, e = ^$ e (5.13) 
leL n 



where ^ is the external magnetic flux threading the 
SQUID loop. The Kirchhoff rule for this circuit takes 
the form 



I r sin < 



0. (5.14) 



While neglecting the Josephson tunneling (I c = 0), this 
equation describes a damped linear oscillator of a con- 
ventional LC-circuit. The resonant frequency is then, 



1 



(5.15) 



and the (weak) damping is 7 = 1/ RC. 

In the absence of dissipation, it is straightforward to 
write down the Lagrangian of the rf-SQUID, 



AE C 



Ej(l 



Er 



0* 



(5.16) 



The last term in this equation corresponds to the energy 
of the persistent current circulating in the loop, 



Er 



4ir 2 L' 



(5.17) 



The potential energy U ((f)) corresponding to the last two 
terms in Eq. (|5.16l) is schetched in Fig. 

For bias flux equal to integer number of flux quanta, 
or (j> e 2irn, the potential energy of the SQUID has one ab- 
solute minimum at <f> = <f> e . For half integer flux quanta 
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FIG. 10: dc SQUID consists of two tunnel junctions included 
in a superconducting loop; arrows indicate the direction of 
the positive Josephson current. 




FIG. 11: Single Cooper pair box (SCB): a small supercon- 
ducting island connected to a bulk superconductor via a tun- 
nel junction; the island potential is controlled by the gate 
voltage V g . 



the potential energy has two degenerate minima, which 
correspond to the two persistent current states circulat- 
ing in the SQUID loop in the opposite directions. This 
configuration of the potential energy provides the basis 
for constructing a persistent-current flux qubit (PCQ). 

C. dc SQUID 

We now consider consider the circuit shown in Fig. I1UI 
consisting of two Josephson junctions coupled in paral- 
lel to a current source. The new physical feature here, 
compared to a single current-biased junction, is the de- 
pendence of the effective Josephson energy of the double 
junction on the magnetic flux threading the SQUID loop. 

Let us evaluate the effective Josephson energy. Now 
the circuit has two dynamical variables, superconducting 
phases, 4>\^ across the two Josephson junctions. Defin- 
ing phases as shown in the figure, and applying consider- 
ations from Sections IV Al and IV Bl we find for the static 
Josephson currents, 

I c i sin 01 - I c2 sin 02 = h, (5.18) 

where I e is the biasing current. Let us further assume 
small inductance of the SQUID loop and neglect the mag- 
netic energy of circulating currents. Then the total volt- 
age drop over the two junctions is zero, V\ + Vi = 0, and 
therefore, cftx + <f>2 = (f> e , where <f> e is the biasing phase 
related to biasing magnetic flux. Introducing new vari- 
ables, 




(5.19) 



and taking into account that 2<f> + — cf> e , we rewrite equa- 
tion l|5.18|l on the form, 

7 c (0 e )sin(0_ +a) =I e , (5.20) 

where 

Ic{<t>e) = +Ic2 + Wc2COS0 e , (5.21) 



and 

, Id — Ic2 , 4>e w , 

tana = —tan — . (5.22) 

Id + Ic2 2 

For a symmetric SQUID with I c \ = I C 2, giving a = 0, 
Eq. (|5.2U|) reduces to the form 

2I C cos(0 e /2) sin 0_ = I e . (5.23) 

The potential energy generated by Eq. I5.20fl has the 
form, 

U ( <f>) =2Ejcos(^j(l- cos 0_ ) - ^I e 0_ , (5.24) 

which indeed is similar to the potential energy of a sin- 
gle current biased junction, Eq. (|5.9(l and Fig. 03 but 
with flux-controlled critical current. This property of 
the SQUID is used in qubit applications for controlling 
the Josephson coupling, and also for measuring the qubit 
flux. 

The kinetic energy of the SQUID can readily be writ- 
ten down noticing that it is associated with the charging 
energy of the two junction capacitances connected in par- 
allel, 

K{ ® = (£) 2(Cl + C2) y- (5 - 25) 

Thus the Lagrangian for the SQUID has a form similar 
to Eq. (533) where E c = (2e) 2 /2(Ci + C 2 ), 

, * . h 2 6 2 ( (b e \ , . n 

L(<f>, 4>) = ^ - 2Ej cos If (1 - cos 0) + -J e 0. 

(5.26) 

D. Single Cooper Pair Box (SCB) 

There is a particularly important Josephson junction 
circuit consisting of a small superconducting island con- 
nected via a Josephson tunnel junction to a large super- 
conducting reservoir (see Fig. Ill|) . 
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FIG. 12: Single Cooper pair transistor (SCT): SCB with loop- 
shape bulk electrodes; charge fluctuations on the island pro- 
duces current fluctuation in the loop. 



The island is capacitively coupled to another massive 
electrode, which may act as an electrostatic gate. The 
voltage source V g controls the gate potential. In the nor- 
mal state, such a circuit is named a Single Electron Box 
(SEB)iii for the following reason: if the junction resis- 
tance exceeds the quantum resistance R q « 26 kfl, and 
the temperature is small compared to the charging en- 
ergy of the island, the system is in a Coulomb blockade 
regimeiii^ where the electrons can only be transferred 
to the island one by one, the number of electrons on 
the island being controlled by the gate voltage. In the 
superconducting state, the same circuit is called a Sin- 
gle Cooper pair Box (SCB) 116117 ; for a review see the 
bookie. An experimental SCB device is shown in Fig. 
121 In this section we consider a classical Lagrangian for 
the circuit. Since the structure now has two capacitances, 
one from the tunnel junction, C, and another one from 



the gate, C g , the electrostatic term in the Hamiltonian 
must be reconsidered. 

Let us first evaluate the electrostatic energy of the 
SCB. It has the form, 

CV 2 , C g {V g ~Vf 



(5.27) 

where V is the voltage over the tunnel junction. Then the 
Lagrangian can be written (omitting the constant term) , 



2 



h 

Ye' 



Cg 



Vg) -£./(!- cos 0), (5.28) 



where C s = C + C g . 

The interferometer effect of two Josephson junctions 
connected in parallel (Fig. I12JI can be used to control the 
Josephson energy of the single Cooper pair box. This 
setup can be viewed as a flux-biased dc SQUID where 
Josephson junctions have very small capacitances and are 
placed very close to each other so that the island confined 
between them has large charging energy. The gate elec- 
trode is connected to the island to control the charge. 



The effect of the gate electrode is only essential for the 
kinetic term, and repeating previous analysis we arrive at 
the following equation for the kinetic energy (assuming 
for simplicity identical junctions), 



K{ 



2 



-0- " ^v g 
C s = 2C + C q 



(5.29) 



Combining this kinetic energy with the potential en- 
ergy derived in previous subsection, we arrive at the La- 
grangian of the SCB (cf. Eq. (|5.28l) ) where both the 
charging energy and the Josephson energy can be con- 
trolled, 



2 l2e 



2Ej cos 



COS( 



(5.30) 



VI. 



QUANTUM SUPERCONDUCTING 
CIRCUITS 



One may look upon the Kirchhoff rules, as well as the 
circuit Lagrangians, as the equations describing the dy- 
namics of electromagnetic field in the presence of the 
electric current. Generally, this electromagnetic field is a 
quantum object, and therefore there must be a quantum 
generalization of the equations in the previous section. 
At first glance, the idea of quantization of an equation de- 
scribing a macroscopic circuit containing a huge amount 
of electrons may seem absurd. To convince ourselves that 
the idea is reasonable, it is useful to recall an early ar- 
gument in favor of the quantization of electron dynamics 
in atoms, and to apply it to the simplest circuit, an rf- 
SQUID: When the current oscillations are excited in the 
SQUID, it works as an antenna radiating electromagnetic 
waves. Since EM waves are quantized, the same should 
apply to the antenna dynamics. 

To quantize the circuit equation, we follow the conven- 
tional way of canonical quantization: first we introduce 
the Hamiltonian and then change the classical momen- 
tum to the momentum operator. The Hamiltonian is 
related to the Lagrangian as 



H(p, <j>)=p4>-L, 



(6.1) 



where p is the canonical momentum conjugated to coor- 
dinate d>. 



P : 



dL 



(6.2) 



For the simplest case of a single junction, Eq. (|5.6|l , the 
momentum reads, 



5.1 <* 



(6.3) 
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The so defined momentum has a simple interpretation: 
it is proportional to the charge q = CV on the junction 
capacitor, p = (h/2e)q, or the number n of electronic 
pairs on the junction capacitor, 



P 



fin. 



(6.4) 



The Hamiltonian for the current-biased junction has 
the form (omitting a constant), 



H = E c n z -Ej COS0-— I, 

2e 



(6.5) 



Similarly, the Hamiltonian for the SQUID circuit has the 
form 



H(n, 4>) = E c n 2 - Ej cos 4> + El 



(6.6) 



The dc SQUID considered in the previous Section IV CI 
has two degrees of freedom. The Hamiltonian can be 
written by generalizing Eqs. I|6.5|) . (|6.6|) for the phases 
<t>±. In the symmetric case we have, 



H = E c n\ + E c n 2 _ - 2Ej cos < 



cos < 



E L 



(20+ - <t>e? 



2e 



For the SCB, Eq. 
the form, 



P 



, the conjugated momentum has 



2e 



h X C "v 



and the Hamiltonian reads, 
H = E c (n - r. 



E j cos < 



(6.8) 



(6.9) 



where now Ec — (2e) 2 /2Cs, and n = p/h has the mean- 
ing of the number of electron pairs (Cooper pairs) on the 
island electrode. n g — ~C g V g /2e is the charge on the 
gate capacitor (in units of Cooper pairs (2e)), which can 
be tuned by the gate potential, and which therefore plays 
the role of external controlling parameter. 

The quantum Hamiltonian results from Eq. (|6.2|l by 
substituting the classical momentum p for the differential 
operator, 



-ih 
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Similarly, one can define the charge operator, 
and the operator of the pair number, 



(6.10) 



(6.11) 



(6.12) 



The commutation relation between the phase operator 
and the pair number operator has a particularly simple 
form, 

[(f), n] = i. (6.13) 

The meaning of the quantization procedure is the fol- 
lowing: the phase and charge dynamical variables can 
not be exactly determined by means of physical measure- 
ments; they are fundamentally random variables with the 
probability of realization of certain values given by the 
modulus square of the wave function of a particular state, 



(0) = / (q) = J PMMWd*. 

(6-14) 

The time evolution of the wave function is given by the 
Schrodinger equation, 



(6.15) 



where H is the circuit quantum Hamiltonian. The ex- 
plicit form of the quantum Hamiltonian for the circuits 
considered above, is the following: 



(6.7) rf-SQUID: 



H = E c h — Ej cos (f> + E L K 7 ' 



Current biased J J: 



H = Ecn — Ej cos 4> H I e 4>\ 

2e 

dc-SQUID: 

H = E c n\ + E c n 2 _ - 2Ej cos 4> + cos < 

2 2e v 

and finally the single Cooper pair box (SCB): 



(6.16) 



(6.17) 



H = Ec(n — n g ) 2 — Ej cos( 



(6.18) 



(6.19) 



For junctions connecting macroscopically large elec- 
trodes, the charge on the junction capacitor is a con- 
tinuous variable. This implies that no specific boundary 
conditions on the wave function are imposed. The sit- 
uation is different for the SCB: in this case one of the 
electrodes, the island, is supposed to be small enough 
to show pronounced charging effects. If tunneling is for- 
bidden, electrons are trapped on the island, and their 
number is always integer (the charge quantization con- 
dition). However, there is a difference between the ener- 
gies of even and odd numbers of electrons on the island: 
while an electron pair belongs to the superconducting 
condensate and has the additional energy Ec, a single 
electron forms an excitation and thus its energy consists 
of the charging energy, Ec/2 plus the excitation energy, 
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A (parity effect To prevent the appearance of indi- 
vidual electrons on the island and to provide the SCB 
regime, the condition A ^> Ec/2 must be fulfilled. Thus 
when the tunneling is switched on, only Josephson tun- 
neling is allowed since it transfers Cooper pairs and the 
number of electrons on the island must change pairwise, 
?i_ = integer. In order to provide such a constraint, peri- 
odic boundary conditions on the SCB wave function are 
imposed, 

ip(<j>) =V# + 2tt). (6.20) 

This implies that arbitrary state of the SCB is a super- 
position of the charge states with integer amount of the 
Cooper pairs, 

V>(0)=J>„e in *. (6.21) 

n 

The uncertainty of the dynamical variables are not im- 
portant as long as the relative mean deviations of dy- 
namical variables are small, i.e. the amplitudes of the 
quantum fluctuations are small. In this case, the particle 
behaves as a classical particle. It is known from quan- 
tum mechanics, that this corresponds to large mass of 
the particle, in our case, to large junction capacitance. 
Thus we conclude that the quantum effects in the circuit 
dynamics are essential when the junction capacitances 
are sufficiently small. The qualitative criterion is that 
the charging energy must be larger than, or comparable 
to, the Josephson energy of the junction, Ec ~ Ej. 

Typical Josephson energies of the tunnel junctions in 
qubit circuits are of the order of a few degrees Kelvin 
or less. Bearing in mind that the insulating layers of 
the junctions have the thickness of few atomic distances, 
and modeling the junction as a planar capacitor, the esti- 
mated junction area should be smaller than a few square 
micrometers to observe the circuit quantum dynamics. 

One of most important consequences of the quantum 
dynamics is quantization of the energy of the circuit. Let 
us consider, for example, the LC circuit. In the classical 
case, the amplitude of the plasma oscillations has con- 
tinuous values. In the quantum case the amplitude of 
oscillation can only have certain discrete values defined 
through the energy spectrum of the oscillator. The linear 
oscillator is well studied in the quantum mechanics, and 
its energy spectrum is very well known, 

E n = hu LC {n + 1/2), n = 0,l,2,... (6.22) 

One may ask, why is quantum dynamics never ob- 
served in ordinary electrical circuits? After all, in high- 
frequency applications, frequencies up to THz are avail- 
able, which corresponds to a distance between the quan- 
tized oscillator levels of order 10K, which can be observed 
at sufficiently low temperature. For an illuminative dis- 
cussion on this issue see the paper by Martinis, Devoret 
and Clarke^. According to Ref. it is the dissipation 
that kills quantum fluctuations: as known from classi- 
cal mechanics, the dissipation (normal resistance) broad- 
ens the resonance, and good resonators must have small 



resonance width, 7 = 1/RC compared to the resonance 
frequency, 7 <C l^lc- It is intuitively clear that the quan- 
tization effect will be destroyed when the level broaden- 
ing exceeds the level spacing. For quantum behavior of 
the circuit, narrow resonances, 7 -C ujlCi are therefore 
essential. However, even in this case it is hard to ob- 
serve the quantum dynamics in linear circuits such as 
LC-resonatorsiS^ because the expectation values of the 
linear oscillator follow the classical time evolution. Thus 
the presence of non-linear circuit elements is essential. 

The linear oscillator provides the simplest example of 
a quantum energy level spectrum: it only consists of dis- 
crete levels with equal distance between the levels. The 
level spectrum of the Josephson junction associated with 
a pendulum potential is more complicated. Firstly, be- 
cause of the non-linearity (non-parabolic potential wells), 
the energy spectrum is non-equidistant (anharmonic), 
the high-energy levels being closer to each other than the 
low-energy levels. Moreover, for energies larger than the 
amplitude of the potential (top of the barrier), Ej, the 
spectrum is continuous. Secondly, one has to take into 
account the possibility of a particle tunneling between 
neighboring potential wells: this will produce broadening 
of the energy levels into energy bands. The level broad- 
ening is determined by the overlap of the wave function 
tails under the potential barriers, and it must be small 
for levels lying very close to the bottom of the potential 
wells. Such a situation may only exist if the level spac- 
ing, given by the plasma frequency of the tunnel junction 
is much smaller than the Josephson energy, tuvj <C Ej, 
i.e. when Ec <C Ej. This almost classical regime with 
Josephson tunneling dominating over charging effects, is 
called the phase regime, because phase fluctuations are 
small and the superconducting phase is well defined. In 
the opposite case, Ec ^S> Ej, the lowest energy level lies 
well above the potential barrier, and this situation corre- 
sponds to wide energy bands separated by small energy 
gaps. In this case, the wave function far from the gap 
edges can be well approximated with a plane wave, 

^,(0=exp^. (6.23) 

This wave function corresponds to an eigenstate of the 
charge operator with well defined value of the charge, q. 
Such a regime with small charge fluctuations is called the 
charge regime. 

VII. BASIC QUBITS 

The quantum superconducting circuits considered 
above contain a large number of energy levels, while for 
qubit operation only two levels are required. Moreover, 
these two qubit levels must be well decoupled from the 
other levels in the sense that transitions between qubit 
levels and the environment must be much less probable 
than the transitions between the qubit levels itself. Typi- 
cally that means that the qubit should involve a low-lying 
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pair of levels, well separated from the spectrum of higher 
levels, and not being close to resonance with any other 
transitions. 



A. Single Josephson Junction (JJ) qubit 

The simplest qubit realization is a current biased JJ 
with large Josephson energy compared to the charging 
energy. In the classical regime, the particle representing 
the phase either rests at the bottom of one of the wells of 
the "washboard" potential (Fig. 0), or oscillates within 
the well. Due to the periodic motion, the average voltage 

across the junction is zero, = 0. Strongly excited states, 
where the particle may escape from the well, correspond 
to the dissipative regime with non-zero average voltage 

across the junction, 7^ 0. 

In the quantum regime described by the Hamiltonian 

H = E c n 2 -Ejcos(b-^-I e (j), (7.1) 

2e 

particle confinement, rigorously speaking, is impossi- 
ble because of macroscopic quantum tunneling (MQT) 
through the potential barrier, see Fig. ^| However, 
the probability of MQT is small and the tunneling may 
be neglected if the particle energy is close to the bot- 
tom of the local potential well, i.e. when E <C Ej. To 
find the conditions for such a regime, it is convenient 
to approximate the potential with a parabolic function, 
U{<t>) ~ (l/2)Ej cos 0o (0— 0o) 2 , where 0o corresponds to 
the potential minimum, Ej sin 0o — (H/2e)I e . Then the 
lowest energy levels, Ek — fiw p (k + 1/2) are determined 
by the plasma frequency, lo p — 2 1 / 4 oj,/(1 — Ie/Ic) 1 ^ 4 . It 
then follows that the levels are close to the bottom of the 
potential if Ec <C Ej, i.e. when the Josephson junction 
is in the phase regime, and moreover, if the bias current 
is not too close to the critical value, I e < I c . 

It is essential for qubit operation that the spectrum in 
the well is not equidistant. Then the two lowest energy 
levels, k = 0, l can be employed for the qubit operation. 
Truncating the full Hilbert space of the junction to the 
subspace spanned by these two states, |0) and |1), we 
may write the qubit Hamiltonian on the form, 

ff, = -~«7„ (7.2) 

where e = E\ — Eo . 

The interlevel distance is controlled by the bias cur- 
rent. When bias current approaches the critical cur- 
rent, level broadening due to MQT starts to play a role, 
E k -> E k +iT k /2. The MQT rate for the lowest level is 
given by22 
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FIG. 13: Quantized energy levels in the potential of a current 
biased Josephson Junction 

where U max = 2v / 2($o/27r)(l - I e /I c ) 3 / 2 is the height of 
the potential barrier at given bias current. 

B. Charge qubits 

1. Single Cooper pair Box - SCB 

An elementary charge qubit can be made with the SCB 
operating in the charge regime, Eq 3> Ej. Neglecting 
the Josephson coupling implies the complete isolation of 
the island of the SCB, with a specific number of Cooper 
pairs trapped on the island. Correspondingly, the eigen- 
functions, 

E c {n-n g ) 2 \n) = E n \n), (7.4) 

correspond to the charge states n = 0, 1, 2..., with the en- 
ergy spectrum E n — Ec(n — n g ) 2 , as shown in Fig. [21 
The ground state energy oscillates with the gate voltage, 
and the number of Cooper pairs in the ground state in- 
creases. There are, however, specific values of the gate 
voltage, e.g. n g = 1/2 where the charge states |0) and 
|1) become degenerate. Switching on a small Josephson 
coupling will then lift the degeneracy, forming a tight 
two-level system. 

The qubit Hamiltonian is derived by projecting the 
full Hamiltonian 1)6.1 9[1 on the two charge states, |0), |1), 
leading to 

Hscb = ™(e<r« + Aa x ), (7.5) 

where e = Ec{l — 2n 9 ), and A = Ej. The qubit level 
energies are then given by the equation 

E^= T \^E 2 {l-2n g y+E 2 , (7.6) 

the interlevel distance being controlled by the gate volt- 
age. At the degeneracy point, n g — 1/2, the diag- 
onal part of the qubit Hamiltonian vanishes, the lev- 
els being separated by the Josephson energy, Ej, and 
the qubit eigenstates corresponding to the cat states, 
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n= -1 1 



FIG. 14: Single Cooper pair box (SCB): charging energy (up- 
per panel), and charge on the island (lower panel) vs gate 
potential. Washed-out onsets of the charge steps indicate 
quantum fluctuations of the charge on the island. 




FIG. 15: Energy spectrum of the SCB (solid lines): it results 
from hybridization of the charge states (dashed lines). 



\Ei),\E 2 ) — |0) =F For theses states, the average 
charge on the island is zero, while it changes to =F2e far 
from the degeneracy point, where the qubit eigenstates 
approach pure charge states. 

The SCB was first experimentally realized by Lafarge 
et aliiii, observing the Coulomb staircase with steps of 
2e and the superposition of the charge states, see alsoii 
Realization of the first charge qubit by manipulation of 
the SCB and observation of Rabi oscillations was done 
by Nakamura et aliii2ii2£, and further investigated the- 
oretically by Choi et aliiSi. 



current is distributed between the two junctions. The 
answer to this question is apparently equivalent to eval- 
uating the persistent current circulating in the SQUID 
loop. This current was neglected so far because the as- 
sociated induced flux <fi = 2<p + — (j) e was assumed to be 
frozen, cf> — 0, in the limit of infinitely small SQUID in- 
ductance, L = 0. Let us now lift this assumption and 
allow fluctuation of the induced flux; the Hamiltonian 
will then take the form l|tj.7fl , in which a gate potential is 
included in the charging term of the SCB, and the term 
containing external current is dropped, 



Hsct = Ec(n. 
-2Ej cos (f>+ cos + 



E, 



n g )U 
(20+ 



E c n\ 



(7.7) 



(?i + = —id/d(j)-\-). For small but non-zero inductance, 
the amplitude of the induced phase is small, 4> = 2</>+ — 
e <C 1, and the cosine term containing <j) + can be ex- 
panded, yielding the equation 



H, 



SCT 



= H 



SCB( 



H osc (4>) + Hi, 



(7.; 



Hscb(<I>-) is the SCB Hamiltonian I|rj.l9f) with the flux 
dependent Josephson energy, Ej((f> e ) = 2Ej cos(</> e /2). 
Hosc{4>) describes the linear oscillator associated with the 
variable d, 



H osc (4>) = AE c n 2 + E L - 



and the interaction term reads, 



Hint = EjSVA \ — 



cos 



(7.9) 



(7.10) 



Thus, the circuit consists of the non-linear oscillator of 
the SCB linearly coupled to the linear oscillator of the 
SQUID loop. This coupling gives the possibility to mea- 
sure the charge state of the SCB by measuring the per- 
sistent currents and the induced flux. 

Truncating Eq. 1)7. 8|) we finally arrive at the Hamil- 
tonian which is formally equivalent to the spin-oscillator 
Hamiltonian. 



Hsct — — - (eer 2 

In this equation, 
Ej sin(0 e /2). 



)a x ) + \4>a x + H 



(7.11) 



= 2Ej cos(0 e /2), and A = 



2. Single Cooper pair Transistor - SCT 

In the SCB, charge fluctuations on the island gener- 
ate fluctuating current between the island and large elec- 
trode. In the two-junction setup discused at the end of 
Section IV Dl an interesting question concerns how the 



C. Flux qubit 

1. Quantum rf-SQUID 

An elementary flux qubit can be constructed from an 
rf-SQUID operating in the phase regime, Ej ^> Eq. Let 
us consider the Hamiltonian (|6.tj|l at <f> e = 7r, i.e. at half 
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FIG. 16: Double-well potential of the rf-SQUID with de- 
generate quantum levels in the wells (black). Macroscopic 
quantum tunneling (MQT) through the potential barrier in- 
troduces a level splitting A, and the lowest level pair forms a 
qubit. 




o 

FIG. 17: Flux qubit: truncation of the junction Hamiltonian; 
dashed lines indicate potentials of the left and right wells with 
ground energy levels. 

integer bias magnetic flux. The potential, [/(</>), shown 
in Fig. 1 161 has two identical wells with equal energy levels 
when MQT between the wells is neglected (phase regime, 
uij <C Ej). These levels are connected with current fluc- 
tuations within each well around averaged values cor- 
responding to clockwise and counterclockwise persistent 
currents circulating in the loop (the flux states). Let 
us consider the lowest, doubly degenerate, energy level. 
When the tunneling is switched on, the levels split, and 
a tight two-level system is formed with the level spacing 
determined by the MQT rate, which is much smaller than 
the level spacing in the well. In the case that the tunnel- 
ing barrier is much smaller than the Josephson energy, 
the potential in Eq. I|5.9[) can be approximated, 

U(</>) = Ej(1-cob<I>)+E l Q^L 

« i^- e f -/^+i±fA (7.12) 

where tf> = <f> — tt, f = 4> e — it, and where 

e = |^ - 1 « 1. (7.13) 

determines the height of the tunnel barrier. 

The qubit Hamiltonian is derived by projecting the 
whole Hilbert space of the full Hamiltonian l|6.6|) on the 



subspace spanned by these two levels. The starting point 
of the truncation procedure is to approximate the double 
well potential with Ui and U r , as shown in Fig. El to 
confine the particle to the left or to the right well, respec- 
tively. The corresponding ground state wave functions |Z) 
and \r) satisfy the stationary Schrodingcr equation, 

Hi\l)=E t \l), H r \r)=E r \r). (7.14) 

The averaged induced flux for these states, and 4> r 
have opposite signs, manifesting opposite directions of 
the circulating persistent currents. Let us allow the bias 
flux to deviate slightly from the half integer value, <f> e = 
7r + /, so that the ground state energies are not equal 
but still close to each other, -E; ~ E r . The tunneling will 
hybridize the levels, and we can approximate the true 
eigenf unction, \E), 

H\E)=E\E), (7.15) 

with a superposition, 

\E) =a\l)+b\r). (7.16) 

The qubit Hamiltonian is given by the matrix elements 
of the full Hamiltonian, Eq. (|7.15|) . with respect to the 
states |Z) and |r), 

H u = Ei + {l\U - Ui\l), 
H rr = E r + (r\U -U r \r), 
H rl =Ei{r\l) + {r\U -Ui\l). (7.17) 

In the diagonal matrix elements, the second terms are 
small because the wave functions are exponentially small 
in the region where the deviation of the approximated 
potential from the true one is appreciable. The off di- 
agonal matrix element is exponentially small because of 
small overlap of the ground state wave functions in the 
left and right wells, and also here the main contribution 
comes from the first term. Since the wave functions can 
be chosen real, the truncated Hamiltonian is symmet- 
ric, Hi r = H r [. Then introducing e = E r — Ei, and 
A/2 = H r i, we arrive at the Hamiltonian of the flux 
qubit, 

# = -i(e<7 z + A<7 x ). (7.18) 

The Hamiltonian of the flux qubit is formally equiva- 
lent to that of the charge qubit, Eq. Ij7.5[l . but the phys- 
ical meaning of the terms is rather different. The flux 
qubit Hamiltonian is written in the flux basis, i.e. the 
basis of the states with certain averaged induced flux, <fii 
and <j) r (rather than the charge basis of the charge qubit). 

The energy spectrum of the flux qubit is obviously the 
same as that of the charge qubit, 

£i, 2 = T^Ve 2 + A 2 , (7.19) 
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FIG. 18: Persistent current flux qubit (PCQ) with 3 junctions. 
The side junctions are identical, while the central junction has 
smaller area. 

as shown in Fig. 1151 However, for the flux qubit the 
dashed lines indicate persistent current states in the ab- 
sence of macroscopic tunneling, and the current degener- 
acy point (e = 0) corresponds to a half-integer bias flux. 
The energy levels are controlled by the bias magnetic flux 
(instead of the gate voltage for the charge qubit). 

At the flux degeneracy point, / = (0 e = it), the 
level spacing is determined by the small amplitude of 
tunneling through macroscopic potential barrier, and the 
wave functions correspond to the cat states, which are 
equally weighted superpositions of the flux states. Far 
from the degeneracy point, the qubit states are almost 
pure flux states. 

The possibility to achieve quantum coherence of 
macroscopic current states in an rf-SQUID with a small 
capacitance Josephson junction was first pointed out in 
1984 by LeggetU^i. However, successful experimental 
observation of the effect was achieved only in 2000 by 
Friedman et al£. 



2. 3-junction SQUID - persistent current qubit (PCQ) 

The main drawback of the flux qubit with a single 
Josephson junction (rf-SQUID) described above concerns 
the large inductance of the qubit loop, the energy of 
which must be comparable to the Josephson energy to 
form the required double-well potential profile. This im- 
plies large size of the qubit loop, which makes the qubit 
vulnerable to dephasing by magnetic fluctuations of the 
environment. One way to overcome this difficulty was 
pointed out by Mooij et al£, replacing the large loop in- 
ductance by the Josephson inductance of an additional 
tunnel junction, as shown in Fig. 1181 The design employs 
three tunnel junctions connected in series in a supercon- 
ducting loop. The inductive energy of the loop is chosen 
to be much smaller than the Josephson energy of the 
junctions. The two junctions are supposed to be identi- 
cal while the third junction is supposed to have smaller 
area, and therefore smaller Josephson and larger charging 
energy. The Hamiltonian has the form, 

H = E c [n\ + n\+ nl/ {1/2 + e)]~ 
Ej [cos 0i + cos 02 + (1/2 + e) cos 3 ] . (7.20) 
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FIG. 19: Potential energy landscape J7(0i,02) = U((/>+y ( l } -) 
of the Delft qubit as a function of the two independent phase 
variables 0i and 02, or equivalently, (j>+ (horizontal axis) and 
<j>- (vertical axis) . Black represents the bottom of the poten- 
tial wells, and white the top of the potential barriers. The 
qubit double-well potential is determined by the potential 
landscape centered around the origin in the horizontal direc- 
tion, and typically has the shape shown in Fig. 16. Courtesy 
of C.H. van der Wal. 



To explain the idea, let us consider the potential energy. 
The three phases are not independent and satisfy the 
relation 01 + 02 + 03 = <fte- Let us suppose that the qubit 
is biased at half integer flux quantum, e = n. Then 
introducing new variables, 0± = (0i ± 02)/2, we have 

C/(0+,0_) = -£, 7 [2cos0_ cos0 + - (1/2 + e) cos 20+]. 

(7.21) 

The potential landscape is shown in Fig. ^3 an d the 
qubit potential consists of the double well structure near 
the points (0+,0_) = (0,0). An approximate form of 
the potential energy is given by 

C/(0+ , 0) » Ej (-2e4>\ + ^ V (7-22) 



Each well in this structure corresponds to clock- and 
counterclockwise currents circulating in the loop. The 
amplitude of the structure is given by the parameter 
eEj, and for e <C 1 the tunneling between these wells 
dominates. Thus this qubit is qualitatively similar to 
the single-junction qubit described above, but the quan- 
titative parameters are different and can be significantly 
optimized. 
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D. Potential qubits 

The superconducting qubits that have been discussed 
in previous sections exploit the fundamental quantum un- 
certainty between electric charge and magnetic flux. This 
uncertainty appears already in the dynamics of a single 
Josephson junction (JJ), which is the basis for elementary 
JJ qubits. There are however other possibilities. One of 
them is to delocalize quantum information in a JJ net- 
work by choosing global quantum states of the network 
as a computational basis. Recently, some rather com- 
plicated JJ networks have been discussed, which have 
the unusual property of degenerate ground state, which 
might be employed for efficient qubit protection against 
decoherence 1 ^ 1 ^ 5 -. 

An alternative possibility of further miniaturization 
of superconducting qubits could be to replace the stan- 
dard Josephson junction by a quantum point contact 
(QPC), using the microscopic conducting modes in the 
JJ QPC, the bound Andreev states, as a computational 
basis, allowing control of intrinsic decoherence inside the 
junction 105 ! 126 . 

To explain the physics of this type of qubit, let us 
consider an rf SQUID (see Fig. [SJ with a junction that 
has such a small cross section that the quantization of 
electronic modes in the (transverse) direction perpendic- 
ular to the current flow becomes pronounced. In such a 
junction, quantum point contact (QPC), the Josephson 
current is carried by a number of independent conduct- 
ing electronic modes, each of which can be considered an 
elementary microscopic Josephson junction characterized 
by its own transparency. The number of modes is propor- 
tional to the ratio of the junction cross section and the 
area of the atomic cell (determined by the Fermi wave- 
length) of the junction material. In atomic-sized QPCs 
with only a few conducting modes, the Josephson current 
can be appreciable if the conducting modes are transpar- 
ent (open modes). If the junction is fully transparent 
(reflectivity R = 0) then current is a well defined quan- 
tity. This will correspond to a persistent current with 
certain direction circulating in the qubit loop. On the 
other hand, for a finite reflectivity (R ^ 0) , the electronic 
back scattering will induce hybridization of the persistent 
current states giving rise to strong quantum fluctuation 
of the current. 

Such a quantum regime is distinctly different from the 
macroscopic quantum coherent regime of the flux qubit 
described in Section IVm where the quantum hybridiza- 
tion of the persistent current states is provided by charge 
fluctuations on the junction capacitor. Clearly, charging 
effects will not play any essential role in quantum point 
contacts, and the leading role belongs to the microscopic 
mechanism of electron back scattering. This mechanism 
is only pronounced in quantum point contacts: in classi- 
cal (large area) contacts, such as junctions of macroscopic 
qubits with areas of several square micrometers, the cur- 
rent is carried by a large number (> 10 4 ) of statistically 
independent conducting modes, and fluctuation of the 
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FIG. 20: Energy spectrum of microscopic bound Andreev lev- 
els; the level splitting is determined by the conact reflectivity. 



net Josephson current is negligibly small. 

In QPCs, the Josephson effect is associated with mi- 
croscopic Andreev levels, localized in the junction area, 
which transport Cooper pairs from one junction electrode 
to the otheipiSLiS. As shown in Fig. I2UI the Andreev 
levels, two levels per conducting mode, lie within the su- 
perconducting gap and have the phase-dependent energy 
spectrum, 



E a = ±A 



)/cos 2 (^/2) + i?sin 2 (</)/2), (7.23) 



(here A is the superconducting order parameter in the 
junction electrodes). For very small reflectivity, R <C 1, 
and phase close to tt (half integer flux bias) the An- 
dreev two-level system is well isolated from the contin- 
uum states. The expectation value for the Josephson 
current carried by the level is determined by the Andreev 
level spectrum, 



I a = 



2e dE a 



(7.24) 



and it has different sign for the upper and lower level. 
Since the state of the Andreev two-level system is deter- 
mined by the phase difference and related to the Joseph- 
son current, the state can be manipulated by driving 
magnetic flux through the SQUID loop, and read out 
by measuring circulating persistent currentji22ii2&. 

This microscopic physics underlines recent proposal for 
Andreev level qubiliiS 5 ^ 6 ". The qubit is similar to the 
macroscopic flux qubits with respect to how it is manip- 
ulated and measured, but the great difference is that the 
quantum information is stored in the microscopic quan- 
tum states. This difference is reflected in the more com- 
plex form of the qubit Hamiltonian, which consists of 
the two-level Hamiltonian of the Andreev levels strongly 
coupled to the quantum oscillator describing phase fluc- 
tuations, 



cos - a z 



R sin - <jy 
2 



Hos 



(7.25) 



23 



H osc [4>] = Ecfi+{Ei J /2){4>~(f) e )' 2 . Comparing this equa- 
tion with e.g. the SCT Hamiltonian 1)7. 8|) . we find that 
the truncated Hamiltonian of the SCB is replaced here 
by the Andreev level Hamiltonian. 



VIII. QUBIT READOUT AND MEASUREMENT 
OF QUANTUM INFORMATION 

In this section we present a number of proposed, and 
realized, schemes for measuring quantum states of vari- 
ous superconducting qubits. 

A. Readout: why, when and how? 

As already mentioned in Section IIII Dl the ultimate 
objective of a qubit readout device is to distinguish the 
eigenstates of a qubit in a single measurement "without 
destroying the qubit", a so called "single-shot" quantum 
non-demolition (QND) projective measurement. This 
objective is essential for several reasons: state prepara- 
tion for computation, readout for error correction during 
the calculation, and readout of results at the end of the 
calculation. Strictly speaking, the QND property is only 
needed if the qubit must be left in an eigenstate after 
the readout. In a broader sense, readout of a specific 
qubit must of course not destroy any other qubits in the 
system. 

It must be carefully noted that one cannot "read out 
the state of a qubit" in a single measurement - this is pro- 
hibited by quantum mechanics. It takes repeated mea- 
surements on a large number of replicas of the quantum 
state to characterize the state of the qubit (Eq. (EQJ) - 
" quantum tomography"^i. 

The measurement connects the qubit with the open 
system of the detector, which collapses the combined sys- 
tem of qubit and measurement device to one of its com- 
mon eigenstates. If the coupling between the qubit and 
the detector is weak, the eigenstates are approximately 
those of the qubit. In general however, one must con- 
sider the eigenstates of the total qubit-detector system 
and manipulate gate voltages and fluxes such that the 
readout measurement is performed in a convenient en- 
ergy eigenbasis (see e.gi 42 i 79 ). 

Even under ideal conditions, a single-shot measure- 
ment can only determine the population of an eigenstate 
if the system is prepared in an eigenstate: then the an- 
swer will always be either "0" or "1". If an ideal single- 
shot measurement is used to read out a qubit superpo- 
sition state, e.g. during Rabi oscillation, then again the 
answer can only be "0" or "1". To determine the qubit 
population (i.e. the |ai| 2 and | (3,2 1 2 probabilities) requires 
repetition of the measurement to obtain the expectation 
value. During the intermediate stages of quantum com- 
putation one must therefore not perform a measurement 
on a qubit unless one knows, because of the design and 
timing of the algorithm, that this qubit is in an energy 




FIG. 21: Measurement of the phase qubit. Long-living levels 
form the qubit, while the dashed line indicates a leaky level 
with large energy. 

eigenstate. Then the value is predetermined and the 
qubit left in the eigenstate (Stern-Gerlach-style) . 

On the other hand, to extract the desired final result it 
may be necessary to create an ensemble of calculations to 
be able to perform a complete measurement to determine 
the expectation values of variables of interest, performing 
quantum state tomographyA^i. 



B. Direct qubit measurement 

Direct destructive measurement of the qubit can be il- 
lustrated with the example of a single J J (phase) qubit, 
Section I VI I Al After the manipulation has been per- 
formed (e.g. Rabi oscillation), the qubit is left in a su- 
perposition of the upper and lower energy states. To 
determine the probability of the upper state, one slowly 
increases the bias current until it reaches such a value 
that the upper energy level equals Or gets close to) the 
top of the potential barrier, see Fig. [2] Then the junc- 
tion, being at the upper energy level, will switch from the 
Josephson branch to the dissipative branch, and this can 
be detected by measuring the finite average voltage ap- 
pearing across the junction (voltage state). If the qubit 
is in the lower energy state the qubit will remain on the 
Josephson branch and a finite voltage will not be de- 
tected (zero- voltage state). An alternative method to 
activate switchingafi is to apply an rf signal with reso- 
nant frequency (instead of tilting the junction potential) 
in order to excite the upper energy level and to induce the 
switching event, see Fig. 1211 (also illustrating a standard 
readout method in atomic physics). 

It is obvious that, in this example, the qubit upper 
energy state is always destroyed by the measurement. 
Single-shot measurement is possible provided the MQT 
rate for the lower energy level is sufficiently small to 
prevent the junction switching during the measurement 
time. It is also essential to keep a sufficiently small rate of 
interlevel transitions induced by fluctuations of the bias 
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FIG. 22: Single Electron Transistor (SET) capacitively cou- 
pled to an SCB. 

current and by the current ramping. 

A similar kind of direct destructive measurement was 
performed by Nakamura et ali to detect the state of the 
charge qubit. The qubit operation was performed at the 
charge degeneracy point, u g = 1, where the level splitting 
is minimal. An applied gate voltage then shifted the SCB 
working point (Fig. I15|) . inducing a large level splitting 
of the pure charge states |0) and |1) (the measurement 
preparation stage). In this process the upper |1) charge 
state went above the threshold for Cooper pair decay, 
creating two quasi-particlcs which immediately tunnel 
out via the probe junction into the leads. These quasi- 
particles were measured as a contribution to the classical 
charge current by repeating the experiment many times. 
Obviously, this type of measurement is also destructive. 



C. Measurement of charge qubit with SET 

Non-destructive measurement of the charge qubit has 
been implemented by connecting the qubit capacitively 
to a SET electrometer—. The idea of this method is to 
use a qubit island as an additional SET gate (Fig. I22fl . 
controlling the dc current through the SET depending 
on the state of the qubit. When the measurement is to 
be performed, a driving voltage is applied to the SET, 
and the dc current is measured. Another version of the 
measurement procedure is to apply rf bias to the SET 
(rf-SET 69 i 133 i 134 i 135 ) in Fig. and to measure the dis- 
sipative or inductive response. In both cases the trans- 
missivity will show two distinct values correlated with 
the two states of the qubit. Yet another version has re- 
cently been developed by the NEC groupie to perform 
single-shot readout: the Cooper pair on the SCB island 
then tunnels out onto a trap island (instead of the leads) 
used as a gate to control the current through the SET. 

The physics of the SET-based readout has been exten- 
sively studied theoretically (see^ii 3 iiS and references 
therein). A similar idea of controlling the transmission 
of a quantum point contact (QPC) (instead of an SET) 
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FIG. 23: Probability distributions P of counted electrons as 
functions of time after the turning on the measurement beam 
of electrons. Courtesy of G. Johansson, Chalmers. 



capacitively coupled to a charge qubit has also been ex- 
tensively discussed in literaturo 132 i 139 i 140 i 141 i 142 i 143 i 144 . 

The induced charge on the SET gate depends on the 
state of the qubit, affecting the SET working point 
and determining the conductivity and the average cur- 
rent. The development of the probability distributions 
of counted electrons with time is shown in Fig. E3 

As the number of counted electrons grows, the distribu- 
tions separate and become distinguishable, the distance 
between the peaks developing as ~ N and the width 
~ \/]V. Detailed investigations^ 4 " show that the two 
electron-number probability distributions correlate with 
the probability of finding the qubit in either of two energy 
levels. The long-time development depends on the inten- 
sity and frequency distribution of the back-action noise 
from the electron current. With very weak detector back 
action, the qubit can relax to |0) during the natural re- 
laxation time T\. With very strong back-action noise at 
the qubit frequency, the qubit may become saturated in 
a 50/50 mixed state. 

D. Measurement via coupled oscillator 

Another method of qubit read out that has attracted 
much attention concerns the measurement of the proper- 
ties of a linear or non- linear oscillator coupled to a qubit. 
This method is employed for the measurement of induced 
magnetic flux and persistent current in the loop of flux 
qubits and charge-phase qubits, as well as for charge mea- 
surement on charge qubits. With this method, the qubit 
affects the characteristics of the coupled oscillator, e.g. 
changes the shape of the oscillator potential, after which 
the oscillator can be probed to detect the changes. There 
are two versions of the method: resonant spectroscopy of 
a linear tank circuit/cavity, and threshold detection us- 
ing biased JJ or SQUID magnetometer. 

The first method uses the fact that the resonance fre- 
quency of a linear oscillator weakly coupled to the qubit 
undergoes a shift depending on the qubit state. The 
effect is most easily explained by considering the SCT 
Hamiltonian, Eq. 

H S ct = ~ (ecr z + A(4> e )o- x ) + \(<j) e )4>o- x 

+4E c ii 2 + ^E L 4> 2 . (8.1) 
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Let us proceed to the qubit energy basis, in which casee 
the qubit Hamiltonian takes the form — (E/2)a z , E = 
A 2 . The interaction term in the qubit eigenbasis 



will consist of two parts, the longitudinal part, X z (f>cr z , 
\ z = (A/E)X, and the transverse part, X x <pa x , X x — 
(e/E)X. In the limit of weak coupling the transverse part 
of interaction is the most essential. In the absence of 
interaction (0 e = 0), the energy spectrum of the qubit + 
oscillator system is 



E 



E 



huj(n - 



= \/8EcEl is the plasma frequency of the 
The effect of weak coupling is enhanced in 



where hui 
oscillator. 

the vicinity of the resonance, when the oscillator plasma 
frequency is close to the qubit level spacing, Jko rs E. Let 
us assume, however, that the coupling energy is smaller 
than the deviation from the resonance, X x <C \Tiuj — E\. 
Then the spectrum of the interacting system in the lowest 
perturbative order will acquire a shift, 



6E n± = ± (n + 1) 



Xz.hu 



E l (Huj - E) ' 



(8.3) 



This shift is proportional to the first power of the oscilla- 
tor quantum number n, which implies that the oscillator 
frequency acquires a shift (the frequency of the qubit is 
also shifted 145 i 146 i 147 i 148 i 149 ). Since the sign of the oscil- 
lator frequency shift is different for the different qubit 
states, it is possible to distinguish the state of the qubit 
by probing this frequency shift. 

In the case of the SCT, the LC oscillator is a generic 
part of the circuit. It is equally possible to use an addi- 
tional LC oscillator inductively coupled to a qubit. This 
type of device has been described by Zorinii for SCT 
readout, and recently implemented for flux qubits by 
Il'ichev et ali 27 i 39 . 

Figure [21 illustrates another case, namely a charge 
qubit capacitively coupled to an oscillator, again pro- 
viding energy resolution for discriminating the two qubit 
levels 1 ^. Analysis of this circuit is similar to the one dis- 
cussed below in the context of qubit coupling via oscilla- 
tors, Section ITXl The resulting Hamiltonian is similar to 
Eq. jCTty . namely, 



H = Hscb + Xa y <t> + H 



(8.4) 

In comparison with the case of the SCT, Eq. 18.4fl has 
a different form of the coupling term, which does not 
change during rotation to the qubit eigenbasis. There- 
fore the coupling constant A directly enters Eq. I|8.3|l . 
Recently, this type of read out has been implemented for 
a charge qubit by capacitively coupling the SCB of the 
qubit to a superconducting strip resonato»iSi*i^2iiS. 



E. Threshold detection 
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FIG. 24: SCT qubit coupled to a readout oscillator. The 
qubit is operated by input pulses u(t). The readout oscillator 
is controlled and driven by ac microwave pulses V g (t). The 
output signal will be ac voltage pulses Vout(t), the amplitude 
or phase of which may discriminate between the qubit "0" 
and "1" states. 




FIG. 25: SCT qubit coupled to a JJ readout quantum oscil- 
lator. The JJ oscillator is controlled by dc/ac current pulses 
h(t) adding to the circulating currents in the loop due to the 
SCT qubit. The output will be dc/ac voltage pulses V ut(t) 
discriminating between the qubit "0" and "1" states. 



inserted in the qubit loop, as shown in Fig. l2~5l 

When the measurement of the qubit state is to be per- 
formed, a bias current is sent through the additional junc- 
tion. This current is then added to the qubit-state depen- 
dent persistent current circulating in the qubit loop. If 
the qubit and readout currents flow in the same direction, 
the critical current of the readout JJ is exceeded, which 
induces the junction switching to the resistive branch, 
sending out a voltage pulse. This effect is used to distin- 
guish the qubit states. The method has been extensively 
used experimentally by Vion et ali2ii22iS^. 

To describe the circuit, we add the Lagrangian of a 
biased J J, Eq. (|5.11|) . 



To illustrate the threshold-detection method, let us 
consider an SCT qubit with a third Josephson junction 



^2 12 



4EV 



■ET] 



2e 



(8.5) 



26 



to the SCT Lagrangian (generalized Eq. I|5.30[l . cf. Eqs. 
fiTl . fTTT) i. 



2eC„ 



LscT = je~ c y- " h<^ V ° 

—2Ej cos <j> + cos 



4£ c 



^L0 2 (8.6) 



(here we have neglected a small contribution of the gate 
capacitance to the qubit charging energy). The phase 
quantization condition will now read: 2<j) + + <p = <p e + <j). 
The measurement junction will be assumed in the phase 



regime, EJ 



^ E™, and moreover, the inductive energy 



will be the largest energy in the circuit, El ^> EJ. The 
latter implies that the induced phase is negligibly small 
and can be dropped from the phase quantization condi- 
tion. We also assume that 4> e = 0, thus 2</> + + <j> = 0. 
Then, after having omitted the variable 0+, the kinetic 
energy term of the qubit can be combined with the much 
larger kinetic energy of the measurement junction lead- 
ing to insignificant renormalization of the measurement 
junction capacitance, C m + C/4 — * C m . As a result, the 
total Hamiltonian of the circuit will take the form, 



H = E c (h- - n g ) - 2Ejcos ( — ) cos< 



n 2 -EJ cos, 



2e 



•7) 



Since the measurement junction is supposed to be almost 
classical, its phase is fairly close to the minimum of the 
junction potential. During qubit operation, the bias cur- 
rent is zero; hence the phase of the measurement junction 
is zero. When the measurement is made, the current is 
ramped to a large value close to the critical current of 
the measurement junction, I e — (2e/h)EJ — SI, tilting 
the junction potential and shifting the minimum towards 
7r/2. Introducing a new variable <p = it + 9, we expand 
the potential with respect to small 9 1 and, truncating 
the qubit part, we obtain 



2 2V 2 



+EVi h 2 



where A = 2\/~Ej. The ramping is supposed to be adi- 
abatic so that the phase remains at the minimum point. 
Let us analyze the behavior of the potential minimum 
by omitting a small kinetic term and diagonalizing the 
Hamiltonian (|8.8|l . The corresponding eigenenergies de- 
pend on 9, 



E±(9) 



E 



E r i 
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2e 4E 



(8.9) 



as shown in Fig. 1261 Then within the interval of the bias 
currents, \SI\ < -(2e/h)(A 2 /AE), the potential energy 
corresponding to the ground state has a local minimum, 
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FIG. 26: Josephson potential energy of the measurement 
junction during the measurement: for the "0" qubit eigen- 
state there is a well (full line) confining a level, while for the 
"1" qubit state there is no well (dashed line). 



while for the excited state it does not. This implies that 
when the junction is in the ground state, no voltage will 
be generated. However, if the junction is in the excited 
state, it will switch to the resistive branch, generating a 
voltage pulse that can be detected. 

With the discussed setup the direction of the persis- 
tent current is measured. It is also possible to arrange 
the measurement of the flux by using a dc SQUID as a 
threshold detector. Such a setup is suitable for the mea- 
surement of flux qubits. Let us consider, for example, the 
three junction flux qubit from Section IVII CI inductively 
coupled to a dc SQUID. Then, under certain assump- 
tions, the Hamiltonian of the system can be reduced to 
the following form: 



H 



(ea z + Aa x ) 



+E S C fi2 - (E} + Acr z ) cos 4> SI . 

2e 



(8.10) 



where Ej is an effective (bias flux dependent) Josephson 
energy of the SQUID, Eq. (|5.24|) . and A is an effective 
coupling constant proportional to the mutual inductance 
of the qubit and the SQUID loops. 



IX. PHYSICAL COUPLING SCHEMES FOR 
TWO QUBITS 

A. General principles 

A generic scheme for coupling qubits is based on 
the physical interaction of linear and non-linear os- 
cillators constituting a superconducting circuit. The 
Hamiltonians for the SCB, rf-SQUID, and plain J J con- 
tain quadratic terms representing the kinetic and poten- 
tial energies plus the non-linear Josephson energy term, 
which is quadratic when expanded to lowest order, 



H = E c (n 



Ej(l 



(9.1) 
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FIG. 27: Fixed inductive (flux) coupling of elementary flux 
qubit. The loops can be separate, or have a common leg like 
in the figure. 



H = E c n 2 +Ej(l - cose 



E, 



0- 



H = E C n 2 + Ej(l - cos <p) + —I, 

2e 



(9.2) 



(9.3) 



In a multi-qubit system the induced gate charge in the 
SCB, or the flux through the SQUID loop, or the phase 
in the Josephson energy, will be a sum of contributions 
from several (in principle, all) qubits. The energy of the 
system can therefore not be described as the sum of two 
independent qubits because of the quadratic dependence, 
and the cross terms represent interaction energies of dif- 
ferent kinds: capacitive, inductive and phase/current. 
Moreover, using JJ circuits as non-linear coupling ele- 
ments we have the advantage that the direct physical 
coupling strength may be controlled, e.g tuning the in- 
ductance via current biased JJs, or tuning the capaci- 
tance by a voltage biased SCB. 



B. Inductive coupling of flux qubits 

A common way of coupling flux qubits is the inductive 
coupling: magnetic flux induced by one qubit threads 
the loop of another qubit, changing the effective external 
flux. This effect is taken into account by introducing the 
inductance matrix L^, which connects flux in the i-th 
loop with the current circulating in the A:-th loop, 



4> 



— ^ Likh- 



(9.4) 



The off-diagonal element of this matrix, L12, is the mu- 
tual inductance which is responsible for the interaction. 
By using the inductance matrix, the magnetic part of the 
potential energy in Eq. I|6.16[) can be generalized to the 
case of two coupled qubits, 



) 2 ^(i- 1 ) 4fc ( 



i)(4>k 



■>ck 



(9.5) 



Then following the truncation procedure explained in 
Section lYl! CI we calculate the matrix elements, 



(l\4>-f\l), (r\4>-f\r), (l\4>-f\r), 



(9.6) 
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FIG. 28: Fixed capacitive coupling of charge qubits 



for each qubit. The last matrix element is exponentially 
small, while the first two ones are approximately equal to 
the minimum points of the potential energy, 4>i and <p ri 
respectively. This implies that the truncated interaction 
basically has the zz-form, 

Hint = A(T 2 i0' 2 2, 



C. Capacitive coupling of charge qubits 

One of the simplest coupling schemes is the capacitive 
coupling of charge qubits. Such a coupling is realized by 
connecting the islands of two SCBs via a small capacitor, 
as illustrated in Fig. |2H1 This will introduce an addi- 
tional term in the Lagrangian of the two non-interacting 
SCBs, Eq. H5.27fl . namely the charging energy of the 
capacitor C3, 



5L 



C 3 V 2 



(9.8) 



The voltage drop V3 over the capacitor is expressed via 
the phase differences across the qubit junctions, 



v 3 = 7r{4>i ~ <fe), 
2e 

and thus the kinetic part of the Lagrangian 
take the form 

1 f Ti \ ^ 

K(<j>i,<j>2) = 77 ( — ) c ik $i§ k 



(9.9) 



will 



2 \2e 



i , k 
2 



^TrX/"'" 1 '''' ' 



(9.10) 



where the capacitance matrix elements are Ca = C^i + 
C3 , and C12 = C3 . Then proceeding to the circuit quan- 
tum Hamiltonian as described in Section IVT1 we find the 
interaction term, 



H int = 2e 2 (C 1 )i2nin 2 . 



(9.11) 
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FIG. 30: Capacitive coupling of single JJ qubits 



FIG. 29: Fixed phase coupling of charge qubits 



This interaction term is diagonal in the charge basis, and 
therefore leads to the zz-interaction after truncation, 



Hi, 



Acr z lCT 2 2, A = — (C 1 )i 



(9.12) 



The qubit Hamiltonians are given by Eq. I|7.5|l with 
charging energies renormalized by the coupling capaci- 
tor. 



L 



SCB 



C 



SCB 



D. J J phase coupling of charge qubits 



FIG. 31: Two charge qubits coupled to a common LC- 
oscillator. 



Instead of the capacitor, the charge qubits can be con- 
nected via a Josephson junction 1 ^, as illustrated in Fig. 

121 

In this case, the Josephson energy of the coupling junc- 
tion Ejs cos(0i — <f>2) must be added to the Lagrangian 
in addition to the charging energy. This interaction term 
is apparently off-diagonal in the charge basis and, after 
truncation, gives rise to xx- and yy-couplings, 

Hint — A(<7 a l<7 x 2 + CTylCT,^), A = — ^p, (9.13) 

or equivalently, 

Hint = 2A(CT +! iCr_ !2 + o-,i(J+a). (9.14) 

E. Capacitive coupling of single JJs 

Capacitive coupling of JJ qubits, illustrated in Fig. PUl 
is described in a way similar to the charge qubit, in terms 
of the Lagrangian Eqs. Ij9.8|l . Ij9.9|) . and the resulting 
interaction Hamiltonian has the form given in Eq. I|9.11|l . 

Generally, in the qubit eigenbasis, |0) and |1), all ma- 
trix elements of t he interaction Hamiltonian are non-zero. 
However, if we adopt a parabolic approximation for the 
Josephson potential, then the diagonal matrix elements 
turn to zero, noo = tin = 0, while the off-diagonal matrix 
elements remain finite, not — — nxo = —i(Ej / 'Ec) 1 ^ 4 ■ 
Then, after truncation, the charge number operator h 



turns to (T y , and the qubit-qubit interaction takes the 
yy-iovxa, 



H int =\* yl a y2 , \ = 2e\ h y^ p2 (C' 1 )^. (9.15) 



F. Coupling via oscillators 

Besides the direct coupling schemes described above, 
several schemes of coupling qubits via auxiliary oscilla- 
tors have been considered^. Such schemes provide more 
flexibility, e.g. to control qubit interaction, to couple two 
remote qubits, and to connect several qubits. Moreover, 
in many advanced qubits, the qubit variables are generi- 
cally connected to the outside world via an oscillator (e.g. 
the Delft and Saclay qubits). To explain the principles 
of such a coupling, we consider the coupling scheme for 
charge qubits suggested by Shnirman et ah 1 - 4 . 



1. Coupling of charge (SCB, SCT) qubits 

In this circuit the island of each SCB is connected to 
ground via a common LC-oscillator, as illustrated in Fig. 
|3*T1 The kinetic energy l|5.27|l of a single qubit should 
now be modified taking into account the additional phase 



29 



difference cf> across the oscillator, 

2 



1 



2C(f> 2 _ t + c g (v g i~(f ) -(j)-, i y 



(9.16) 

The cross term in this equation can be made to vanish 
by a change of qubit variable, 



= 0i - a<< 



Cy 



(9.17) 



The kinetic energy will then split into two independent 
parts, the kinetic energy of the qubit in Eq. (|5.28|l , and 
an additional quadratic term, 



CC, 



9 A? 



(9.18) 



which should be combined with the kinetic energy of the 
oscillator, leading to renormalization of oscillator capac- 
itance. 

Expanding the Josephson energy, after the change of 
variable, gives 

Eji cos(4>i — a<j>) « Eji cos & — Ejia6sm6i . (9.19) 

provided the amplitude of the oscillations of <fi is small. 
The last term in this equation describes the linear cou- 
pling of the qubit to the LC-oscillator. 

Collecting all the terms in the Lagrangian and perform- 
ing quantization and truncation procedures, we arrive at 
the following Hamiltonian of the qubits coupled to the 
oscillator (this is similar to Eq. (|8.4|) for the SCT) , 



H — 2_j {HscB,i + Xi<J y i(f>) + H osc , 



(9.20) 



i=l,2 

where H^ CB is given by Eq. I|7.5|l . and 

EjiCg 



A, 



is the coupling strength, and 

Hose = Ecoscfi 2 + E L <t> 2 /2, 



(9.21) 



(9.22) 



is the oscillator Hamiltonian where the term in Eq. (|9.18(l 
has been included. 

The physics of the qubit coupling in this scheme is 
the following: quantum fluctuation of the charge of one 
qubit produces a displacement of the oscillator, which 
perturbs the other qubit. If the plasma frequency of the 
LC oscillator is much larger than the frequencies of all 
qubits, then virtual excitation of the oscillator will pro- 
duce a direct effective qubit-qubit coupling, the oscillator 
staying in the ground state during all qubit operations. 
To provide a small amplitude of the zero-point fluctu- 
ations, the oscillator plasma frequency should be small 
compared to the inductive energy, or Ecosc *C Er, . Then 
the fast fluctuations can be averaged out. Noticing that 



ISI 



JJ 



SCT 



SCT 



FIG. 32: Charge (charge-phase) qubits coupled via a com- 
mon Josephson junction providing phase coupling of the two 
circuits 



the displacement does not change the oscillator ground 
state energy, which then drops out after the averaging, 
we finally arrive at the Hamiltonian of the direct effective 
qubit coupling, 



Hint 



Ai A 



1^2 



Ej 



(9.23) 



for the oscillator-coupled charge qubits in Fig. |3*T1 



2. Current coupling of SCT qubits 

Charge qubits based on SCTs can be coupled by con- 
necting loops of neighboring qubits by a large Josephson 
junction in the common linli 155 i 156 i 15 V 58 i 15 V 60 i 161 , as il- 
lustrated in Fig. 123 

The idea is similar to the previous one: to couple qubit 
variables to a new variable, the phase of the coupling 
Josephson junction, then to arrange the phase regime 
for the junction with large plasma frequency [Eccou-pi <C 
Ejcoupi), and then to average out the additional phase. 
Technically, the circuit is described using the SCT Hamil- 
tonian, Eqs. (|7.7|l . (|7.8|) . for each qubit, 



H S ct = E c {h- - n g ) 2 + E c n 2 + 
(26+ - 6 e ) 2 

- 2Ej cos 6 + cos 6_ + E L v Y+ - reJ , 



(9.24) 



and adding the Hamiltonian of the coupling junction, 



H c = E c ,cK - Ej^ c cos 



(9.25) 



The phase 6 C across the coupling junction must be added 
to the flux quantization condition in each qubit loop; 
e.g., for the first qubit 26 +t i + 6 C = + 6\ (for the 
second qubit the sign of 4> c will be minus). Assuming 
small inductive energy, El <^ Ej tC , we may neglect 6; 
then assuming the flux regime for the coupling Joseph- 
son junction we adopt a parabolic approximation for the 
junction potential, Ej tC & 2 /2. 
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With these approximations, the Hamiltonian of the 
first qubit plus coupling junction will a take form sim- 
ilar to Eq. 19.24J where Ej^ c will substitute for El, and 
<p c will substitute for 20+ — 4> e . Finally assuming the am- 
plitude of the </> c -oscillations to be small, we proceed as 
in the previous subsection, i.e. expand the cosine term 
obtaining linear coupling between the SCB and the oscil- 
lator, truncate the full Hamiltonian, and average out the 
oscillator. This will yield the following interaction term, 

Hi n t = ~~~ (7x1^x2, Aj = -Ejsin^- (9.26) 

Ej, c * 

This coupling scheme also applies to flux qubits: in this 
case, the coupling will have the same form as in Eq. I|9.7|l . 
but the strength will be determined by the Josephson 
energy of the coupling junction, cf. Eq. (|9.26|) . rather 
than by the mutual inductance. 
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FIG. 33: Flux transformer with variable coupling controlled 
by a SQUID. 



on the controlling flux is given by the equatio 
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G. Variable coupling schemes 

Computing with quantum gate networks basically as- 
sumes that one-and two-qubit gates can be turned on 
and off at will. This can be achieved by tuning qubits 
with fixed, finite coupling in and out of resonance, in 
NMR-style computing^. 

Here we shall discuss an alternative way, namely 
to vary the strength of the physical coupling between 
nearest-neighbour qubits, as discussed in a number of re- 
cent papcr a 155 - 156 - 158 - 159 - 160 - 163 - 164 - 165 - 166 



1. Variable inductive coupling 

To achieve variable inductive coupling of flux qubits 
one has to be able to the control the mutual inductance 
of the qubit loops. This can be done by different kinds 
of controllable switches (SQUIDS, transistors )iS in the 
circuit. In a recent experiment, a variable flux trans- 
former was implemented as a coupling element (see Fig. 
13311 by controlling the transforming ratioi££. The flux 
transformer is a superconducting loop strongly induc- 
tively coupled to the qubit loops, which are distant from 
each other so that the direct mutual qubit inductance is 
negligibly small. Because of the effect of quantization of 
magnetic flux in the transformer loopAi^, the local vari- 
ation of magnetic flux <!>i induced by one qubit will af- 
fect a local magnetic flux $2 in the vicinity of the other 
qubit creating effective qubit-qubit coupling. When a dc 
SQUID is inserted in the transformer loop, as shown in 
Fig. 1331 it will shortcircuit the transformer loop, and the 
transformer ratio $2/$! will change. The effect depends 
on the current flowing through the SQUID, and is propor- 
tional to the critical current of the SQUID. The latter is 
controlled by applying a magnetic flux $ ca; to the SQUID 
loop, as explained in Section IV CI and shown in Fig. 13*31 
Quantitatively, the dependence of the transformer ratio 



$2 A , Ej TT$ CX 
T~ = 1 + 7T COS — 

$1 V E L $0 



(9.27) 



where Ej is the Josephson energy of the SQUID junction, 
and El is the inductive energy of the transformer. 



2. Variable Josephson coupling 

A variable Josephson coupling is obtained when a sin- 
gle Josephson junction is substituted by a symmetric dc 
SQUID whose effective Josephson energy 2Ej cos(<j> e /2) 
depends on the magnetic flux threading the SQUID loop 
(see the discussion in Section IV CJI) . This property is 
commonly used to control level spacing in both flux and 
charge qubits introduced in Section IV Dl and it can also 
be used to switch on and off qubit-qubit couplings. For 
example, the coupling of the charge-phase qubits via 
Josephson junction in Fig. 1321 can be made variable 
by substituting the single coupling junction with a dc 
SQUID 155 ! 156 . 

The coupling scheme discussed in Section IIXF II is 
made controllable by using a dc SQUID design for the 
SCB as explained in Section lVD|) . Indeed, since the cou- 
pling strength depends on the Josephson energy of the 
qubit junction, Eq. (|9.21() . this solution provides vari- 
able coupling of the qubits. Similarly, the coupling of 
the SCTs considered in Section liX F 21 can be made con- 
trollable by employing a dc SQUID as a coupling ele- 
ment. A disadvantage of this solution is that the qubit 
parameters will vary simultaneously with varying of the 
coupling strength. A more general drawback of the dc 
SQUID-based controllable coupling is the necessity to 
apply magnetic field locally, which might be difficult to 
achieve without disturbing other elements of the circuit. 
This is however an experimental question, and what are 
practical solutions in the long run remains to be seen. 
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FIG. 34: Coupled charge qubits with current-controlled phase 
coupling: the arrow indicates the direction of the controlling 
bias current. 



3. Variable phase coupling 

An alternative solution for varying the coupling 
is based on the idea of controlling the properties 
of the Josephson junction by applying external dc 
current 1 ^! 1 ^' 1 ^, as illustrated in Fig. OH 

Let us consider the coupling scheme of Section IIX F 21 
the coupling strength here depends on the plasma fre- 
quency of the coupling Josephson junction, which in turn 
depends on the form of the local minimum of the junction 
potential energy. This form can be changed by tilting the 
junction potential by applying external bias current (Fig. 
I34f> . as discussed in Section IV Al The role of the exter- 
nal phase bias, e , will now be played by the minimum 
point (f>o of the tilted potential determined by the applied 
bias current, Ej tC sin (j) = (h/2e)I e . Then the interaction 
term will read, 



Hi, 



£ 2 7 sin 2 (0 o /2) 
E j c cos 4> 



(9.28) 



and local magnetic field biasing is not required. 



4- Variable capacitive coupling 



FIG. 35: 
SCB. 



Variable capacitance tuned by a voltage-controlled 



n 2 , and it also depends on the gate voltages of the qubits 
and the coupling junction. In contrast to the previous 
scheme, here the coupling junction is not supposed to be 
in the phase regime; however, it is still supposed to be 
fast, Ej Eji. Then the energy gap in the spectrum 
of the coupling junction is much bigger than the qubit 
energy, and the junction will stay in the ground state 
during qubit operations. Then after truncation, and av- 
eraging out the coupling junction, the Hamiltonian of the 
circuit will take the form, 



H — 'Y2 Hscb,i 



+ eo (c 2 i + CT z2 ) , 



(9.30) 



where the qubit Hamiltonian is given by Eq. I|7.5|) . and 
the function eo is the ground state energy of the coupling 
junction. The latter can be generally presented as a linear 
combination of terms proportional to <7 z \a z2 and <y z \ + 

£o + cr z2 ) = a + va zl a z2 + /3 (a zl + a z2 ) ■ (9.31) 

with coefficients depending on the gate potentials. The 
second term in this expression gives the zz-coupling (in 
the charge basis), and the coupling constant v may, ac- 
cording to the analysis of RefiiS^, take on both positive 
and negative values depending on the coupling junction 
gate voltage. In particular it may turn to zero, implying 
qubit decoupling. 



Variable capacitive coupling of charge qubits based on 
a quite different physical mechanism of interacting SCB 
charges has been proposed in Refii§£. The SCBs are then 
connected via the circuit presented in Fig. 1351 

The Hamiltonian of this circuit, including the charge 
qubits, has the form 

H = J2 HscB^+Ec{n-q{ni+h 2 )) 2 -E., cos^, (9.29) 

i 

where Ec and Ej ~ E c are the charging and Josephson 
energies of the coupling junction, and h and cj> are the 
charge and the phase of the coupling junction. The func- 
tion q is a linear function of the qubit charges, n±, and 



H. Two qubits coupled via a resonator 

In the previous discussion, the coupling oscillator plays 
a passive role, being enslaved by the qubit dynamics. 
However, if the oscillator is tuned into resonance with a 
qubit, then the oscillator dynamics will become essential, 
leading to qubit-oscillator entanglement. In this case, 
the approximation of direct qubit-qubit coupling is not 
appropriate; instead, manipulations explicitly involving 
the oscillator must be considered. 

Let us consider, as an example, operations with two 
charge qubits capacitively coupled to the oscillator. As- 
suming the qubits to be biased at the degeneracy point 
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and proceeding to the qubit eigenbasis (phase basis in 
this case) , we write the Hamiltonian on the form (cf. Eq. 

<rrm . 

& = - £ \ Y azi ~ Xl(Jx ^ + ^osM- (9.32) 

Let us consider the following manipulation involving the 
variation of the oscillator frequency (cf. Refi§4): at time 
t = 0, the oscillator frequency is off-resonance with both 
qubits, 

MO) < Ai < A 2 - (9.33) 

Then the frequency is rapidly ramped so that the oscil- 
lator becomes resonant with the first qubit, 

M*i) = Ai, (9.34) 

the frequency remaining constant for a while. Then the 
frequency is ramped again and brought into resonance 
with the second qubit, 

Rw(ta) = A 2 . (9.35) 

Finally, after a certain time it is ramped further so that 
the oscillator gets out of resonance with both qubits at 
the end, 

hto{t > t 3 ) > A 2 . (9.36) 

When passing through the resonance, the oscillator is 
hybridized with the corresponding qubit, and after pass- 
ing the resonance, the oscillator and qubit have become 
entangled. For example, let us prepare our system at 
t = in the excited state ip(Q) = 1 100) = |1)|0)|0), 
where the first number denotes the state of the oscilla- 
tor (first excited level), and the last numbers denote the 
(ground) states of the first and second qubits, respec- 
tively. After the first operation, the oscillator will be 
entangled with the first qubit, 

^{t 1 <t<t 2 ) = (cos 6>i 1 10) +sin0 ie iQ |Ol)) |0). 

(9.37) 

After the second manipulation, the state 1 100) will be 
entangled with state |001), 

tl)(t>t 3 ) = cos6»i (cos6» 2 |100) + sin6> 2 e i/3 |001)) 

+ sin(9ie M |010). (9.38) 

To ensure that there are no more resonances during 
the described manipulations, it is sufficient to require 
M0) > A 2 - Ai. 

If the controlling pulses are chosen so that 6* 2 = n/2, 
then the initial excited state will be eliminated form the 
final superposition, and we'll get entangled states of the 
qubits, while the oscillator will return to the ground 
state, 



The manipulation should not necessarily be step-like, it 
is sufficient to pass the resonance rapidly enough to pro- 
vide the Landau-Zener transition, i.e. the speed of the 
frequency ramping should be comparable to the qubit 
level splittings. 

References to recent work on the entanglement of 
qubits and oscillators will be given in Section IV C. 

X. DYNAMICS OF MULTI-QUBIT SYSTEMS 

A. General N-qubit formulation 

A general N-qubit Hamiltonian with general qubit- 
qubit coupling can be written on the form 

\ — / Ai \ 1 V — > 

a=- 2_Jy^o zl + —o xi ) + - OviOvj 

(10.1) 

To solve the Schrodinger equation 

H\tf>) = E\il>) (10.2) 

we expand the N-qubit state in a complete basis, e.g. the 
2 N basis states of the (<j z ) n operator, 

= o |0..00) + ai|0..01) + a 2 |0..10) + ...a (2 ^_ 13 

(10.3) 

and project onto the basis states 

HJ2\P)(PW = EM (10.4) 
p 

obtaining the ususal matrix equation 

(q\H\p)a p = Ea q (10.5) 

p 

where a q = (q\ip), with typical matrix elements given by 

H qp = (q\H\p) = (1..01|tf|0..11) (10.6) 

B. Two qubits, longitudinal (diagonal) coupling 

The first case is an Ising-type model Hamiltonian, rel- 
evant for capacitively or inductively coupled flux qubits, 

H = - ^ (tJ-CTz* + ~Y a xi) + ^12 Vzl&z2 (10-7) 
i=l,2 

We expand a general 2-qubit state in the the a z basis 
q = {«}, 



^(t>t 3 ) = |0) (cos6»ie l/3 |01)) + sin(9ie ia |10)) . 

(9.39) 



|^) = ai|00) + o 2 |01) + a 3 |10) + a 4 \ll) (10.8) 
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With q = {kl} and p = {ran}, we get 



H qp = -{kl\ ^2 ~^ azl + ~^ a xi) + ^12 Ozl(Tz2\mn) 
i=l,2 

ei Ai 
= -(y ( fc Kil m > + — (fc|crxi|m))(/|n) 

-(fc|m)(|(/|a z2 |n)) + ^(/|a K2 |n)) 

+ A12 (fc|cr»i|m)(Z|tr z2 |n) 

(10.9) 

The longitudinal zz-coupling only connects basis states 
with the same indices (diagonal terms), \kl) <-> |fcZ). Eval- 
uation of the matrix elements results in the Hamiltonian 
matrix (setting A12 = A, e\ + e 2 = 2e, t\ — e 2 = 2Ae) 



/ e + X ~IA 2 



H = 



-±A 2 Ae 

W 



\ 



V 



4A a 



-±A X 

-lAi 
-Ae — A -|A 2 



4a 2 



(10.10) 



+ A/ 



The one-body operators can only connect single- 
particle states, and therefore the rest of the state must be 
unchanged (unit overlap). Zero overlap matrix elements 



dealing with pure charge states. The secular equation 
factorises according to 

det (H -E) = 

= [(A 2 - E 2 ) + e(A + E)}[(\ 2 - E 2 ) - e(A + E)] = 

(10.13) 

giving the eigenvalues 

E 1= e + A, E 2)3 = -A, E 4 = -e + X (10.14) 

The corresponding eigenvectors and 2-qubit states are 
given by 



E 1 = e + A 



.E 2 3 = —A 



a 2 = a.3 = 04 = 



= |00) 



ai = 04 = 



(10.15) 
(10.16) 

(10.17) 



|^a) - |01> , (03 = 0) ; = |10) , (02 = 0) (10.18) 



1) = S km , (l\n) = 8i n 



(10.11) 



therefore make some matrix elements zero. 

The two-body operators describing the qubit-qubit 
coupling connect two-particle states with up to two dif- 
ferent indices (single-particle basis states). In the present 
case, the qubit coupling Hamiltonian H2 = A a z \a Z 2 has 
the charge basis as eigenstates and cannot couple differ- 
ent charge states. It therefore only contributes to shifting 
the energy eigenvalues on the diagonal by ±A. In con- 
trast, qubit coupling Hamiltonians involving a x \<J X 2 or 
CyiCy 2 couplings contribute to the off-diagonal matrix 
elements (in the charge basis), and specifically remove 
the zero matrix elements above. 

In the examples below we will specifically consider the 
case when t\ — e 2 = 2Ae = 0, which can be achieved 
by tuning the charging energies via fabrication or gate 
voltage bias. The eigenvalues of the matrix Schrodinger 
equation with the Hamiltonian given by (|10.1Ufl are then 
determined by 

det (H-E) = 0= [(A 2 - E 2 ) + ~(A X + A 2 ) 2 ] 

x [(A 2 - E 2 ) + i(Ai - A 2 ) 2 ] - e 2 (A + E) 2 

(10.12) 



1. Biasing far away from the degeneracy point 

Far away from the degeneracy points in the limit e 3> 
Ai, A 2 , the Hamiltonian matrix is diagonal and we are 



E A = -e + A 



ai = a 2 = 03 = 



N>4 



111) 



(10.19) 



(10.20) 



2. Biasing at the degeneracy point 

Parking both qubits at the degeneracy point, e = 0, 
the secular equation factorises in a different way, 

det (H-E)=Q= [(A 2 - E 2 ) + i(Ai + A 2 ) 2 ] 



x[(A 2 -i? 2 ) + i (A 1 -A 2 ) 2 ] 



again giving a set of exact energy eigenvalues 



£ M = ±\l\ 2 + 7 (Ai + A 



(10.21) 



(10.22) 



^3 = ±yA 2 + -(A!-A 2 ) 2 (10.23) 

where E\ < E 2 < E3 < E4. For E = E\ Al the corre- 
sponding eigenvectors and 2-qubit states are given by 

Ai + A 2 

ai=a 4 , a 2 = a 3 ; , a 2 = -a x — — — — - (10.24) 

A + in, 
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\E X ) = oi(|00) + |11)) + a 2 (|01) + |10)) (10.25) 

\E 4 ) = oa(|00) + |H)) - oi(|01) + |10>) (10.26) 
After normalisation 



"i = ^/l + T^4-r ; a 2 = j^n (10.27) 



\2EA 



\2E X \ 



For E = -©2,3, the corresponding eigenvectors and 2- 
qubit states are given by 

h = -h , h = -b 3 ; , b 2 = -h A ^~^ (10.28) 



\E 2 ) = &i(|00) - |11)) + fiadOl) - |10» (10.29) 
\E 3 ) = 6a(|00> - |H)) - 6i(|01) - |10>) (10.30) 
After normalisation 




3. Two-qubit dynamics under dc-pulse excitation 

We now investigate the dynamics of the two-qubit state 
in the specific case that the dc-pulse is applied equally 
to both gates, e x (t) = e 2 (t), where e x — e x {n gX ) — 
E c (l - 2n gX ), e 2 = e 2 (n g2 ) = E c (l - 2n g2 ), taking the 
system from the the pure charge-state region to the co- 
degeneracy point and back. This is the precise analogy 
of the single-qubit dc-pulse scheme discussed in Section 

wu 

Since the charge state |00) of the starting point does 
not have time to evolve during the steep rise of the dc- 
pulse, it remains frozen and forms the initial state 1 00) 
at time t — at the co-degeneracy point, where it can be 
expanded in the energy eigenbasis, 

|00) = |V(0)> =c 1 \E 1 ) +c 2 \E 2 ) + c 3 \E 3 ) +c 4 |£i) 

(10.32) 

To find the coefficients we project onto the charge basis, 
kl=0,l 

<W|00) = c x (kl\E x ) + c 2 (kl\E 2 ) + c 3 (kl\E 3 ) + c±{kl\E A ) 

(10.33) 

and use the explict results for the energy eigenstates to 
calculate the matrix elements, obtaining 



|00) = oil^i) + h\E 2 ) - b 2 \E 3 ) - a 2 \E 4 ) (10.34) 

For t > this stationary state then develops in time 
governed by the constant Hamiltonian as 



|^(f)) = a x e 



-^lE^+he-^^E? 



Inserting the energy eigenstates from expressions 

frnrnirnrnmrminnni we obtain the time evoiu- 

tion in the charge basis, 

\m) = 



|00) [a\e 


-iE x t _ 


f a \e E ^ 


+ b\e- iE ^ + b\e E - l \ 


|01) [a x a 2 {e 


-iE t t _| 


-e lElt ) + 


hb 2 (e- iE2t + e iE ^)] 


|10) [a x a 2 {e 


-iE t t _| 


_ e ^i*)_ 


hb 2 (e- iE2t + e iE ^)] 


+ |H) [ale 


-iE t t _ 


f a 2 e lBlt 


_ b 2 e ^E 2 t _ 6 2 e iB 2 tj 



(10.36) 

This state could be the input state for another gate oper- 
ation, where some other part of the Hamiltonian is sud- 
denly varied. It might also be that this is the state to 
be measured on. For that purpose one needs in principle 
the probabilities 



Pl (t) = |(00|V(t))| 2 

P2 (i)H(oihA(t))| 2 
P3 (i)H(io|v^))| 2 
p 4 (*) = |(ii|V>(*)>l 2 



(10.37) 
(10.38) 
(10.39) 
(10.40) 



As one example, we will calculate p 3 (t) + pi(t), which 
represents the probability of finding the qubit 1 in the 
upper state |1), independently of the state of qubit 1 
(the states of which are summed over). 

p 3 (t) = | [(a\ + a\) cosE x t ~ (b 2 x + b\) cosE 2 t] 

-i[{a\ - al)smE x t - {bl-bl)smE 2 t) | 2 (10.41) 
p 4 (i) = | [2ia x a 2 sinE x t - 2ib x b 2 sinE 2 t] | 2 (10.42) 

Adding p 3 (t) and pi(t) we obtain 

P3(t)+P4{t) = 

1 1 1 

= — — — cos E x t cos E 2 t + —x sin E x t sin E 2 t 
= i-i(l + x)cos£+t+(l-x)cosi;_i (10.43) 
where E + = E 2 + E x , E- = E 2 — E x , and where x = 

(1 _ 2|£i£ 2 |)- 



C. Two qubits, transverse x-x coupling 

This model is relevant for charge qubits in current- 
coupled loops. The Hamiltonian is now given by 



i=l,2 



-a xl + A a xX <7 x2 (10.44) 



-b 2 e+ iE ^\E 3 ) - a^+^E^) 



(10.35) 



Proceeding as before we have 



H qp = (M\ £ ~^ azi + ^ o- zX a- z2 \mn) 

i=l,2 
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= Cj(kWi 1) \m)-^-(k\a xl \m))(l\n) 

+{k\m)( e j(l\a z2 \n))^(l\a x2 \n)) 

+ A {k\a xl \m){l\a x2 \n) 

(10.45) 

This coupling only connects basis states which differ in 
both indices ("anti-diagonal" terms), i.e. 1 00) <-*■ |11) and 
1 01) 1 10) , meaning that interaction will only appear 
on the anti-diagonal. Evaluation of the matrix elements 
result in the Hamiltonian matrix (again setting A12 = 
A, ei +e 2 = 2e, ei - e 2 = 2Ae) 



H = 



-|Ai 
A 



4A 2 
Aei 
A 

-|Ar 



•iAi A ^ 



A 
-Ae 



■iA 2 -e 



(10.46) 
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FIG. 36: Control pulse sequences involved in quantum state 
manipulations and measurement. Top: microwave voltage 
pulses u(t) are applied to the control gate for state manipula- 
tion. Middle: a readout dc ac pulse (DCP) or ac pulse (ACP) 
Ib(t) is applied to the threshold detector/discriminator a time 
td after the last microwave pulse. Bottom: output signal V(t) 
from the detector. The occurence of a output pulse depends 
on the occupation probabilities of the energy eigenstates. A 
discriminator with threshold Vth converts V(t) into a boolean 
0/1 output for statistical analysis. 



XI. EXPERIMENTS WITH SINGLE QUBITS 
AND READOUT DEVICES 

In this section we shall describe a few experiments 
with single-qubits that represent the current state-of- 
the-art and quite likely will be central components in 
the developmement of multi-qubit systems during the 
next five to ten years. The first experiment presents 
Rabi oscillations induced and observed in the elemen- 
tary phase qubit and readout oscillator formed by a 
single JJ-junction 30 i 31 i 32 i 33 i 34 . The next example de- 
scribes a series of recent experiments with a flux qubili 2 ^ 
coupled to different kinds of SQUID oscillator read- 
out devices^SLiS. A further example will discuss the 
charge-phase qubit coupled to a JJ-junction oscillator 21 
and the recent demonstration of extensive NMR-style op- 
eration of this qubil— . The last example will present the 
case of a charge qubit (single Cooper pair box, SCB) 
coupled to a microwave stripline oscillatoni^Li^SiiSiiiS, 
representing a solid-state analogue of "cavity QED". 

Before describing experiments and results, however, we 
will discuss in some detail the measurement procedures 
that give information about resonance line profiles, Rabi 
oscillations, and relaxation and decohrence times. The 
illustrations will be chosen from Vion et al— and related 
work for the case of the charge-phase qubit, but the ex- 
amples are relevant for all types of qubits, representing 
fundmental procedures for studying quantum systems. 

A. Readout detectors 

Before discussing some of the actual experiments, it 
is convenient to describe some of basic readout-detector 
principles which more or less the same for the SET, rf- 
SET, JJ and SQUID devices. A typical pulse scheme for 
exciting a qubit and reading out the response is shown 
in Fig. IH1 



In Fig. OHIthe readout control pulse can be a dc pulse 
(DCP) or ac pulse (ACP). A DCP readout most often 
leads to an output voltage pulse, which may be quite 
destructive for the quantum system. An ACP readout 
presents a much weaker perturbation by probing the 
ac-response of an oscillator coupled to the qubit, creating 
much less back action, at best representing QND readout. 



1. Spectroscopic detection of Rabi oscillation 

In the simplest use of the classical oscillator, it does 
not discriminate between the two different qubit states, 
but only between energies of radiation emitted by a lossy 
resonator coupled to the qubit. In this way it is possible 
to detect the "low- frequency" Rabi oscillation of a qubit 
driven by continuous (i.e. not pulsed) high-frequency 
radiation tuned in the vicinity of the qubit transition en- 
ergy. If the oscillator is tunable, the resonance window 
can be swept past the Rabi line. Alternatively, the Rabi 
frequency can be tuned and swept past the oscillator win- 
dow by changing the qubit pumping power—. 



2. Charge qubit energy level occupation from counting 
electrons: rf-SET 

In this case, the charge qubit is interacting with a beam 
of electrons passing through a single-electron transistor 
(SET) coupled to a charge qubit (e.g. the rf-SET^), as 
discussed in Section IVIIII and illustrated in Fig. [22 In 
these cases the transmissivity of the electrons will show 
two distinct values correlated with the two states of the 
qubit. 
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FIG. 37: Current-voltage characteristic of tunnel junction 
(solid line) consists of the Josephson branch - vertical line at 
V = 0, and the dissipative branch - curve at eV > 2A. When 
the current is ramped, the junction stays at the Josephson 
branch and when the current approaches the critical value 7 C , 
the junction switches to the dissipative branch (dashed line). 



Coupled qubit-classical-oscillator system: switching 
detectors with dc-pulse (DCP) output 



In Section IVlIII we analyzed the case of an SCT qubit 
current-coupled to a JJ-oscillator (Fig. I25|) and discussed 
the Hamiltonian of the coupled qubit- JJ-oscillator sys- 
tem. The effect of the qubit was to deform the oscillator 
potential in different ways depending on the state of the 
qubit. The effect can then be probed in a number of 
ways, by input and output dc and ac voltage and current 
pulses, to determine the occupation of the qubit energy 
levels. 

Using non-linear oscillators like single JJs or SQUIDS 
one can achieve threshold and switching behaviour where 
the JJ/SQUID switches out of the zero- voltage state, re- 
sulting in an output dc-voltage pulse. 

Switching J J: The method is based on the dependence 
of the critical current of the J J on the state of the qubit, 
and consists of applying a short current DCP to the JJ 
at a value lb during a time At, so that the JJ will switch 
out of its zero- voltage state with a probability P sw (Ib). 
For well-chosen parameters, the detection efficiency can 
approach unity. The switching probability then directly 
measures the qubit 's energy level population. 

Switching SQUID: In the experiments on flux qubits 
by the Delft group, two kinds of physical coupling of 
the SQUID to the qubit have been implemented, namely 
inductive coupling (Fig. [33] (left))Z»iSS and direct cou- 
pling (Fig. Ei (right)) i 24 i 25 i 26 The critical current of 
the SQUID depends on the flux threading the loop, and 
therefore is different for different qubit states. The prob- 
lem is to detect a two percent variation in the SQUID 
critical current associated with a transition between the 
qubit states in a time shorter than the qubit energy relax- 
ation time T\. The SQUID behaves as an oscillator with a 
characteristic plasma frequency lo p = [(L + L^Csh]^ 1 ^ 2 ■ 
This frequency depends on the bias current and on 
the critical current Ic via the Josephson inductance 
Lj = $0 / 277/(7 \/ 1 — lb /Ic (t ne shunt capacitor with ca- 





FIG. 38: Schematics of readout dc SQUID coupled to flux 
qubit; left - inductive coupling, right - direct phase coupling 



pacitance C s h and lead inductance L is used to "tune" 
Ljp). Thus, the plasma frequency takes different values 
^ depending on the state of qubit, 



(0) 



or u> 



p 



representing 

two different shapes of the SQUID oscillator potential. 

In the dc-pulse-triggered switching SQUI D 7 i 24 i 25 , a 
dc-current readout-pulse is applied after the operation 
pulse (s) (Fig. I36fl . setting a switching threshold for the 
critical current. The circulating qubit current for one 
qubit state will then add to the critical current and make 
the SQUID switch to the voltage state, while the other 
qubit state will reduce the current and leave the SQUID 
in the zero-voltage state. 

In an application of ac-pulse-triggered switching 
SQUID 2 -, readout relies on resonant activation by a mi- 
crowave pulse at a frequency close to u> p , adjusting the 
power so that the SQUID switches to the finite voltage 
state by resonant activation if the qubit is in state |0), 
whereas it stays in the zero- voltage state if it is in state 
1 1) . The resonant activation scheme is similar to the read- 
out scheme used by Martinis et al^Sj^i^SiS. (see Section 

ixm . 



Coupled qubit-classical-oscillator system: ac-pulse 
(ACP) non-switching detectors 



This implementation of ACP readout uses the qubit- 
SQUID combination 7 shown in Fig. 1381 (left), but with 
ACP instead of DCP readout, implementing a nonde- 
structive dispersive method for the readout of the flux 
qubrfci§£. The detection is based on the measurement 
of the Josephson inductance of a dc-SQUID inductively 
coupled to the qubit. Using this method, Lupascu et 
alJ^ measured the spectrum of the qubit resonance line 
and obtained relaxation times around 80 fj,s, much longer 
than observed with DCP. 

A related readout scheme was recently implemented by 
Siddiqi et al. 1 — using two different oscillation states of 
the non-linear JJ in the zero-voltage state. 



37 




50 



V (MHz) 



FIG. 39: Left: qubit energy level scheme. The vertical dashed 
line marks the qubit working point and transition energy. The 
arrow marks the detuned microwave excitation. Right: pop- 
ulation of the upper level as a function of the detuning; the 
inverse of the half width (FWHM) of the resonance line gives 
the total decoherence time Ti. Courtesy of D. Esteve, CEA- 
Saclay. 



B. Operation and measurement procedures 

A number of operation and readout pulses can be ap- 
plied to a qubit circuit in order to measure various prop- 
erties. The number of applied microwave pulses can vary 
depending on what quantities are to be measured: reso- 
nance line profile, relaxation time, Rabi oscillation, Ram- 
sey interference or Spin Echo, as discussed below. 



1. Resonance line profiles and T2 decoherence times 

To study the resonance line profile, one applies a sin- 
gle long weak microwave pulse with given frequency, fol- 
lowed by a readout pulse. The procedure is then repeated 
for a spectrum of frequencies. The Rabi oscillation am- 
plitude, the upper state population, and the detector 
switching probability p(t) will depend on the detuning 
and will grow towards resonance. The linewidth gives 
directly the total inverse decoherence lifetime I/T2 = 
l/2Ti + The decoherence-time contributions from 

relaxation (1/Ti) and dephasing (l/T^) can be (approx- 
imately) separately measured, as discussed below. 



2. Ti relaxation times 

To determine the T\ relaxation time one measures the 
decacy of the population of the upper |1) state after a 
long microwave pulse saturating the transition, varying 
the delay time td of the detector readout pulse. The 
measured T\ =1.8 microseconds is so far the best value 
for the Quantronium charge-phase qubit. 
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FIG. 40: Decay of the switching probability of the charge- 
qubit readout junction as a function of the delay time td 
between the excitation and readout pulses. Courtesy of D. 
Esteve, CEA-Saclay. 
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FIG. 41: Left: Rabi oscillations of the switching probability 
measured just after a resonant microwave pulse of duration. 
Right: Measured Rabi frequency (dots) varies linearly with 
microwave amplitude (voltage) as expected. Courtesy of D. 
Esteve, CEA-Saclay. 



3. Rabi oscillations and T^RaU decoherence time 

To study Rabi oscillations (frequency Q ~ u, the am- 
plitude of driving field) one turns on a resonant mi- 
crowave pulse for a given time t^ w and measures the up- 
per |1) state population (probability) pi(i) after a given 
(short) delay time td- If the systems is perfectly coher- 
ent, the state vector will develop as cosftt |0)+sinf2t |1), 
and the population of the upper state will then oscil- 
late as sin 2 Sit between and 1 . In the presence of de- 
coherence, the amplitude of the oscillation of pi (t) will 
decay on a time scale towards the average value 

Pi(t = 00) = 0.5. This corresponds to incoherent satura- 
tion of the to 1 transition. 
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FIG. 42: Ramsey fringes of the switching probability af- 
ter two phase-coherent microwave 7r/2 pulses separated by 
the time delay t. The continuous line represents a fit by 
exponentially damped cosine function with time constant 
I2 = Trj, = 0.5/is. The oscillation period coincides with the 
inverse of the detuning frequency <5 (here S — v — vq\ = 20.6 
MHz). Courtesy of D. Esteve, CEA-Saclay. 
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FIG. 43: Spin-echo experiment. The left part shows the 
basic Ramsey oscillation. The right part shows the echo signal 
appearing in the time window around twice the time delay 
between the first 7r/2-pulse and the 7r-pulse. Courtesy of D. 
Esteve, CEA-Saclay. 



4- Ramsey interference, dephasing and T2, Ramsey 
decoherence time 

The Ramsey interference experiment measures the de- 
coherence time of the non-driven, freely precessing, qubit. 
In this experiment a n/2 microwave pulse around the x- 
axis induces Rabi oscillation that tips the spin from the 
north pole down to the equator. The spin vector rotates 
in the x-y plane, and after a given time At, another tt/2 
microwave pulse is applied, immediately followed by a 
readout pulse. 

Since the 7r/2 pulses are detuned by S from the qubit 
|0) — * |1) transition frequency, the qubit will precess 
with frequency 6 relative to the rotating frame of the 
driving field. Since the second microwave pulse will be 
applied in the plane of the rotating frame, it will have 
a projection cos<5t on the qubit vector and will drive 
the qubit towards the north or south poles, resulting 
in a specific time-independent final superposition state 
cos St |0) + sin 5t |1) of the qubit at the end of the last 
7r/2 pulse. The readout pulse then catches the qubit 
in this superposition state and forces it to decay if the 
qubit is in the upper |1) state. The probability will 
oscillate with the detuning frequency, and a single-shot 
experiment will then detect the upper state with this 
probability. Repeating the experiment many times for 
different pi/2 pulse separation At will then give |0) or 
|1) with probabilities cos 2 St and sin 2 St. Taking the 
average, and then varying the pulse separation, will 
trace out the Ramsey interference oscillatory signal. 
Dephasing will make the signa decay on the timescale T$. 



5. Spin-echo 

The spin-echo and Ramsey pulse sequences differ in 
that a 7r-pulse around the x-axis is added in between the 
two 7r/2-pulses in the spin-echo experiment, as shown in 
FigEHl As in the Ramsey experiment, the first 7r/2-pulse 
makes the Bloch vector start rotating in the equatorial 
x-y plane with frequency E/h = Vq\. The effect of the 
7r-pulsc is now to flip the entire x-y plane with the rotat- 
ing Bloch vector around the x-axis, reflecting the Bloch 
vector in the x-z plane. The Bloch vector then continues 
to rotate in the x-y plane in the same direction. Finally 
a second 7r/2-pulse is applied to project the state on the 
z-axis. 

If two Bloch vectors with slightly different frequency 
start rotating at the same time in the x-y plane, they will 
move with different angular speeds. The effect of the 7r- 
pulse at time At will be to permute the Bloch vectors, 
and then let the motion continue in the same direction. 
This is similar to reversing the motion and letting the 
Bloch vectors back-trace. The net result is that the two 
Bloch vectors re-align after time 2 At. 

In NMR experiments, the different Bloch vectors cor- 
respond to different spins in the ensemble. In the case 
of a single qubit, the implication is that in a series of re- 
peated experiments, the result will be insensitive to small 
variations SE of the qubit energy between measurements, 
as long as the energy (rotation frequency) is constant 
during one and the same measurement. If fluctuations 
occur during one measurement, then this cannot be cor- 
rected for. The spin-echo procedure can therefore remove 
the measurement-related line-broadening associated with 
slow fluctuations of the qubit precession, and allow ob- 
servation of the intrinsic coherence time of the qubit. 
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C. NIST Current-biased Josephson Junction Qubit 



Several experimental groups have realized the Joseph- 
son Junction (JJ) qubit 29 i 30 i 31 i 32 i 33 i 34 . Here we describe 
the experiment performed at NIST 3 ^— In this experi- 
ment the junction parameters and bias current were cho- 
sen such that a small number of well defined levels were 
formed in the potential well (Fig. I21|) . with the inter- 
level frequencies, lo 01 /2it = 6.9GHz and Wi2/27r6.28GHz, 
the quality factor was Q=380. The experiment was per- 
formed at very low temperature, T= 25 mK. 

The qubit was driven from the ground state, |0), to 
the upper state, |1), by the resonance rf pulse with fre- 
quency woi, and then the occupation of the upper qubit 
level was measured. The measurement was performed by 
exciting qubit further from the upper level to the auxil- 
iary level with higher escape rate by applying the second 
resonance pulse with frequency W12. During the whole 
operation, across the junction only oscillating voltage de- 
velops with zero average value over the period. When the 
tunneling event occurred, the junction switched to the 
dissipative regime, and finite dc voltage appeared across 
the junction, which was detected. Alternatively, post- 
measurement classical states "0" and "1" differ in flux 
by $0, which is readily measured by a readout SQUID. 



The relaxation rate was measured in the standard way 
by applying a Rabi pumping pulse followed by a mea- 
suring pulse with a certain delay. Non-exponential relax- 
ation was observed, first rapid, ~ 20 nsec, and then more 
slow, ~ 300 nsec. By reducing the length of the pumping 
pulse down to 25 nsec, i.e. below the relaxation time, 
Rabi oscillations were observed. In this experiment the 
amplitude of the pumping pulse rather than duration was 
varied, which affected the Rabi frequency and allowed the 
observation of oscillation of the level population for fixed 
duration of the pumping pulse. 



The Grenoble groups* has observed coherent oscil- 
lations in a multi-level quantum system, formed by a 
current-biased dc SQUID. These oscillations have been 
induced by applying resonant microwave pulses of flux. 
Quantum measurement is performed by a nanosecond 
flux pulse that projects the final state onto one of two 
different voltage states of the dc SQUID, which can be 
read out. The number of quantum states involved in the 
coherent oscillations increases with increasing microwave 
power. The dependence of the oscillation frequency on 
microwave power deviates strongly from the linear regime 
expected for a two-level system and can be very well ex- 
plained by a theoretical model taking into account the 
anharmonicity of the multi-level system. 




FIG. 44: Upper panel: Scanning electron micrograph of a 
small-loop flux-qubit with three Josephson junctions of crit- 
ical current ~ 0.5 mA, and an attached large-loop SQUID 
with two big Josephson junctions of critical current ~ 2.2 mA. 
Arrows indicate the two directions of the persistent current 
in the qubit. Lower panels: Schematic of the on-chip cir- 
cuit; crosses represent the Josephson junctions. The SQUID 
is shunted by two capacitors to reduce the SQUID plasma fre- 
quency and biased through a small resistor to avoid parasitic 
resonances in the leads. Symmetry of the circuit is introduced 
to suppress excitation of the SQUID from the qubit-control 
pulses. The MW line provides microwave current bursts in- 
ducing oscillating magnetic fields in the qubit loop. The cur- 
rent line provides the measuring pulse It, and the voltage line 
allows the readout of the switching pulse Vout. Adapted from 
Chiorescu et alw^ 



D. Flux qubits 

1. Delft 3-junction persistent current qubit with de-pulse 
(DCP) readout 

The original design of the 3-junction qubit (Fig. ^was 
published by Mooij et al. in 1999&, and the first spec- 
troscopic measurements by van der Wal et al. in 2000 7 
(and simultaneously by Friedman et al£ for a single- 
junction rf-SQUID). Recently the Delft group has also 
investigated designs where the 3-junction qubit is shar- 
ing a loop with the measurement SQUID to increase the 
coupling strength—, as shown in Fig. 1441 

To observe and study Rabi oscillations, the qubit was 
biased at the degeneracy point and the qubit |0) — > |1) 
transition excited by a pulse of 5.71 GHz microwave ra- 
diation of variable length, followed by a bias-current (/&) 
readout pulse applied to the SQUID (Fig. lower left 
panel) . The first (high) part of the readout pulse (about 
10 ns) has two functions: It displaces the qubit away from 
the degeneracy point so that the qubit eigenstates carry 
finite current, and it tilts the SQUID potential so that 
the SQUID can escape to the voltage state if the qubit is 
in its upper state. The purpose of the long lower plateau 
of the the readout pulse is to prevent the SQUID from 
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FIG. 45: AFM picture of the charge-phase qubit ("quantro- 
nium") corresponding to the circuit scheme in Fig. 1251 The 
working point is controlled by two external "knobs" , a voltage 
gate controlling the induced charge (n g ) on the SCT island, 
and a current gate controlling the total phase across the SCT 
via the external flux (<f> e ) threading the loop. The large read- 
out JJ is seen in the left part of the figure. Courtesy of D. 
Esteve, CEA-Saclay. 



returning (" retrapping" ) to the zero- voltage state. With 
these operation and readout techniques, Rabi (driven) 
oscillations, Ramsey (free) oscillations and spin-echos of 
the qubit were observed^*, giving a Ramsey free oscil- 
lation dephasing time of 20 ns and spin echo dephasing 
time of 30 ns. 

As mentioned above, relaxation times around 80 /is 
have been measured with ACP readouliiS, demonstrat- 
ing that the properties of readout devices are criti- 
cally important for observing intrinsic qubit decoherence 
times. Recent further improvements have resulted in 
Ramsey decoherence times Ti^amsey of 200 ns, Rabi 
(driven) decoherence times T^ ^abi of 5 us, and relax- 
ation times Ti of more than 100 /is. 



E. Charge-phase qubit 

1. General considerations 

As described in Section VII, the charge-phase qubit 
circuit consists of a single- Cooper-pair transistor (SCT) 
in a superconducting loop, The Hamiltonian for the SCT 
part of this circuit is given by 

H S CT = -^[e(n g )<T z +A((t> e )<T x ] (11.1) 

where the charging and tunneling parameters are them- 
selves functions of external control parameters, gate 
charge n g and loop flux <f> e , e(n g ) — Ec(l — 2n g ) and 




FIG. 46: Charge- phase qubit energy surface. (N g = n g ; 
8/2n — 4> e .) Courtesy of D. Esteve, CEA-Saclay. 

A((j) e ) — 2Ej cos(0 e /2). The qubit energy levels 

Si^T^K^ + A^ (11.2) 

then form a 2-dimensional landscape as functions of the 
gate charge n g (V) and phase 4> e (<&), which are functions 
of the gate voltage V — C g 2en g and gate magnetic flux 
<3> = (h/2e)(f> e . The energy level surfaces are therefore 
functions of two parameters, gate voltage and flux, giv- 
ing us two independent knobs for controlling the working 
point of the (charge-phase) qubit. Expanding the energy 
in Taylor series around some bias working point (14, $b), 
one obtains 

SE(V, $) = E(V, $) - E(Vb, $&) 

= f^+S'* + 5?V + ^V (11-3) 
SV <5$ 2 SV 2 <5$ 2 

The derivatives are response functions for charge, cur- 
rent, capacitance and inductance (omitting the cross 
term) (i — 1,2), 

8Ei(V, *) = Qt 5V + Ii 6$ + Q SV 2 + U 6<S> 2 (11.4) 

or, equivalcntly, 

5Ei(n g , <f> e ) = Qi 5n g + Ii 6(f> e + Ci 8n g 2 + L, ; 5(j) e 2 

(11.5) 

On the energy level surfaces (Fig. 14611 . the special point 
(n g ,4> e ) — (0.5,0) is an extreme point with zero first 
derivative. This means that the energies of the |0) and 
|1) states will be invariant to first order to small varia- 
tions of charge and phase, which will minimize the qubit 
sensitivity to fluctuations of the working point caused by 
noise. 
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FIG. 47: The charge-phase qubit frequency surface. (N g 
n g ; 8/2-k — 4> e .) Courtesy of D. Esteve, CEA-Saclay. 



8 100 



o 
c 



o 
u 



• Linewidth 

A Ramsey 

■ Detuning pulse 

O Echo 



N =1/2 I s = 




_i_ 



-0.3 



-0.2 
- I 5/2* I 



-0.1 



0.0 0.1 
I Ng- 1/2 I 



FIG. 48: Free-evolution decoherence times for the quantron- 
ium charge-phase SCT qubitiiSS (N g = n 9 ; 5/27T = e .) Full 
and dashed curves represent results of theoretical modeling. 
Courtesy of D. Esteve, CEA-Saclay. 



The point (n g ,4> e ) = (0.5,0) is often referred to as 
the "degeneracy" point because the charging energy is 
zero, e(n g ) = Eq{1 — 2n ff ) = 0. The level splitting at 
this point is determined by the Josephson tunneling in- 
teraction Ej, and is a function of the external bias <f) e ) 
A(<£ e ) = 2Ej cos(0 e /2). 

Figure 07| shows the frequency surface Uqi = AE/h, 
AE = E% — E\. During operation, the qubit is preferably 
parked at the degeneracy point (n g ,<f> e ) = (0.5,0) where 
Q = I = 0, in order to minimize the decohering influence 
of noise. In order to induce a qubit response, for gate 
operation or readout, one must therefore either (i) move 
the bias point away from point (0.5, 0) to have finite first- 
order response with Q ^ or I ^ 0, or (ii) stay at (0.5, 0) 
and apply a perturbing ac-field that makes the second- 
order response significant. 



2. The CEA-Saclay " quantronium" charge-phase qubit 

The "quantronium" charge-phase circuit is given by 
Fig. I45I adding a large readout Josephson junction in 
the loop (cf. Fig. I25fr . The minimum linewidth (Fig. 
corresponds to a Q-value of 20000 and a decoherence time 
T-2 = 0.5 /zs. With the measured relaxation time Ti=1.8 
/zs (Fig. I40|) . the dephasing time can be estimated to 
T^0.8 /zs. 

A set of results for Rabi oscillations of the Saclay 
quantronium qubit is shown in Fig. I41I The decoher- 
ence time represents decoherence under driving condi- 
tions. Another, more fundamental, measure of decoher- 
ence is the free precession dephasing time T 2 when the 
qubit is left to itself. This is measured in in the Ramsey 
two-pulse experiment, as already described above in Fig. 
l4*2l showing the CEA-Saclay data 21 i 22 , giving T 2 = 0.5 /zs. 



The spin-echo results in Ref^ gave a lifetime as long as 
T 2 = lzzs. 

In a recent systematic experimental and theoretical in- 
vestigation of a specific charge-phase device, investigat- 
ing the effects of relaxation and dephasing on Rabi os- 
cillation, Ramsey fringes and spin-echos^^, one obtains 
the following picture of different coherence times for the 
quantronium charge-phase qubit: 

Fluctuations of charge 8n g and flux 5<j) e will shift 
energy levels and make the qubit transition energy 
yj e(n g + Sn g ) 2 + A(<j) e + 54> e ) 2 fluctuate. However, the 
charge degeneracy point is a saddle point, which means 
that at that working point the qubit transition frequency 
is insensitive to low-frequency noise to first order, allow- 
ing long coherence times. As can be seen in Fig. 021 
the coherence is sharply peaked ound the "magic" de- 
generacy point. The spin-echo experiment indicates the 
presence of slow charge fluctuations from perturbing im- 
purity two- level systems (TLS) (1/f noise) and that long 
coherence time can be recovered by spin-echo techniques 
until the decoherence becomes too rapid at significant 
distances from the magic point along the charge axis. In 
contrast, in the experiments with this device the decoher- 
ence due to phase fluctuations could not be significantly 
compensated for by spin-echo techniques, indicating that 
the phase noise (current and flux noise) in this device is 
high-frequency noise. 
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XII. EXPERIMENTS WITH QUBITS COUPLED 
TO QUANTUM OSCILLATORS 

A. General discussion 

The present development of quantum information pro- 
cessing with Josephson Junctions (JJ-QIP) goes in the 
direction of coupling qubits with quantum oscillators, for 
operation, readout and memory. In Section IVIIII we dis- 
cussed the SCT, i.e. a single Cooper-pair transistor in 
a superconducting loopii, providing one typical form of 
the Hamiltonian for a qubit-oscillator coupled system. 
We also showed with perturbation theory how the Hamil- 
tonian gave rise to qubit-dependent deformed oscillator 
potential and shifted oscillator frequency. In I V 1111 in 
addition we discussed the SCT charge-coupled to an LC- 
oscillator, which is also representative for a flux qubit 
coupled to a SQUID oscillator, and which describes a 
charge qubit in a microwave cavity. 

To connect to the language of quantum optics and cav- 
ity QED and the current work on solid-state applications, 
we now explicitly introduce quantization of the oscillator. 
Quantizing the oscillator, ~ (a + + a), a representative 
form of the qubit-oscillator Hamiltonian reads: 



H — H q + Hi n t + H osc (12-1) 

H q = - l -Eo z (12.2) 

Hint = g <?x (a + + a) (12-3) 

H osc = huj (a + a + h (12.4) 



Introducing the step operators a± = a x ± i cr y , the inter- 
action term can be written as 

Hint = 9 (o+ a + c- a + ) + g (a + a + + <r_ a) (12.5) 

The first term 

Hint — g (o+ a + er_ a + ) (12.6) 

represents the resonant (co-rotating) part of the interac- 
tion, while the second term represents the non-resonant 
counter- rotating part. In the rotating- wave approxima- 
tion (RWA) one only keeps the first term, which gives the 
Jaynes-Cummings modeii£iiiZ2iiZ2iiZi. Diagonalizing the 
Jaynes-Cummings Hamiltonian to second order by a uni- 
tary transformation gives 

H = -\(E + ^-)a z + (Twj+^-a z )a+a (12.7) 
2 o o 

where 5 = E — tko is the so called detuning. The re- 
sult illustrates what we have already discussed in detail, 
namely that (i) the qubit transition energy E is shifted 
(renormalized) by the coupling to the oscillator, and (ii) 
the oscillator energy hui is shifted by the qubit in differ- 
ent directions depending on the state of the qubit, which 
allows discriminating the two qubit states. 



I11> 

. I01>-I10> 




FIG. 49: Qubit-oscillator level structure. The notation for 
the states is: \qubit; oscillator) = 1 0/1 ; n — 0, 1, ..); (a) E w 
2huj: very large "detuning", weak coupling, (b) E ~ hw: 
resonance, strong coupling and hybridization, level ("vacuum 
Rabi") splitting. 



Figure 021 shows the basic level structure of the qubit- 
oscillator system in the cases of (a) weak and (b) strong 
coupling. 

What is weak or strong coupling is determined by the 
strength of the level hybridization, which in the end de- 
pends on qubit-oscillator detuning and level degeneracies. 
The perturbative Hamiltonian in Eq. (|12.7|) is valid for 
large detuning (g <C 5 qr ) and allows us to approach the 
case of strong coupling between the qubit and the oscil- 
lator. Close to resonance the degenerate states have to 
be treated non-perturbatively. 

Figure EJ a ) corresponds to the weak-coupling (non- 
resonant) case where the levels of the two subsystems are 
only weakly perturbed by the coupling, adding red and 
blue sideband transitions 1 01) — ► 1 10) and 1 00) — > |11) to 
the main zero-photon transition 1 00) — * |10). 

Figure I49f b) corresponds to the resonant strong- 
coupling case when the 1 01 > and 1 10 > states are de- 
generate, in which case the qubit-oscillator coupling pro- 
duces two "bonding-antibonding" states in the usual way, 
and the main line becomes split into two lines ("vacuum 
Rabi splitting" . Of major importance are the linewidths 
of the qubit and the oscillator relative to the splittings 
caused by the interaction. To discriminate between the 
two qubit states, the oscillator shift must be larger than 
the average level width. 

Since the qubit and the oscillator form a coherent 
multi- level system, as described before fSection IX B 3|l . 
the time evolution can be written as ci (i) |00) +C2 (£) |01) + 
C3(t)|10)+C4(i)|ll), which in general does not reduce to a 
product of qubit and oscillator states, and therefore rep- 
resents (time-dependent) entanglement. This provides 
opportunities for implementing quantum gate operation 
involving qubits and oscillators. 

Generation and control of entangled states can be 
achieved by using microwave pulses to induce Rabi oscil- 
lation between specific transitions of the coupled qubit- 
oscillator system, and the result can be studied by spec- 
troscopy on suitable transitions or by time-resolved de- 
tection of suitable Rabi oscillations. 
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FIG. 50: Resonant frequencies indicated by peaks in the 
SQUID switching probability when a long (300 ns) microwave 
pulse excites (saturates) the system before the readout pulse. 
Data are represented as a function of the external flux & e xt 
through the qubit area away from the qubit symmetry point. 
The blue |00> -» |11> and red |01) -> 1 10) sidebands are shown 
by down- and up-triangles, respectively; continuous lines are 
obtained by adding 2.96 GHz and -2.90 GHz, respectively, 
to the central continuous line (numerical fit). These values 
are close to the oscillator resonance frequency v v at 2.91 GHz 
(solid circles). Courtesy of J.E. Mooij, TU Delft. 



B. Delft persistent current flux qubit coupled to a 
quantum oscillator 

The Delft experiment demonstrates entanglement 
between a superconducting flux qubit and a SQUID 
quantum oscillator. The SQUID provides the measure- 
ment system for detecting the quantum states (threshold 
switching detector, Fig. |4~4"|) . It is also presents an effec- 
tive inductance that, in parallel with an external shunt 
capacitance, acts as a low-frequency harmonic oscillator; 
the qubit and oscillator frequencies are approximately 
hv q = hi>Qi ps 6 GHz and hv r w 3 GHz, corresponding 
to Fig. |l|a). 

In the Delft experiment^, performing spectroscopy on 
the coupled qubit-oscillator multi-level system reveals the 
variation of the main and sideband transitions with flux 
bias Q ex t, as shown in Fig. 1201 The presence of visible 
sideband transitions demonstrates the qubit-oscillator in- 
teraction and (weak) level hybridization). 

Short microwave pulses can now be used to induce Rabi 
oscillations between the various qubit-oscillator transi- 
tions. In the Delft experiment*^, microwave pulses with 
frequency v q — uqi w A/h « 5.9 GHz (qubit symme- 
try point) were first used to induce Rabi oscillations 
with ~ 25 ns decay time at the main qubit transition 
1 00) -> 1 10) (and |01) -> |11)) for different values of 
the microwave power. The Rabi frequency as function 
of the amplitude of the microwave voltage demonstrated 
qubit-oscillator hybridization (avoided crossings) at the 
oscillator v p (y r ) and Larmor A/h (!) frequencies. 

The dynamics of the coupled qubit-oscillator system 
can be studied by inducing microwave-driven Rabi oscil- 
lation between the blue |00) -> ] 11) and red |01) -> | T0> 
sidebands. In particular one can study the conditional 
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FIG. 51: Generation and control of entangled states via n and 
27T Rabi pulses on the qubit transition |00) — * |10), followed 
by Rabi driving of blue sideband transitions 1 00) — >. A n pulse 
excites the system from 1 00) to 1 10) , which suppresses the blue 
sideband transitions 1 00) — > |11) (second curve from the top). 
On the other hand, with a 2n pulse the system returns to 1 00) , 
which allows Rabi oscillation on the blue sideband. Adapted 
from Chiorescu et aim— 



dynamics in microwave multi-pulse experiments: by co- 
herently (de)populating selected levels via proper timing 
of Rabi oscillations induced by a first pulse, a second 
pulse can induce Rabi oscillations on another transition 
connected to the levels controlled by the first pulse. Al- 
ternatively, such Rabi oscillations can instead be blocked 
by the first excitation. 

This is shown experimentally in Fig. ??: by preparing 
the initial state with initial n and 27r pulses, the side- 
band Rabi oscillations could be turned off and on again. 
The corresponding Rabi oscillations are shown in the left 
part of Fig. ??(b), demonstrating rapid decay due to the 
strong damping of the SQUID oscillator (lifetime ~ 3 ns; 
Q=100-150). 

In the previous example, the control pulse (first pulse) 
was applied to the main 1 00) — > 1 10) transition, control- 
ling the populations of the "0" and "1" levels. In the 
next example, microwave control pulse is applied to the 
1 00) — > 1 01) transition, inducing Rabi oscillations which 
populate the first excited state of oscillator. A second 
microwave pulse (in principle, with different frequency) 
can then induce Rabi oscillations on the 1 01) — * | TO) 
transition (red-sideband). The experimental result is 
shown Fig. Clearly, for sufficiently long qubit and os- 
cillator lifetimes, one can prepare entangled Bell states, 
i(|00) ± |01)) and |(|01) ± 1 10)) by applying tt/2 mi- 
crowave pulses to the 1 00) — * 1 1 1 ) and |01) — > | TO) tran- 
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MW pulse length <ns) 

FIG. 52: Generation and control of entangled states via n and 
2ir Rabi pulses on the qubit transition 1 00) — > |10), followed 
by Rabi driving of red sideband transitions 1 1 0) — > |01). A n 
pulse excites the system from 1 00) to 1 10) , which In the right 
panel, a n pulse excites the system from 1 00) to 1 10) , which 
populates the 1 10) state and allows Rabi oscillation on the red 
sideband transitions 1 01) — » |10). Adapted from Chiorescu et 
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FIG. 53: Generation and control of entangled states and 
study of decay and lifetimes. Here a Rabi tt pulse on the os- 
cillator transition |00) — > 1 1 ) populates the state 1 01) , which 
then allows Rabi oscillation on the red sideband transition 
01) — » 1 10). The decay of both the Rabi oscillation and the 
average probability gives evidence for short oscillator life time 
(~ 3 ns; Q=100-150). Adapted from Chiorescu et alt— 



sitions. 




FIG. 54: Yale SCB charge-phase qubit coupled to an oscil- 
lator in the form of a superconducting microwave stripline 
resonator. Adapted from Wallraff et alw— 



C. Yale charge-phase qubit coupled to a strip-line 
resonator 

In the Yale experimentsi^iii^ a coherent qubit- 
quantum oscillator system is realized in the form of a 
Cooper pair box capacitively coupled to a superconduct- 
ing stripline resonator (Fig. I54|l forming one section of a 
microwave transmission line. The stripline resonator, a 
finite length (24 mm) of planar wave guide, is a solid-state 
analogue of the cavity electromagnetic resonator used in 
quantum optics to study strong atom-photon interaction 
and entanglement. In the Yale experiments, the qubit is 
placed in the transverse field in the gap between the res- 
onator strip lines, i.e. inside the microwave cavity. The 
Ec and Ej parameters are such that the SCT is effec- 
tively in the charge-phase region, and it is operated by 
controlling both the charge and phase (flux) ports. 

The large effective electric dipole moment d of Cooper 
pair box and the large vacuum electric field Eq in the 
transmission line cavity lead to a large vacuum Rabi fre- 
quency VRabi = 2dEo/h, which allows reaching the strong 
coupling limit of cavity QED in the circuit. 

In the Yale experiment, spectroscopic measurements 
are performed by driving the qubit with resonant mi- 
crowave pulses and simultaneously detecting the fre- 
quency and intensity-dependent amplitude and phase of 
probe pulses of microwave radiation sent through a trans- 
mission line coupled to the stripline resonator via in- 
put/out capacitors (Fig. 154^1 . The oscillator frequency 
is fixed at hv r 6 GHz. The qubit level separation, on 
the other hand, is tunable over a wide frequency range 
around 6 GHz in two independent ways: (i) by vary- 
ing the magnetic flux (f> e through the loop, forming the 
phase (flux) gat e l V Dl or (ii) by varying the charge n g via 
the voltage gate, primarily around the charge degeneracy 
point, n g —l/2, to minimize the effect of charge fluctua- 
tions. 

The flux bias is used to tune the qubit transition fre- 
quency at n g =l/2 to values larger or smaller than the res- 
onator frequency. The tuning the qubit frequency with 
the charge gate will provide two distinct cases: the qubit 
always detuned from the oscillator, and the qubit being 
degenerate with the resonator at certain values of n g . 



A central resullii^i is the evidence for the qui 
oscillator hybridization and splitting of the degener; 
1 01) and 1 10) states, z/oi — * v ±, as illustrated in F 
021 (right). This splitting is often called "vacuum Ral 
splitting because the ocillator is in its ground (vacuu 
state. 

The 1 01) , 1 10) level splitting is observed through i 
crowave excitation of the — > 1 resonance transit: 
and observing the splitting of the resonance line as 1 
qubit and the oscillator are tuned into resonance (by ti 
ing the qubit frequency to zero qubit-oscillator detunii 

S = V q — V r = 0). 

Another central resullii^ is the evidence for a lc 
qubit dephasing time of T 2 > 200 ns under optimal c< 
ditions: qubit parked at the charge degeneracy poi 
and weak driving field (low photon occupation numl 
in the resonator cavity). In the experimenfclSl the qut 
oscillator detuning is large, the qubit resonance trar 
tion v q — voi is scanned by microwave radiation, and 1 
dispersive shift of the resonator frequency v r seen by a 
microwave probe beam is used to determine the occupa- 
tion of the qubit levels. Scanning the qubit |00) — + 1 10) 
resonance line profile allowed to determine the lineshape 
and linewith as a function of microwave power, giving the 
best value of T 2 > 200 ns. Moreover, observation of the 
postion of the resonance as a function of microwave power 
allowed the determination of the ac-Stark shift, i.e. the 
Rabi frequency as a function of the photon occupation 
(electric field) of the cavity. The measurement induces 
an ac-Stark shift of 0.6 MHz per photon in the qubit level 
separation. Fluctuations in the photon number (shot 
noise) induce level fluctuations in the qubit leading to 
dephasing which is the characteristic back-action of the 
measurement. A cross-over from Lorentzian to Gaus- 
sian resonance line shape with increasing measurement 
power is observed and explained by dependence of the 
resonance linewidth on the cavity occupation number, 
exceeding the linewidth due to intrinsic decoherence at 
high rf-power. 



D. Comparison of the Delft and Yale approaches 

In comparison, the Delft experiment^ corresponds to 
the case of very large detuning (v q = 6 GHz, v r = 3 GHz; 
6 = Vg — v r ~ v r ; Fig. 1491 (left)), so that one will observe 
a main resonance line 1 00) — > 1 10) and two sidebands, 
1 00) -> 1 11) (blue) and |01) -> 1 10) (red). Decreasing the 
detuning (5, the blue and red sidebands will move away 
to higher resp. lower frequencies, and one will arrive at 
the case of zero detuning and qubit-oscillator degeneracy. 
If the qubit-oscillator coupling is larger than the average 
linewidth, one will then observe a splitting of the main 
line (Fig. |H (right)). 
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FIG. 55: The NEC 2-qubit system: two capacitively cou- 
pled charge qubits. The left qubit is a single Cooper pair 
box (SCB) and the right qubit is an SCT with flux-tunable 
Josephson energy. Courtesy of J.S. Tsai, NEC, Tsukuba. 
Japan. 



XIII. EXPERIMENTAL MANIPULATION OF 
COUPLED TWO-QUBIT SYSTEMS 

A. Capacitively coupled charge qubits 

An AFM picture of the NEC SCB 2-qubit system^ is 
shown in Fig. 1551 and the corresponding circuit J J circuit 
in Fig. EHl (note that in the NEC circuit the left SCT is 
replaced by a simple SCB). 

Two coupled qubits constitute a 4-level system. The 
Hamiltonian for the NEC system of two capacitively cou- 
pled charge qubits (SCBs) was analyzed in the Appendix. 
The four energy eigenvalues ^1^,3,4(^1,^2) are plot- 
ted in Fig. [23 as a functions of the gate charges. At 
each point in the parameter space, an arbitrary two-qubit 
state can be written as a superposition of the four two- 
qubit energy eigenstates. 

As described in Sections IV and X, the general proce- 
dure for operating with dc-pulses is to initialize the sys- 
tem to its ground state 1 00) at the chosen starting point 
(n ff io,n 9 2o) and then suddenly change the Hamiltonian 
at time t = to the gate bias (n g \, n g2 ). If the change is 
strictly sudden, then the state at (n g i, n g2 ) at time t = 
is \ip(t = 0)) = 1 00), which can be expanded in the energy 
basis of the Hamiltonian, 

|V(0)) = |00) =c 1 \E 1 ) +c 2 \E 2 ) + C3 \E 3 )+c i \E 4 ) (13.1) 

This stationary state then develops in time governed by 
the constant Hamiltonian as 

|^(i)) =cie- iBl *|i? 1 )+c 2e -^ t |£ 2 ) + 
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(i) Input preparation 



FIG. 56: The NEC 2-qubit system: Energy-level structure 
as a function of the gate charges n g i and n g i, independently 
controlled by the gate voltages V g i and V g 2- Courtesy of J.S. 
Tsai, NEC, Tsukuba, Japan. 



+c 3 e- iE3t \E 3 ) + Ci e- iEit \E 4 ) (13.2) 

If one re-expands this state in the charge basis of the 
starting point A (Fig. I56|l . then one obtains a system 
with periodic coefficients a, (t) , 

\i/>(t)) = oi(t)|0Q) + a 2 (4)|01) + a 8 (t)|10) + a 4 (4)|ll) 

(13.3) 

developing in time through all of the charge states. 

To perform a two-qubit conditional gate operation, one 
performs a series of pulses moving the system around 
in parameter space. Specifically, the NEC scheme is to 
apply sequential dc-pulses to the charge pulse gates of 
each of the two qubits in Fig. [S3] Two cases have been 
studied so far: 

(1) V g i(t) = V g 2(t): This is the case with common 
control dc-pulses for studied by Pashkin et al« Since 
ei(t) = £2(4), the plane of operation is at 45 degrees to 
the axes in Fig. |3S1 As a result, Pashkin et al^ observed 
interference effects and beating oscillations between the 
two qubits, as described by Eq. I|10.43|) . 

(2) V gl {t) = 0;V g2 {t), V gl (t);V g2 (t) = 0: This is the 
case with separate sequential dc pulses on separate gates 
recently studied by Yamamoto et alw The scheme is 
illustrated in Fig. [SBI Starting at point A, first the sys- 
tem is pulsed in the n g i direction, putting the system 
in a superposition of the states at points A and B; then 
the system is pulsed in the n g2 direction, allowing condi- 
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FIG. 57: The NEC 2-qubit system: pulse sequences for a 
CNOT gate (actually, the gate is a NOT - CNOT gate). 
Courtesy of J.S. Tsai, NEC, Tsukuba, Japan. 



tional gate operation due to the different time evolution 
of the states departing from A or B. 

Specifically, an entangling gate of controlled-NOT 
(CNOT) type was demonstrated by Yamamoto et aliSi 
using the pulse sequences shown in Fig. I57I First one 
applies a dc-pulse to gate 1 of the control qubit, mov- 
ing out (down left) left from |00) (point A in Fig. I56|) 
in the n g \ direction to the single-qubit degeneracy point 
(point E), staying for a certain time, and then moving 
back (turning off the pulse) , putting the system in a su- 
perposition a|00) + /3| 10) (points A and B). Next one 
applies a dc-pulse to gate 2 of the target qubit, moving 
out (right) in the n g2 direction and back. The pulse is 
designed to reach the first degeneracy point (point G), 
allowing the state to develop into a superposition of 00 
(point A) and 01 (point C) if the control state was 00. If 
instead the control state was 10 (point B), the dc-pulse 
on gate 2 will not reach the two-qubit degeneracy point 
(point H) and the development will be roughly adiabatic, 
taking the state back to 10 (point B), never reaching 11 
(point D). 

With timing such (7r/2-pulse) that 00 -> 01 (A -> G -> 
C), the control gate leads to 00+10, and the target gate 
only modifies the first component (00 — > 01), resulting in 
01+10, i.e. one of the Bell states. This is shown in the 
top panel of Fig. |CT| 

If instead one first applies a 7r/2-pulse in the n g i di- 
rection to the target qubit, inducing 00 — * 01 (C), and 
then applies a 7r/2-pulse in the n g \ direction to this state, 
reaching 01+11, and then again applies a 7r/2-pulse to the 
target qubit 2, resulting 00+11, one obtains the other 
Bell state. This is shown in the bottom panel of Fig. |CT| 
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Flux bias and readout pulse 



Microwave-pulse qubit operation 





FIG. 58: Scanning microscope image of two inductively cou- 
pled qubits surrounded by a DC-SQUID. Courtesy of J. Ma- 
jer, TU Delft. 



B. Inductively coupled flux qubits 

Figure EH shows a circuit with two inductively coupled 
flux qubits forming a four-level quantum system, excited 
by a single microwave line and surrounded by a single 
measurement SQUI D 37 ' 38 . With this circuit Majer et 
ali 37 ' 38 have spectroscopically mapped large portions of 
the level structure and determined the device parameters 
entering in the Hamiltonian matrix Eq. (|10.10|l . finding 
good agreement with the design parameters. Majer et al. 
found clear manifestations of qubit-qubit interaction and 
hybridization in the level structure as a function of bias 
flux. Presently ter Haar et a\M^ are investigating a more 
strongly coupled system with a J J in the common leg (cf. 
Fig. 132(1 . inducing Rabi oscillations and performing con- 
ditional spectroscopy along the lines described in Section 
XII in connection with the coupled qubit-oscillator sys- 
tem. 

A similar device with two flux qubits inside a cou- 
pling/measurement SQUID has recently been investi- 
gated by et aliiiSiii 7 ", so far demonstrating Rabi oscilla- 
tion of individual qubits. Finally, Izmalkov et ali2& have 
spectroscopically demonstrated effects of qubit-qubit in- 
teraction and hybridization for two inductively coupled 
flux qubits inside a low- frequency tank circuit. 



C. Two capacitively coupled J J qubits 

Capacitive coupling of two JJ qubits (Section IX 
E, Fig. I30[) has recently been investigated by several 
groups^SiiZ 8 ", showing indirect^ and direcliiZ 8 " evidence 
for qubit entanglement. Figure ISH1 shows the circuit used 
by McDermott et al^Z&, The potential-well and level 
structure of each J J qubit under operation and measure- 
ment conditions are shown in Fig. 1601 

The 2-qubit circuit behaves as two dipole-coupled 
pseudo-2-level systems, illustrated in Fig^Hl With "iden- 
tical" qubits, the 1 01) and 1 10) states are degenerate and 



FIG. 59: Circuit scheme for two capacitively coupled current 
(flux) biased JJ qubits with rf and dc control/readout lines. 
Adapted from McDermott et alw^ 




FIG. 60: Potential energies and quantized energy levels the 
Josephson phase qubit: left, during qubit operation; right, 
during state measurement, in which case the qubit well is 
much shallower and state |1) rapidly tunnels to the right hand 
well. Adapted from McDermott et aim— 

become hybridized and split by the interaction. A mi- 
crowave 7r-pulse tuned to the |0) — ► |1) transition of one of 
the isolated qubits will "suddenly" induce a 1 00) — > 1 10) 
transition, populating the \E+) = 1 10) + 1 1) and \E—) = 
1 10) — 1 01) states with equal weights. This will leave 
the system oscillating between the two qubits, between 
the 1 10) and |01) states, \ip(t)) = \E+) + e~ lSEt \E_) = 
cos(£Et/2)|10) + sin(£Ei/2)|01) where SE = E- - E+. 
This means that the two-qubit system oscillates between 
non-entangled and entangled states. With one ideal de- 
tector for each qubit we can measure the state of each 
qubit with any prescribed time delay between measure- 
ments. With simultaneous measurements we can deter- 
mine all the probabilities p^ = \(ij\ip(t))\ 2 , ideally giv- 
ing p w = |(1 + cos(SEt)), poi = |(1 — cos(SEt)), and 
Poo = Pn = in the absence of relaxation, decohcrence 
and perturbations caused by the detectors. 

Figure HH] shows the actual experimental results of 
Ref^i 8 " for the probabilities pij (t): The pio and poi prob- 
abilities oscillate out-of-phase, as expected. In addition, 
the average probabilities and oscillation amplitudes all 
decay with time. The results are compatibl o 178 with the 
the finite rise time of the initial 7r-pulse (5 ns), the single- 
qubit readout fidelity (70 percent) and the single-qubit 
relaxation time T\ (25 ns), and the limited two-qubit 
readout fidelity as used in the numerical simulations. Of 
particular interest is that for simultaneous (within 2 ns) 
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FIG. 61: Interaction of the coupled qubits in the time do- 
main. The qubits are tuned into resonance and the system 
is suddenly prepared in state (10) by a microwave 7r-pulse 
applied to qubit 1. Simultaneous single-shot measurement 
of each of the qubits 1 and 2 reads out the probabilities 
poo,pio,poi,JJii for finding the 1+2 system in states 00, 10, 
01, and 11, respectively. The full lines represent results of 
numerical simulations. Adapted from McDermott et aim— 



readout of the qubits, the experiments show that readout 
of one qubit only leads to small perturbation of the other 
qubrfci2£, which is promising for multi-qubit applications. 



XIV. QUANTUM STATE ENGINEERING WITH 
MULTI-QUBIT JJ SYSTEMS 



Due to the recent experimental progress, protocols 
and algorithms for multi-qubit JJ systems can soon 
begin to be implemented in order to test the perfor- 
mance of JJ qubit and readout circuits and to study 
the full dynamics of the quantum systems. The gen- 
eral principles are well known (see e.g4 9 ' 50 ) and have 
very recently been successfully applied in several other 
systems to achieve interesting and significant results 
in well-controlled quantum systems: ion-trap technol- 
ogy has been used to entangle 4 qubitsi^i, implement 
2-qubit gates and test Bell's inequalitie9iSi*i24, per- 
form the Deutsch-Josza algorithmic, achieve telepor- 
tation (within the system)ASii2£ii2£, and perform er- 
ror correctioniS; quantum optics has recently demon- 
strated long-distance teleportationiSS as well as 4-photon 
cntanglementi22ii2i. There are presently a considerable 
theoretical literature on how to implement these and sim- 
ilar protocols and algorithms in J J circits. Below we will 
describe a few examples to illustrate what may need to 
be done experimentally. 
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FIG. 62: EPR anticorrelation in the singlet state 



A. Bell measurements 

The first essential step is to study the quantum dy- 
namics of a two-qubit circuit and to perform a test of 
Bell's inequalities by creating entangled two-qubit Bell 
states (Section IIIII Refi^ 9 .) and performing simultaneous 
projective measurements on the two qubits, similarly to 
the ion trap experimentsi^iiiSi. Clearly the experiment 
will be a test of the quantum properties of the circuit 
and the measurement process rather than a test of a Bell 
inequality. 

The general principle is to (a) entangle the two qubits, 
(b) measure the projection along different detector axes 
(" polarization directions" ) , (c) perform four independent 
measurements, and finally (d) analyze the correlations. If 
the detector axes are fixed, as is the case for the JJ read- 
out (measuring charge or flux), one can instead rotate 
each of the qubits. 

Figure 1621 shows the quantum network for creating a 
Bell state perform single qubit rotations R x (<f>i), 
Rx{.^2)i and finally perform a projective measurement on 
each of the qubits along the same common fixed quan- 
tization axis. For each setting of the angles, (0i,02), 
on then performs a series of measurements obtaining the 
probabilities Py = \(ij\4>)\ 2 = (V'luXulV')- These can 
be combined into the results for finding the two qubits 
in the same state, P sam e = Poo + Pii- or m different 
states, Pdiff = Poi + -PtOj an( l finally into the difference 
q(<fii, <h) = Psame- Pdiff - This quantity q(<p 1 ,<p 2 ) is eval- 
uated in four experiments for two independent settings 
of the two angles, 4>i — (a,<5), (f>2 — (A 7), constructing 
the function 



\q(a,(3) 
k(a,7) 



q(S,0)\ 
•?(<*, 7) I 



(14.1) 



The Clauser-Horne-Shimony-Holt (CSCH) 

conditioniS for violation of classical physics is then 
B > 2 (maximum value 2^2). 

The application to JJ charge-phase qubits has been 
discussed by Refsi 193 i 194 . Experimentally, a first step 
in this direction has been taken by Martinis et al»i££ 
who directly detected the anticorrelation in the oscillat- 
ing superposition of Bell states \tp{t)) — cos(S Et /2)\10) + 
sm{6Et/2)\01). 



49 



lTp>=alO>+bl1> 



I0> [ 



I0> 



Classical channel 

/ 

» k • 



if* 



FIG. 63: Teleportation in a 3-qubit system. The unknown 
qubit (1) to be teleported is \tp) — a\0) + b\l). The result 
of the measurement is sent by Alice via classical channels to 
Bob who performs the appropriate unitary transformations to 
restore the original single-qubit state. As a result of the tele- 
portation, a specific but unknown state has been transferred 
from one qubit to another. 



B. Teleportation 

In the simplest form of teleportation an unknown 
single-qubit quantum state is transferred from one part 
of the system to another, i.e. from one qubit to another, 
as illustrated in Fig. ESI In Fig- EUD the initial state is 
given by 



( O |0) + 6|l»(|00)) 



(14.2) 



Next, applying CNOT and Hadamard gates entangles 
qubits 2 and 3 into a Bell state, 



( O |0)+6|l))(l/*2)(|00) + |ll)) 



(14.3) 



which is the resource needed for teleportation (in quan- 
tum optics this corresponds to the entangled photon pair 
shared between Alice and Bob) . One member of the Bell 
pair (qubit 2) is now "sent" to Alice and entangled with 
the unknown qubit to be teleported, creating a 3-qubit 
entangled state, 

|00)(o|0)+6|l)) (14.4) 
+ |01>(6|0)+a|l» 
+ |10)(a|0) - b\l)) 
+ |ll)(-6|0)+o|l» 

At this point, qubit 3 is sent to Bob, meaning that a 3- 
qubit entangled state is shared between Alice and Bob. 
Moreover, at this point, in each component of this 3-qubit 
entangled state in Eq. I|14.4f) . the two upper qubits are in 
eigenstates. This means that a projective measurement 
of the these two qubits by Alice will collapse the total 
state to one of the four components. If the result of 
Alice's measurement is (ij), the first two qubits are in 
the state and the 3-qubit state is known, given by 
the corresponding component in the previous equation. 
Any of these 3-qubit products can be transformed by a 



|ip>= al0>+bl1> 



(IO>+l1>)/2 




FIG. 64: Teleportation in a 3-qubit system without measure- 
ment and classical transmission. The unknown qubit (1) to be 
transferred is again \ip) — a|0)+6|l). The state of qubits 1 and 
2 are now used to control CNOT gates to restore the original 
single-qubit state on qubit 3. As a result of the teleportation, 
a specific but unknown state has been transferred from qubit 
1 to qubit 3, leaving qubits 1 and 2 in superposition states. 



unitary transformation back to the original state: 

l(a\0)+b\l)) = \i>) (14.5) 

«7 B (6|0)+a|l)) = |^) (14.6) 

er,(o|0>-&|l» = M (14.7) 

a z a x (-b\0) + a\l)) = (14.8) 

corresponding to resp. no change, bit flip, phase flip 
and bit-plus-phase flip of the original unknown state 
Alice's measurement causes instantaneous collapse, after 
which Alice by classical means (e.g. e-mail!) can tell 
Bob which unitary transformation to apply to recover 
the original single-qubit state. 

An alternative approach, without measurement and 
classical transmission, is shown in Fig. H2] In this case, 
the final 3-qubit state is still a disentangled product state 
with the correct state \ip) = a\0) + b\l) on qubit 3, 
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(|o> 



|l»(a|0> + &|l» (14.9) 



while qubits 1 and 2 are now in superposition states. If 
desired, these states can be rotated back to |00) by single 
qubit gates. 

The teleportation protocol has been implemented in 



ion trap experiment: 



,186.187 



(and of course in quantum 



opticsiSi). A proposal for a setup for implementing a 
teleportation scheme in JJ circuitry is shown in Fig. 1651 
describing a 3-qubit chain of charge-phase qubits coupled 
by current-biased large JJ oscillators, in principle allow- 
ing controllable nearest-neighbour qubit couplings iSiiS: 
In the absence of bias currents, to first order the qubits 
are non-interacting and isolated from each other. The 
basic two-qubit gates are achieved by switching on the 
bias-current-controllcd qubit-qubit interaction for a given 
time. CNOT gates can be achieved in combination with 
single-qubit Hadamard and phase gate o 158 i 159 i 160 . By ap- 
plying the sequence of gates shown in Fig. 1641 (or an 
equivalent sequence), the unknown state will be physi- 
cally moved from the left end to the right end of the chain. 
Moreover, by applying coupling pulses simultaneously to 
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FIG. 65: Network of loop-shaped SCT charge qubits, coupled 
by large Josephson junctions. The interaction of the qubits 
(i) and (i + 1) is controlled by the current bias 1^ Individ- 
ual qubits are controlled by voltage gates, V g i. Single-qubit 
readout is performed by applying an ac current pulse to a par- 
ticular J J readout junction (not shown), or using an RF-SET 
capacitively coupled to the island [27]. 



several qubits, one can in principle create multi-qubit 
entangled states in fewer operations than with sequential 
two-qubit gate»i2&ii2& as well as perform operations in 
parallel on different parts of the system. 

Extending the system to a five-qubit chain one can in 
a similar way transfer an entangled two-qubit state from 
the left to the right end of the chain. 



C. Qubit buses and entanglement transfer 

With entangled "flying qubits" like photons, quan- 
tum correlations can be shared in a spatially highly ex- 
tended states. However, with solid-state circuitry the 
issue becomes how to transfer entanglement among spa- 
tially fixed qubits. 

The standard answer is to apply entangling two-qubit 
gates between distant qubits. The qubit-qubit interac- 
tion can be direct, e.g. dipole-dipole-type interaction be- 
tween distant qubits, or mediated by excitations in the 
system. The transfer can be mediated via virtual or real 
excitations in a passive polarizable medium, a "system 
bus" , or via a protocol for coupling qubits and bus oscil- 
lators. 

A classic example of protocol-controlled oscillator- 
mediated coupling is the Cirac-Zoller gatei^I between 
two ions sequentially entangled via exchange coupling 
with the lowest vibrational mode of the ions in the trap. 
Another example is the Molmer-Sorensen g a tei22ii22i222i 
which creates qubit-qubit interaction via virtual ex- 
citation of ion-trap modes. Similar concepts for 
controlling entanglement have recently been theoreti- 
cally investigated in applications to JJ-qubit-oscillator 

^^.[. ,,164. 201. 202. 203. 204. 205. 206. 207. 208 

An different approach is to allow the bus (" spin-chain" ) 
states to develop in time governed by the fixed bus Hamil- 
tonian and to tailor the interactions and the initital con- 
ditions such that useful dynamics and information trans- 
fer is achieved222i2i2i2iii2i2. Also these concepts have 
been applied to JJ-circuits in a number of recent the- 



FIG. 66: Teleportation with error correction in a 5-qubit 
system without measurement. 
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FIG. 67: Coding, decoding and correcting a bit flip error in 
a 3-qubit logic qubit. The first gate on the left represents two 
sequential CNOT gates. The last gate befor the garbage can 
is a control-control-NOT (Toffoli) gate. 
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D. Qubit encoding and quantum error correction 

Quantum error correction (QEC) (see e.g. 
Refs, 49 ! 220 ! 221 ! 222 ! 223 ! 224 ! 225 ! 226 ) will most certainly 
be necessary for successful operation of solid state quan- 
tum information devices in order to fight decoherence. 
The algorithmic approach to QEC follows the principles 
of classical error correction, by encoding bits to create 
redundancy, and by devising procedures for identifying 
and correcting the errors based on specific models for 
the errors. 

The quantum circuit in Fig. 1661 illustrates QEC in 
terms of five-qubit teleportation with bit errors, illus- 
trating restoration of the state including error correction 
by measurement-controlled unitary transformation. 

The quantum circuit in Fig. 1671 demonstrates some 
essential steps of QEC in the case of one logical qubit 
encoded in three physical qubits: The first step is to 
encode the physical qubit in logical qubit basis states 
1 000) and | XIX) by applying two CNOT gates to create a 
3-qubit entangled state: 

(a\0) + b\l)) 1 00) -> (a|000) + 6jXXX» (14.10) 

Next, there may occur a bit-flip error in one of the qubits: 

a|000) + £>|1H) (no bit flip) (14.11) 
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FIG. 68: Coding, decoding and correcting a bit flip error in 
3-qubit logic qubit. 



a|100) + 6|011) 
a|010> + 6|101> 
o|001) + 6|110) 



(bit flip in qubit 1) (14.12) 
(bit flip in qubit 2) (14.13) 
(bit flip in qubit 3) (14.14) 



The next step applies a number of disentangling CNOT 
gates to check for the type of error. For the four possible 
states above we obtain: 



o|000) + 6|111> 
a|100) +6|011) 
a|000) +6|111) 
o|000> + 6|111) 



(o|0)+6|l))|00) (14.15) 

(a|l)+6|0» |11> (14.16) 

(o|0)+6|l»|10) (14.17) 

(o|0)+6|l))|01) (14.18) 



At this stage, qubits 2 and 3 have become independent 
eigenstates, just as in tcleportation. The correspond- 
ing eigenvalues are called syndrome, and indicate which 
corrective operations should be implemented. Remark- 
ably, in two of the above bit-flipped cases, with error 
syndromes ((01) and (10), qubit 1 is now correct, the er- 
ror residing in qubits 2 or 3 in the workspace (" ancillas" ). 
The only case where a transformation is needed is when 
the error syndrom is (11), which corresponds to a bit 
flip in qubit 1, to be corrected with a CCNOT (control- 
control-NOT, or Toffoli, gate), controlled by the truth 
table of an AND gate. In Fig. [ft7|we have used a com- 
pact single-gate notation for CCNOT, while in reality it 
must be implemented through a sequence of eight two- 
qubit gates 4 ^ (there is no direct three-qubit interaction 
in the Hamiltonian). At this stage, the 3-qubit state is 
completely disentangled into a product state. 

The final step consists in re-encoding the physical qubit 
1 to the logical qubit a|000) +6|111). However, although 
the physical qubit 1 is always correct at this stage, the 
workspace is not. This can be handled in a number of 
ways. Figure dumps the "hot" qubits in the garbage 
can (i.e. leaves them, or forces them, to relax), and 
re-encodes with fresh qubits from a larger workspace. 
Alternatively, Fig. HJB1 describes a variation where the old 
qubits 2 and 3 are re-initialized by entanglement with a 
measurement device which then dissipates the heat from 
the bitflip and the error correction procedure. 

In both of these cases, as described, we need a total 
workspace with five qubits. In principle, however, one 
can do with three qubits if we can rapidly re-initialize 



qubits 2 and 3 without inducing new errors. Recently 
an error correcting procedure with cooling of the system 
without measurement has been proposed 2 ^. 

Phase flips (sign changes), e.g. a|0) + 6|l) — >a|0) — 
can be handled by similar 3-qubit encoding, decoding and 
correction. Combining these two approaches gives the 9- 
qubit Shor cod o 220 for correcting also combined bit and 
phase flips, e.g. a|0) + 6|1) — > a| 1) — &|0). The minimum 
code to achieve the same thing is a 5-qubit code^ 21 . and 
there are more efficient codes with 7 qubits 2 ^ 2 ^ 2 ^. 

A related approach to fighting decoherence is to en- 
code the quantum information in noiseless subsystems, 
so called Decoherence Free Subspaces (DFS) 124 i 12 M 28 i 229 

A different approach is the so called "bang-bang" 
dynamic correction method 230 i 231 i 232 i 233 i 234 i 235 , basically 
corresponding to very frequent application of the spin- 
echo technique. This is related to the quantum Zeno 
effect (see e.g. Refi 2 ^, describing situations where the 
quantum system via very strong interaction with the en- 
vironment is forced into a subspace from where it cannot 
decay. 

In this Section we have restricted the discussion to 
a few different protocols for quantum state engineering, 
representing basic steps in algorithms for solving specific 
computational problems. For a discussion of quantum 
algorithms and computational complexity we refer to the 
books by Nielsen and Chuang 4 ^ and Gruska-S, and to 
the original papers. We would however like to mention a 
few papers discussing how to implement a few well-know 
algorithms in JJ circuitry. The basics of quantum gates 
in JJ-circuits can be found in e.g. Refs. 42 - 237 - 238 . The 
Deutsch-Josza (DJ) and related algorithms have been 
discussed by Siewert and FazirjS 3 ^ describing a 3-qubit 
DJ implementation with three charge-phase qubits con- 
nected in a ring via phase-coupling JJs with variable 
Josephson coupling energy (SQUIDS), and by Schuch 
and Siewert analyzing a 4-qubit implementation 2 ^. For 
the Grover search algorithm there seems to be nothing 
published on the implementation in JJ circuitry. Con- 
cerning Shor's factorization algorithm there is a recent 
paper by Vartiainen et al^a^. There is also recent work 
on optimization of two-qubit gates 2 ^ 3 -, and relations be- 
tween error correction and entanglements 2 ^ 

Finally there are a number of papers connecting JJ- 
ciruits wih geometric phases and holononomic quantum 



computing 



246,247.248.249 



and on adiabatic computatio 



252 



and Cooper pair pumps2£2i25i. For a recent paper dis- 
cussing the universality of adiabatic quantum computing, 
see Aharonov et ali^S. 



XV. CONCLUSION AND PERSPECTIVES 

Within 5 years, engineered J J quantum systems with 
5-10 qubits will most likely begin seriously to test the 
scalability of solid state QI processors. 

For this to happen, a few decisive initial steps and 
breakthroughs are needed and expected: The first essen- 
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tial step is to develop JJ-hardware with long coherence 
time to study the quantum dynamics of a two-qubit cir- 
cuit and to perform a "test" of Bell's inequalities (or 
rather the JJ-ciruitry) by creating entangled two-qubit 
Bell states and performing simultaneous projective mea- 
surements on the two qubits. 

A first breakthrough would be to perform a significant 
number of single- and two-qubit gates on a 3-qubit clus- 
ter to entangle three qubits. Combined with simultane- 
ous projective readout of individual qubits, not disturb- 
ing unmeasured qubits, this would form a basis for the 
first solid-state experiments with teleportation, quantum 
error correction (QEC), and elementary quantum algo- 
rithms. This will provide a platform for scaling up the 
system to 10 qubits. 

This may not look very impressive but nevertheless 
would be an achievement far beyond expectations only 
a decade back. The NMR successes, e.g. running Shor- 
type algorithms using a molecule with 7 qubits^, are 
based on technologies developed during 50 years using 
natural systems with naturally long coherence times. 
Similarly, semiconductor technologies have developed for 
50 years to reach today's scale and performance of clas- 
sical computers. It is therefore to be expected that QI 
technologies will need several decades to develop truly 
significant potential. Moreover, in the same way as for 
the classical technologies, QI technologies will most prob- 
ably develop slowly step by step, "qubit by qubit" , which 
in itself will be an exponential development. 

Moreover, in future scalable information processors, 
different physical realizations and technologies might be 
combined into hybrid systems to achieve fast processing 
in one system and long coherence and long-time informa- 
tion storage in another system. In this way, solid state 
technologies might be combined with ion trap physics to 
build large microtrap systems^ 5 ^, which in turn might be 
coupled to superconducting Josephson junctions proces- 
sors via microwave transmission lines 255 . 
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GLOSSARY 

Adiabatic evolution - Development of a quantum sys- 
tem without transitions among the quantum levels. 
Algorithm - Finite sequence of logical operations, which 
produces a solution for a given problem. 
Level crossing - Degeneracy of quantum levels appear- 
ing at a certain value of a controlling system parameter 
(e.g. gate voltage, bias flux, etc.). 



Anticrossing - Lifting of degeneracy (level crossing) of 
quantum levels during variation of the system parameters 
when an interaction is switched on 

Average measurement - Measurement of an expecta- 
tion value of a dynamic variable in a certain state 
Bloch sphere - Geometrical representation of the man- 
ifold of quantum states of a two-level system as points 
on the unit sphere. 

Bloch vector - Normalized state vector of a two-level 
system represented by a radial unit vector of the Bloch 
sphere. 

Charging energy - Electrostatic energy of a capacitor 
charged with a single electron (e) (or a single Cooper pair 
(2e)). 

Charge qubit - Superconducting qubit based on a a 
single Cooper pair box (SCB), whose computational basis 
consists of the two charge states of the superconducting 
island. 

Charge-phase qubit - Superconducting qubit based on 
a single Cooper pair box (SCB) whose charging energy is 
of the order of the Josephson energy of tunnel junctions 
CNOT gate - Controlled-NOT gate: two-qubit gate 
which changes or does not change the state of a target 
qubit depending on the state of a controlling qubit. 
Cooper pair - Bound state of two electrons (2e), the 
elementary charge carrier in superconducting equilibrium 
state. 

Coulomb blockade - Suppression of current through 
a tunnel junction or small metallic island due to large 
charging energy associated with a passage of a single elec- 
tron. 

CPHASE gate - Controlled-phase gate: two-qubit gate 
which changes or does not change the phase of a target 
qubit depending on the state of a controlling qubit. 
Decoherence - Evolution of a quantum system, inter- 
acting with its environment; cannot be described with 
a unitary operator; consists of decay of phase coherence 
(dephasing) and/or changing of level population (relax- 
ation). 

Density matrix - Characteristics of a quantum system, 
which contains full statistical information about the state 
of the system. 

Dephasing - Decay of phase coherence of a superpo- 
sition state, represented by decreasing off-diagonal ele- 
ments of the density matrix. 

Entanglement - Specific non-local coupling of quantum 
systems when the wave function of whole system cannot 
be presented as a product of partial wave functions 
Flux qubit - Superconducting qubit based on a SQUID, 
whose computational basis consists of the two states of 
the SQUID having opposite directions of the induced 
flux. 

Hadamard gate - Transformation of computational ba- 
sis states of a single qubit to equally weighted superpo- 
sitions of the basis states (cat states). 
Holonomic quantum computation - Using the geo- 
metric phases when a quantum system is taken around a 
closed circuit in the space of control parameters. 
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Gate operation - Controlled transformation of the state 
of one or several qubits; a basic element of an algorithm. 
Josephson effect - Non-dissipative current flow be- 
tween two superconductors separated by a non- 
superconducting material (insulator, normal metal, etc.). 
Josephson junction - Junction of two superconductors, 
which exhibits the Josephson effect. 
Josephson critical current - Maximal value of the 
Josephson current maintained by a particular junction. 
Josephson energy - Inductive energy of a Josephson 
junction proportional to the critical Josephson current. 
7r pulse - High frequency control pulse with a specific 
duration applied to a qubit, producing inversion of the 
qubit level populations (tt rotation; qubit flip) 
it/2 pulse - High frequency control pulse with a specific 
duration applied to a qubit, typically tipping the Bloch 
vector from a pole to the equator, or from the equator to 
a pole, on the Bloch sphere. 

Phase gate - Single qubit gate, transforms a superpo- 
sition of two quantum states into another superposition 
with different relative phase of the states 
Precession - Dynamic evolution of a two-level system in 
a superposition state, i.e. linear combination of energy 
eigenstates. 

QND measurement - Quantum Non-Demolition Mea- 
surement: measurement of a state of a quantum system, 
which docs not destroy the quantum state and makes 
possible repeated measurements of the same state. 
QPC - Quantum Point Contact: a constriction in a con- 
ductor with ballistic transport through a small number 
of conduction channels. 

Qubit - Quantum two-level system; basic element of a 
quantum processor. 

PCQ - Persistent Current Qubit: synonymous with flux 
qubit. 

Rabi oscillation - Dynamics of two-level system under 



resonant driving perturbation, consists of periodic oscil- 
lation of the level populations with the frequency propor- 
tional to the amplitude of the perturbation. 
Readout - Measurement of a qubit state. 
Relaxation - Change of population of the energy eigen- 
states resulting in approaching the equilibrium popula- 
tion. 

SCB - Single Cooper pair Box: superconducting analog 
of SEB, where it is energetically favorable to have only 
paired electrons on the island. 

SCT - Single Cooper pair Transistor: a superconducting 
device containing a small island whose charging energy is 
controlled by an electrostatic gate electrode to increase 
or decrease current flowing through the island from one 
large electrode (source) to another (drain). 
SEB - Single Electron Box: small metallic island con- 
nected to a large electrode via resistive tunnel junction, 
whose charging energy hence amount of trapped electrons 
is controlled by an electrostatic gate. 
SET - Single Electron Transistor: a device containing a 
small island whose charging energy is controlled by an 
electrostatic gate electrode to increase or decrease cur- 
rent flowing through the island from one large electrode 
(source) to another (drain). 

rf-SET - SET driven by an rf signal, is used as an ul- 
tra sensitive electrometer by monitoring a linear response 
function of the SET, which is sensitive to the electrostatic 
gate potential. 

SQUID - Superconducting Quantum Interferometer De- 
vice: a device consisting of a one or more Josephson junc- 
tions included in a superconducting loop. 
dc-SQUID - SQUID containing two Josephson junc- 
tions. 

rf-SQUID - SQUID containing one Josephson junction. 
Single-shot measurement - A measurement which 
gives an "up/down" answer in one single detection event. 
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